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SECTION  1 


INTRODUCTION 


The  Air  Transport  of  Radiation  (ATR)  Code  has 
been  developed  under  the  sponsorship  of  the  De¬ 
fense  Nuclear  Agency  (DNA)  as  its  standard 
model  of  initial  radiation  free  field  dose  (kerma 
free-in-air)  from  nuclear  weapons  in  uniform  aii 
and  air-over-ground  geometries.  The  ATR  code 
is  capable  of  predicting  free  field  dose  from 
x-ray,  prompt  gamma-ray,  neutron,  neutron- 
induced  (secondary)  gamma-ray  and  fission 
product  gamma-ray  components  of  initial  radi¬ 
ation.  It  makes  such  predictions  with  a  minimum 
of  user  input,  thus  making  it  ideal  for  use  by 
weapon  ejects  analysts  who  are  unfamiliar  with 
the  details  of  radiation  transport  methodology.  As 
such,  it  has  been  used  to  aid  in  the  establishment 
of  targeting  doctrine  and  troop  safety  exclusion 
radii.  Portions  of  the  code  have  been  extracted 
for  use  as  modules  in  larger  codes.  The  code  has 
also  been  used  to  supply  radiation  environment 
data  sets  for  use  in  a  variety  of  applications,  such 
as  systems  analysis  and  engagement  codes,  the 
ONA  Effects  Manual  One  (EM-1)  and  the  ONA 
Electronic  Handbook  Series. 

The  development  of  ATR  was  begun  in  1970  and 
continued  over  the  next  seven  years.  A  bibliogra¬ 
phy  of  ATR-related  reports  is  provided  in  Table 
1 .  In  its  initial  form  ATR  was  intended  to  provide 
transport  data  for  prompt  radiation  components 
for  sources  in  uniform  air  and  in  the  vicinity  of  the 
air-ground  interface.  By  the  time  version  4  was 
released  in  1976,  the  basic  code  had  been  ex¬ 
panded  to  include  delayed  gamma  rays  and  to 
treat  sources  near  the  top  of  the  atmosphere, 
among  other,  additional  capabilities. 

This  period  of  ATR  development  coincided  with 
major  improvements  in  the  state-of-the-art  in  ra¬ 
diation  transport  codes  and  cross  sections  and 
the  capabilities  of  computers  to  cop®  with  large 
computational  tasks.  Thus,  in  1 975  th* entire  ATR 
data  base  was  replaced  with  results  calculated 
using  new  DNA-sponsored  evaluations  of  nitro¬ 
gen  and  oxygen  cross  sections  (ATR3).  At  that 
time  new  corrections  to  the  uniform  air  data  for 
the  proximity  of  the  source  to  the  ground  were 
also  incorpora. ed.  These  were  derived  from 
DNA-sponsored  calculations  of  fission  and  14 
MeV  neutron  sources  at  heights  to  300  meters 
and  ranges  to  1500  meters,  such  problems  being 
the  largest  which  could  be  run  at  that  time. 


In  1977,  ATR  underwent  another  modification. 
The  effort,  sponsored  not  by  DNA  but  by  Ballistic 
Research  Laooratory  (BRL),  accomplished  the  in¬ 
stallation  of  source  energy-dependent  air-over¬ 
ground  correction  factors  for  neutron  and  gamma 
radiation.  The  data  from  which  those  correction 
factors  were  derived  were  calculated  under  DNA 
sponsorship  and  were  the  results  of  the  first  effort 
to  apply  DOT  in  the  adjoint  mode  for  such  a  prob¬ 
lem,  i.e.,  long  ranges  in  air-over-ground  geome¬ 
try.  That  version  of  the  ATR  code  is  known  as 
ATR4. 1 .  Subsequent  to  the  production  of  ATR4.  • 
it  was  found  that  the  two-dimensional  transport 
calculations  on  which  its  air-over-ground  correc¬ 
tion  factors  were  based  were  performed  using  a 
ground  moisture  level  significantly  oelow  normal 
levels.  This  caused  ATR4.1  to  overestimate  the 
neutron  dose  and  overestimate  the  secondary 
gamma  ray  dose  in  the  vicinity  of  the  ground.  As 
a  result  of  this  finding,  the  most  accurate  applica¬ 
tion  of  ATR,  using  version  4. 1  and  those  previous, 
may  be  obtained  by  using  ATR4  for  neutron  and 
secondary  gamma  ray  doie  data  and  ATR4.1  for 
prompt  and  delayed  gamma  ray  dose  data. 

The  above  paragraphs  provide  a  very  brief  sum¬ 
mary  of  the  complex  development  process  of  the 
ATR  code.  A  more  comprehensive  history  of  that 
development  up  to  1977  is  contained  in  Appendix 
A  of  this  report.  From  1977  to  the  present  ATR 
has  been  used  by  many  agencies  in  its  various 
versions  and  has  been  examined  tor  internal  con¬ 
sistency  and  tested  against  independent  data. 
This  report  describes  the  results  of  a  program  be¬ 
gun  in  1985  to  revise  ATR  based  on  insights 
gained  from  that  collective  experience.  The  pro¬ 
gram  includes  revision  or  replacement  of  the  fol¬ 
lowing  portions  of  ATR: 

o  Uniform  Air  Transport  Data  Base 
Neutron 
Gamma  Hay 
Secondary  Gamma  Ray 
Humidity  Correction 

o  Air-Ground  Interface  Transport  Correction 


Table  1.  ATR  documentation  summary. 


REPORT 

CONTENT 

SA I- 7 1-565- LJ 
November  1971 

Paper  given  at  RSIC  Workshop  on  Radiation  Transport 
in  Air.  Describes  logic  of  code  and  some  data  base 
development. 

DNA2803I 

Hay  1972 

First  report  describing  basic  concepts  of  ATR  and 
the  first  distributed  version  of  the  code;  describes 
neutron  and  secondary  gamma-ray  data  base  genera¬ 
tion.  Includes  air/ground  and  exponential  air 
correction  factors. 

DNA3144A 

April  1973 

Users  manual  for  ATR-2  version  of  the  code;  does 
not  describe  data  base  generation. 

DNA327ST 

August  1974 

Describes  data  base  generation  for  photons  (prompt 
gamma  rays  and  x-rays)  and  the  prompt  gamna-ray 
air/ground  correction  factors. 

DNA3362Z 

August  1974 

Summary  of  the  capabilities  of  the  ATR  code  with 
updates  to  ATR-2. 

DNA3819F 

July  1975 

Describes  ATR-3  including  new  data  base  using  DNA 
cross  section  library,  new  air/ground  correction 
factors,  low  energy  x-rays,  and  new  REGROUP  routine. 

DNA4061 

January  1976 

Describes  TDATR,  the  time-dependent  prompt  photon 
and  secondary  gamma-ray  version  of  ATR. 

DNA3395F 

January  1976 

Describes  fission  product  model  and  summarizes  total 
capability  of  the  ATR-4  code. 

BRL  CR  343 

August  1977 

Describes  ATR-4. 1  with  energy-dependent  air/ground 
correction  factors.  Work  supported  by  BRL. 
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o  Delayed  Gamma  Ray  Model 

Source 

Transport 

o  User  Features 

Air  density/ humidity  Input 

Fluence-to-dose  conversion  factors 

In  addition  the  revised  version  of  ATR.  which 
should  be  referred  to  as  ATR5.  has  been  exten¬ 
sively  tested  against  dose  measurements  made 
on  a  number  of  atmospheric  tests  carried  out  in 
Nevada  and  the  Pacific. 

1.1  UNIFORM  AIR  TRANSPORT  DATA. 

The  models  of  neutron  and  gamma-ray  transport 
in  uniform,  dry  air.  implemented  in  ATR,  through 
version  4.1,  are  based  on  calculations  using 
one-dimensional  discrete  ordina*es  methodology 
with  quadrature  and  multigroup  coupled 
cross  sections  composed  of  22  neutron  and  18 
gamma-ray  groups  with  a  P3  Legendre  polynomi- 
nal  scattering  approximation.  Such  calculational 
limitations  are  typical  of  those  necessitated  by 
available  computer  and  data  resources  at  the 
outset  of  ATR  development.  Now.  more  than  a 
decade  later,  modern  evaluated  cross  sections 
can  be  used  in  one-dimensional  calculations  hav¬ 
ing  an  S40  quadrature,  a  1 74  neutron  -  38  gamma 
ray  group  cross  section  set  and  a  Legendre  scat¬ 
tering  approximation  as  high  as  Pj.  Data  calcu¬ 
lated  using  such  methodology  are  coplotted  with 
those  from  ATR4/4.1  for  neutron,  secondary 
gamma  ray  and  prompt  gamma  ray  kerma  from 
a  fission  weapon  source  in  Figures  1.  2  and  3. 
respectively.  The  data  for  transport  through  dry 
air  vary  between  ATR4  and  newer  methods  differ 
by  twenty  percent  or  more,  depending  on  dis¬ 
tance.  More  importantly,  significant  variations 
occur  when  moisture  is  added  to  the  atmo¬ 
sphere,  particularly  for  neutron  kerma.  There¬ 
fore,  the  first  task  in  the  program  described 
herein  is  to  produce  a  new  uniform  air  data  set 
for  transport  from  neutron  and  gamma-ray 
sources  and  to  parameterize  those  data  for  inclu¬ 
sion  m  ATR.  The  second  task  is  to  develop  a  data 
set  of  radiation  transport  in  air  of  various  moisture 
contents  on  which  to  base  a  model  of  air  humidity 
effects  on  free  field  dose  and  to  implement  that 
model  in  ATR. 


1.2  AIR-GROUND  INTERFACE  TRANSPORT 
CORRECTION. 

ATR  versions  3  and  4  contain  models  for  the  per¬ 
turbation  of  spatial  dose  distributions  near  the 
air-ground  interface  based  on  two-dimensional 
discrete  ordinate  calculations  using  fission  and 
fusion  sources  having  heights  extending  to  300 
meters.  The  useful  horizontal  distance  limit  of 
those  data  is  1500  m.  The  simple  four-element 
ground  used  in  those  calculations  lacks  the  trace 
elements  necessary  to  accurately  depict  the  sec¬ 
ondary  gamma-ray  environment,  particularly  at 
short  distances,  as  shown  in  Figure  4. 

ATR4. 1  contains  an  air/ground  perturbation  mod¬ 
el  based  on  adjoint  two-dimensional  discrete  or¬ 
dinates  calculations,  which  provide  the  import¬ 
ance  of  each  source  particle  energy  group  (dose 
per  source  neutron  or  gamma  ray)  for  producing 
free  field  dose  one  meter  above  the  ground.  The 
useful  limits  of  these  calculations  are  approxi¬ 
mately  1000  meters  source  height  and  1200  me¬ 
ters  horizontal  distance.  These  calculations  were 
performed  for  a  mean  Western  European  soil 
containing  virtually  all  trace  elements  of  interest, 
but  with  an  abnormally  low  moisture  content.  The 
result  was  a  set  of  fully  source  energy-differential 
air-ground  correction  factors  'or  ATR4.1.  but 
ones  which  resulted  in  the  overprediction  of  neu¬ 
tron  dose  and  under  prediction  of  secondary 
gamma-ray  dose,  as  shewn  in  Figures  5  and  6. 
respectively. 

The  third  task  in  the  present  program  is  to  devel¬ 
op  and  implement  a  set  of  air-ground  interface 
dose  correction  factors  based  on  a  contempo¬ 
rary  adjoint,  two-dimensional  discrete  ordinates 
calculation  which  incorporates  a  240  angle  quad¬ 
rature.  a  37  neutron  •  21  gamma-ray  group  cross 
section  set  with  P3  Legendre  scattering,  ex¬ 
tended  spatial  limits  ( 1 000  m  vertical  2000  m  hor¬ 
izontal)  and  mean,  moist  Western  European  soil. 

1,3  DELAYED  GAMMA-RAY  TRANSPORT. 

The  delayed  (fission  product)  gamma-ray  model 
implemented  in  ATR4/4 ,1  is  an  appendage  having 
its  own  transport  data  set.  rather  than  being  fully 
integ.ated  into  ATR.  using  the  same  transport 
data  sets  as  the  prompt  radiation  component. 
Thus,  ATR4/4.1  does  not  display  a  total  initial 
gamma-ray  dose,  but  rather  a  total  comprised 
of  prompt  and  secondary  gamma  rays  and  a  sep¬ 
arate  tabulation  of  delayed  gamma-ray  dose, 
which  must  be  combined  by  hand  to  obtain  the 
true  total.  In  addition  to  this  inconvenience,  the 
dose  calculated  using  the  ATR4/4.1  delayed 
gamma-ray  model  does  not  agree  well  with  test 
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Secondary  gaaaa  ray  keraa  at  the  ground  fro*  a  boosted 
fission  source  (HOB  «  100  aeters)  calculated  using 
ATR5  (ooo)  and  ATR4  ( — ) . 


TISSUE  KERMA  PER  SOURCE  NEUTRON  ATR4.1/ATR5 


cU’ire  6.  Neutron  ker*»  fro*  a  fission  source  at  100  »eter» 
burst  heikht,  ATU4.1  and  ATT15. 


TISSUE  KERMA  PER  SOURCE  NEUTRON  ATR4.1/ATR8 


HORIZONTAL  DISTANCE  (m) 


Figure  8.  Secondary  gaaaa  ray  keraa  fro*  a  fission  source  at 
1QQ  aeters  burst  height,  ATR4.1  and  ATH5. 


measurements,  as  will  be  described  later. 

Task  four  of  the  present  effort  is  to  revise  tne 
delayed  radiation  model,  fully  integrating  it  into 
the  ATR  structure,  to  enabie  it  to  calculate  dose 
frcrr.  fission  products  of  the  three  rr'am  hssi'e  e»- 
ment«.  U-235,  IJ-238  and  Pu- 1 39.  and  to  test  -<■  o 
resulting  mode'  against  relevant  .veapon  test  data 
and  results  of  other  caiculationai  models. 

1.4  USER  FEATURES. 

In  addition  to  the  obvious  need  to  include  provi¬ 
sion  for  input  of  atmospheric  moisture  content 
and  weapon  fissile  components  associated  with 
each  weapon  yield  fraction,  as  required  by  the 
ATR  modifications  described  above,  a  fifth  task 
of  this  program  is  to  improve  input  commands 
and  other  user  features  of  the  cone,  such  as  the 
fluence-to-dose  conversion  factor  library.  Specif¬ 
ically,  the  input  commands  are  to  be  modified  to 
include  direct  specifications  of  mean  air  density 
between  source  and  detector  location 

1.5  TECHNICAL  APPROACH. 

The  aoproach  to  accomplishing  previous  revi¬ 
sions  of  ATR  was  to  maintain  a  comoiete  history 
of  the  evolution  of  ATR  imbedded  n  the  code  it¬ 
self.  i  e. .  to  make  the  code  capable  of  reproduc¬ 
ing  results  from  all  previous  versions  of  the  code 
as  well  as  that  from  the  current  version  tor  a  given 
problem.  That  was  achieved  either  by  overlaying 
correction  factors  upon  existing  models  or  by 
simoly  retaining  both  the  new  and  old  models  or 
data  sets  within  the  code.  The  overlay  or  additive 
modification  approach  was  used  primarily  be¬ 
cause  't  was  found  to  be  easier  to  parameterize 
corrections  to  the  existing  data  base  than  to  para¬ 
meterize  a  data  base  starting  from  scratch  Fur 
ther.  maintenance  of  previous  .ersiors  of  the 


code,  considered  necessary  for  the  sake  of  con¬ 
sistency  and  historical  perspective,  was  implicit 
m  tnct  approach. 

Aher  three  overlaid  revisions  to  ATR  it  was  de- 
c  fif'd  mat  the  previous  approach  Kad  reached  a 

cc-rt  of  d'mm.sh.ng  returns  tv"'-  n  'erms  of 
o  ~  v-i  ~.iar'ty  and  crecs'cn  ’“-e  p-esect  mortift- 
ca:  ons  to  ATR  are  based  on  caicu  at-onal  results 
derived  from  much  higher  quality  cross  section 
reoresentations  and  sophisticated  application  of 
radiation  transport  calculation  methods  than 
tncse  available  in  the  oast.  !n  addition,  parame¬ 
terization  of  the  new  data  can  be  accomplished 
usmg  approaches  not  tried  m  the  original  mudel. 
Therefora.  the  philosophy  of  retention  of  histori¬ 
cal  data  m  ATR  has  been  abandoned  m  favor  of 
simply  replacing  the  present  data  sets  with  re w 
ones.  The  command  structure  cl  ATR  has  been 
retained  with  a  minimum  o*  mod;''Caticns  in  order 
to  make  ATR  revision  as  transparent  to  the  user 
as  possible.  The  programming  style  has  been 
kept  to  a  level  common  to  most  computers,  re¬ 
taining  the  machine-independent  quality  of  the 
code  The  release  of  the  revised  version  of  the 
code,  known  as  ATR5.  ;s  being  held  in  abeyance 
pending  an  additional  modification  to  provide 
realistic  energy-  and  ang'e-rfiHerentiai  fluence  at 
the  air-ground  interlace  rnd  fmal  testing  Release 
'S  expected  during  calendar  year  1988. 

The  remainder  of  this  report  describes  (Section 
2)  me  suDerficiat  modifications  to  ATR.  including 
revised  input  commands,  energy  group  architec¬ 
ture.  internal  source  spectra,  and  tluence- 
to  dose  conversion  (actors.  (Section  3)  the 
revised  ATR  uniform  air  trarsoort  data  base, 
i Section  4)  the  revised  ATR  air  ground  interface 
perturbation  model,  and  i Section  5)  the  revised 
ATR  delayed  gamma-ra,  model,  including  tests 
of  mat  model  aga.nst  ■  nevant  weapon  test  data. 
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SECTION  2 


AIR  INPUT  FORMAT.  aND  DEFAULT  MODIFICATIONS 


s  section  describes  :Ke  major  sucer'c  ai 

C*"'3''-'^05  *0  A  *53  CCC0  3$  3  r©Su;t  Oi  ’^0 

Z'-'es©-"'’  cfograrr.  -''c*^ce 

c  ^put  Commands 

Meteorological  Data 

F:3Sfie  Material  Specification 

o  Energy  Group  Architecture 

o  internal  Source  Scectrum 

o  Fiuence-to-Oose  Converson  Factors 

2.1  METEOROLOGICAL  DATA  INPUT. 

specification  of  atmospheric  conditions  in 
has  been  modified  to  prcvde  for  direct  input 
of  tne  mean  density  between  source  and  detector 
and  ‘or  tno  input  of  atmosphere  moisture  coo- 
tsrt.  Previously,  the  air  dens  ty  used  in  ATP  cai- 
c»  a’ens  was  speeded  by  socatmg  tne  source 
ar  p  detector  withe  an  atmcscneric  density  pro- 
f-  e  cased  on  the  U  S.  Standard  Atmosphere  cry 
a  r i  That  was  acco  npr-sneo  bv  specifying  ground 
e  evat.on  relative  to  sea  eve)  and  source  and  de- 
•erdcr  height  relative  to  the  ground,  in  pr-nocs*. 


an  ATR  user  could  obtam  the  desired  average  air 
denscy  by  teat  and  error  plough  the  variation  of 
t"e  ground  e’evat-pn.  but  me  procedure  was 
bumperseme.  To  aiiev  a*e  m  s  problem  ana  to 

make  ATP  more  flexible  in  application,  two  modi¬ 
fications  have  been  mace  to  the  code.  First, 
while  the  standard  atmospheric  density  profile 
has  been  retained  withe  ATR.  allowing  atmo¬ 
spheric  conditions  to  be  compietety  specified  by 
the  input  of  source  height,  detector  height  and 
ground  elevation,  those  atmospheric  specifes- 
tions  can  be  overridden  through  the  use  of  a  new 
command  card.  ’MET.  Second,  the  specification 
of  atmospheric  conditions  thought  the  use  ot  the 
'MET  command  can  now  include  moisture  con¬ 
tent.  The  objective  of  this  modification  ip  to  per¬ 
mit  the  ATR  user  tospecify  a  set  of  meteorological 
parameters  to  obtain  a  mean  air  density  and 
mo  sture  content  which  exactly  corresponds  to  a 
problem  of  interest. 

The  format  of  the  input  command  is: 

'MET(i)  Va'uss 

where  the  index  l  determines  the  meaning  of  the 
sequence  of  values  following  the  parenthesis,  as 

fc'iows: 


Value 

Value 

Value 

1 

—2 — 

_2_ 

i 

T  Hirpfti  ature 

Pressure 

%Pe*ative  Humidity 

7 

Temperature 

Pressure 

Dew  Pomt  Temperature 

3 

Density  of  Mo  st  A.r 

Wt  s  of  wa’er 

Not  Used 

4 

0‘i'isity  of  MO'St  A.' 

Density  of  Water  Vapor 

Noi  Used 

s 

Density  of  Dry  Air 

Wt  N  of  Water 

Not  Used 

6 

Density  rf  Dry  A  if 

Density  of  Water  Vapor 

Not  Used 

t  i 


As  with  past  ATR  commands,  the  values  may  be 
entered  in  free  format  separated  only  by  a  space. 
Units  of  input  values  are  restricted  to: 


The  equations  ,nco<  porated  in  ATR  ‘or  the  calcu¬ 
lation  of  atmospheric  density  components  ar9 
provided  in  Appendix  B  of  this  repci: 


Temperature:  Decrees  Centigrade  t°C) 

pressure  Miil'bars  (mb) 

Relative  Hum.dity:  Percent  iG) 

Density:  Grams  per  cubic  centimeter  (g.cc) 

If  a  ’  MET  command  is  not  used.  ATR  automatical¬ 
ly  determines  the  average  air  density  between 
source  and  detector  locations  as  :n  previous  ver¬ 
sions  of  ATP.  Also,  the  moisture  content  of  the 
air  ts  set  to  a  dafajlt  va/ue  of  0  0  percent  oeni  fy 
(dry  air). 

Examples  of  ’MET  input  commands  are  a3  fol¬ 
lows: 

Example  1.  *MET(1)  10  0  970  0  80  0 
Represents  atmosphere  having  a  temoeraturs  t0 
C.  a  barometric  pressure  of  970  mb  and  a  rela¬ 
tive  humidity  of  80%. 


vVhen  a  'MET  command  is  encountered  in  ATR, 
the  coda  catenates  a  r.ew  grojnd  e>e  .tion  and 
ovemdes  any  othe*  speef  ed  : — :-.,yo  me 
* GROUND  command.  ”  c  re*  y  me. van 
s  se  (acted  such  th-  :  toe  mean  ?>  -  cans.ty  •:otre- 
soonds  to  that  midway  between  source  and  de¬ 
tector.  This  is  done  to  insure  tna!  a  realistic 
atmospheric  density  profile  exists  b3tween 
source  and  detector.  The  existence  of  such  a 
prctiie  s  imoortant,  because  the  density  at  the 
source  location  is  used  m  A'R  to  calculate  the 
fireoall  radius,  which  ts  pertinent  to  delayed  radi¬ 
ation  intensity. 

The  revised  ground  elevation  is  calculated 
using  the  expression: 

Z,  =  (  -13800  x  In  (  ,i  /  (,0)  i1  ’  04  .  ,  z,  ♦ 
Zd  )/ 2 


(D 


whe^e 


Example  2  ‘MET (4)  1  1889E-3  7  507E-8 
Represents  an  atmosphere  having  a  total  ( moist t 
air  density  of  1  889  *  10  3  g/cc  and  a  water  vapor 
density  of  7  507  *  10  9  g/cc  and  >s  the  eourvaiem 
of  that  atmosphere  specified  m  Exar.ipie  i 

Some  useful  relationships  between  ihe  various 
atmospheric  parameters  are  provided  below  and 
may  be  considered  typical  of  a  range  of  condi¬ 
tions  to  be  found  at  an  elevation  of  1  im  above 

sea  level 

MO'St  Air 


:i.£i 

PlrrJu 

Density 

U.  UP3. 

Wt% 

iVaiai 

30 

899 

50 

1  024E  3 

1  47 

25 

899 

50 

1  044E  4 

1  10 

20 

899 

50 

1  060E-3 

0  812 

15 

899 

50 

1  083E  3 

0  593 

10 

899 

50 

1  103F  3 

0  4?7 

5 

899 

50 

u, 

C* 

0  303 

0 

899 

50 

1  1 45f  3 

0  212 

Z,  s  the  ground  e»evation  above  sea  tr-vet  in 

me’ers 

P  is  the  nean  aw  density  ig  cc) 

p0  is  tli*  sea  levef  air  ctensity  (1.225E-03 

g  cc) 

Z,  is  the  source  height  above  the  ground  m 

meters 

Zd  s  the  deteetty  height  above  the  grcimd 
in  meters 

Ground  'eveis  calculated  using  this  expression 
r  orrespond  to  appropriate  elevation*,  as  speci- 
hod  by  me  U  S  Standard  Atmosphere.  within  a 
few  percent  tor  source  detector  mean  elevations 
above  sea  level  to  5000  meters  input  combina¬ 
tions  of  source  and  detector  heights  and  mean 
air  densities  which  result  in  ground  elevations  be¬ 
low  sea  level  are  not  allowed  Heyond  10.000  me¬ 
ters  the  accuracy  of  this  expression  is  worse  than 
a  lew  pm  cent,  but  that  is  not  significant  m  relation 
in  its  use  m  ATR.  which  i$  to  determine  whether 
the  firehell  touches  the  ground 

2  2  FISSILE  MATERIAL  SPECIFICATION. 

Section  5  of  this  report  describes  revisions  msde 
m  the  ATR  delayed  radiation  model  formerly,  all 


fission  product  radia-ioo  was  baaed  on  photon 
emission  from  the  fission  products  of  a  single  fis¬ 
sile  material  U-225.  However,  as  wilt  be  de¬ 
scribed  later  in  this  ropod,  the  model  has  been 
revised  tc  accommodate  delayed  photon  emis¬ 
sion  from  the  fission  products  oi  three  fissile  ma¬ 
terials.  U-235,  U-228,  and  Pu-239.  The  inclusion 
of  such  detail  is  important  be  cause,  over  much 
of  the  time  regime  of  interest  to  initial  radiation, 
i he  photon  energy  errostwon  rate  of  U-238  is 
double  that  of  U-235.  which  in  tun  is  twenty-five 
percent  more  than  that  of  Pu-239. 

The  quantity  required  for  spec-tying  the  contribu¬ 
tion  of  each  nuclide  to  the  delayed  gamma-ray 
free  field  is  the  yield  or  the  fraction  of  the  total 
yield  due  to  the  fission  of  each  nuclide.  The  com¬ 
mand  used  tc  accomplish  the  input  of  these 
quantities  is  the  'FP  command.  Fmmerfy.  speci¬ 
fication  of  the  fission  product  source  require  up 
to  two  input  statements: 

"FP-Y  value. 

which  specified  the  total  yield  verue,  in  kiiofons, 
and 

'FP-F  value, 

which  specified  the  fraction  of  thi>  total  yield  due 
to  fusion,  ‘n  the  event  the  "FP-F  command  was 
omitted  a  value  of  1 .0  was  adopted  j  default 
value. 

As  .  evsed,  a  third  "FP  Inp.i  command  has  be«i 
included,  as  follows: 

"FP-1  values. 

which  spec  Ties  the  respective  fraction*  ot  the  f/s- 
?.on  yield  due  fo  U-235.  U-2v8  and  Pu  239  fission, 
IN  THAT  O HDCR  Omission  of  the  'FP-t  command 
results  m  U-235  fi»sion  by  default  F»ampi<*s  of 
revised  fission  product  source  input  commands 
are  as  follow*  ■ 

Example  1 ,  "FP-Y  21 
"FP-F  1  00 
'FP  I  0  0  0  2  0  8 

Represents  «  tcfsl  weapon  yield  N  21  kt.  ICON 
of  which  results  from  lission,  20%  U-238  and 
80%  Pu-239. 


Example  2.  "FP-Y  1000 
"FP-F  0.70 
•FP-I  0.05  0.35  0  0 

Represents  a  total  weapon  yield  of  10C0  kt  (one 
megaton)  of  which  70%  is  produced  by  fission, 
5%  from  U-235  and  95%  from  U-238. 

As  stated  earlier,  this  ra-son  of  ATR  has  inte¬ 
grated  the  delayed  radiat-oo  computat-on  into  'he 
code.  Thus,  delayed  radiation  close  is  now  in¬ 
cluded  automatically  as  all  or  part  cf  the  output 
resulting  from  the  following  commands: 

•  DOSE/FPf 

which  specifies  the  fission  product  dose  ver¬ 
sus  running  geometry  coordinate 

"DOSE/G Gi 

which  specifies  the  total  gamma  ray  dose 
versus  running  geometry  coordinate. 

"DOSE/T/ 

which  specifies  the  total  neutron  plus  gam¬ 
ma-ray  dose  \  srsus  running  geometry  coor¬ 
dinate. 

However,  as  of  this  time  the  "PRINT  and  "WRITE 
command  have  not  been  modified  to  include  *is- 
s*oo  product  gamma  ray  ou’put.  The  primary  rea¬ 
son  for  this  is  that  full  incorporation  of  that 
componen*  would  require  calculation  of  time-m- 
tegrsl.  angle-differential  dux  This  cannot  be 
readily  accomplished  in  the  one-dimensional  'or- 
mat  common  to  »he  other  radiation  components 
because  the  ftssixt  product  source  ooes  not  re¬ 
main  fixed  at  the  original  burst  location  but  rather 
uses  with  the  ascending  fireball  Incorporation  of 
fission  product  gamma  rays  into  the  '  PRINT  and 
"WRITE  commands  <9  expected  to  occur  when 
ATR  is  mod  Fed  to  produce  duence  data  m  a 
two-dimensional  format  That  modification  is  part 
of  the  ongoing  program  to  modify  ATR  to  account 
for  energy-  and  angle  differential  duence  pertur¬ 
bations  st  the  air-ground  interlace,  as  mentioned 
previously 

2.3  ENERGY  GROUP  ARCHITECTURE. 

The  energy  boundaries  of  the  multigroup  neutron 
and  gamma-ray  duence  modeled  in  ATR  have 
been  changed,  although  ;he  22  neutron  -  10  gam¬ 
ma  ray  group  tormat  has  been  retained  The  n«w 
grouo  structure  is  the  compeuble  with  all  major 
ci  oss  section  sets  sponsored  by  the  Defense  Nu- 
cieer  Agency  .Those  include  the  ONA  F  ew-Qroup 
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Gross  Section  Set  (DLC-31 )  fR@f  .4) ,  the  revised 
DNA  Faw-Group  set  now  being  tested 
iDLC- 1 30)  (Pet  . 51 ) .  and  Vitamin  E  iRef .  60) .  The 
new  AIR  energy  bounds  W  neutrons  3nd  gamma 
rays  are  provided  in  Tables  2  and  3,  respectively 
along  with  those  of  CLC-31  and  DLC-130  N c'.e 
mat  me  -p  ^os t  gamma  ray  e-"e' jy  g'ot.c  Poun- 
a','  n  i*j  ev-sea  aTm  ,s  '2  *.-e»  wr-ftie  as  tra* 
of  CH..C  1  i  s  14  MeV.  That  change  was  made  to 
provide  a  more  realistic  differential  spectrum 
(photons  Me V)  in  the  highest  energy  gamma  ray 
group,  since  the  highest  energy  gamma  ray  tc 
be  encountered  ir.  air-ground  transport  is  that  of 
10  3  MeV.  produced  by  neutron  caoture  m  nitro¬ 
gen. 

2.4  INTERNAL  SOURCE  SPECTRA. 

The  source  spectra  traditionally  contained  in  ATR, 
unclassified  fission  and  thermonuclear  weapon 
neutron  spectra  origmal'y  devetooed  by  W  Rig¬ 
gers  of  Los  Alamos  Nat.onal  Laboratory,  and  the 
generic  prompt  gamma-ray  source  spectrum 
have  been  regrouped  into  the  new  energy  bounds 
described  above  Those  sources,  accessible  m 
ATR  through  the  "Z-SOURCE  command.  are  pro¬ 
vided  m  Tanies  4  ana  5  lor  neutrons  and  gamma 
rays,  respectively 

2  3  FLUENCE-TO-DOSE  CONVERSION 
FACTORS. 

The  most  !f»  quently  used  portion  of  the  Outout 
from  an  ATh  calculation  is  that  which  provides  to¬ 
tal  dose.  i<erma.  or  a  "imilar  integral  quantity  as 
a  function  of  range  New  fluence-to-dose  conver¬ 
sion  factors  have  been  incorporated  »n  ATR  .n  an 
eifort  to  more  completely  satisfy  diverse  user 
needs 

Tables  6  and  7  list  the  revised  AIR  flu- 
er.ce-to-dose  conversion  factors  lor  neutrons 
and  gamma  rays  r  esfrectively ,  t  /  energy  group 
Those  conversion  factors  and  their  origin  am  Oe 
smbed  in  detail  in  Appendices  U  and  T)  ol  this 
eport.  beef  summaries  of  which  are  provoed 
hr  tow. 

o  ANSI  Standard  Oose  For  the  Calculation 
cf  Human  Dose  Equivalent.  i  m  .  quati 
t\  factor- weighted  0  <se  these  values 
are  relevant  to  radiation  protection  and 
risK  of  long  term  illness,  such  as  <  am  m , 
from  radiation  enposu re  Dose  eiio'vatenl 
has  no  Known  relevance  to  early  «<lo<  is. 
such  a*  prodromal  syoc.ome  or  early 
mortality. 


o  Soft  Tissue  Kerma  -  Kinetic  energy  of 
charged  particles  releasad  n  soft  tissue 
by  ionizing  radiation  This  quantity  is  used 
to  calculate  bee  ‘^d  tissue  dose  and  re¬ 
places  both  the  soft  hssue  Kerma  and 
-enderson  dose  con  •  ersic'1  'actors  n 
»wous  versions  oi  Arp  T-e  -onoerson 
lose  convers'cn  fact..;  was  scttu  that 
‘or  a  tour  element  tissue  mode'  -owron- 
r  dered  obsolete. 

o  Mio-head  and  Mid-phantom  Dose  -  F!u- 
ence-to-soft  tissue  dose  conversion  fac¬ 
tors  which  include  the  effect  of 
transmission  into  the  mid-head  arid 
mid-torso  locations  in  a  70  Kg  anthropo¬ 
morphic  phantom. 

o  Concrete  Kerma  -  Kinetic  energy  of 
charged  particles  relec  ed  in  concrete  by 
■onizing  radiation. 

o  Dry  Air  Kerma  -  Kinetic  energy  o‘  charged 
particles  released  m  dry  air  by  ionizing  ra¬ 
diation  Dry  Air  Keima  <n  units  of  rad3 
may  be  converted  to  units  of  roentgens 
by  dividing  the  rad  values  bv  0  37 

o  Non-Ionizing  Silicon  Kerma  -  Kmeiic  ener¬ 
gy  released  m  silicon  crystal  by  neutrons 
which  is  associated  with  crysta'ime  struc¬ 
ture  changes  rntner  than  rreat-nn  of  free 
electrons. 

r  1  MnV  Equivalent  Damage  F'uence  -  The 
fiuence  of  I  MeV  neutron  Past  reactor 
teenage  neutrons,  having  a  mean  energy 
if  app'oumately  t  MeV!  required  to  pro¬ 
duce  the  same  amount  of  non-ionizing 
S'licon  Kerma-'«iated  darraqe 

n  lom/ing  Silicon  Kerma  k  inane  energy  i  e 
leaser/  m  silir  nn  crystal  by  neutrons  and 
■gamma  rays  m  the  tnrm  of  electron  hole 
pairs 

The  'rXTSf  command  will  r;ause  ATR  to  calculate 
and  display  value*  for  all  the  above  (loan  quanti¬ 
ties  in  addition  the  'FlUXWT  Command  allows 
the  user  to  provide  up  10  5  sets  of  additional  Hu 
wim  to  dose  conversion  factors  as  nput  to  A  T R 
The  '  -st  set  of  fluence-to-dose  factors  entered 
in  (tie  revised  ATR  via  the  ’FLUXWT  command 
wilt  produce  integral  vatuea  displayed  as  pad  of 
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Table  2.  ATR  detector  energy  group  boundaries  for  neutrons  (MeV)  with 
DLC-31  and  DLC-130  group  boundaries  for  comparison. 


Group 


ATR5 

Lower  Upper 

Energy  Energy 


Group 


DLC-31 

Lower  Upper 
Energy  Energy 


DLC-130 

Lower  Upper 

Energy  Energy 


1. 

1.0C( 

-n 

-  4.14( 

;-07) 

1. 

1.00 

-n 

2. 

4.14| 

-07 

-  1.13 

-06 

2. 

4.14 

-07) 

3. 

1.13 

-06 

-  3.06 

-06 

3. 

1.13 

-06) 

4. 

3.06 

-06 

-  1.07 

-05 

4. 

3.06 

-06 

5. 

1.07 

-05 

-  2.90 

-05 

5. 

1.07 

-05) 

6. 

2.90 

-05 

-  1.01 

-04 

8. 

2.90 

-05) 

7. 

1.01 

("04; 

-  1.23 

-03 

7. 

1.01 

(-04) 

4.14 

-07 

1.13 

-06 

3.06 

-06 

1.07 

-05 

2.90 

-05 

1.01 

-04 

5.83 

-03 

8. 

1.23 (-03) 

-  2.19(-02) 

9. 

2. 19 (-02) 

-  1.11 (-01) 

10. 

1.11 (-01) 

-  1.58(-01) 

11. 

1.58(-0l) 

-  5.50(-0l) 

12. 

5.50(-01) 

-  1.11 (+00) 

13. 

l.ll(+00) 

1  -  1.83(400) 

14. 

1.83) 

>0°] 

-  2.311 

;4°°) 

15. 

2.311 

(-00) 

-  3.011 

l+oo) 

16. 

3.011 

1*00] 

-  4.071 

>00) 

17. 

4.071 

(400’ 

-  4.971 

(+00) 

18. 

4.971 

(4°°; 

-  6.381 

( 4°0) 

19. 

6.381 

(-00) 

-  8.191 

(400) 

20. 

8. 19(+00)  -  1.00(+01) 

21. 

1.00(+01)  -  1.22(+01) 

22. 

1.22(401)  -  1.49(401) 

8. 

5.83 

(-03) 

-  1.23 

(-03) 

9. 

1.23 

—03 

-  3.35 

-03) 

10. 

3.35 

-03 

-  1.03 

-02) 

11. 

1.03 

-02 

-  2.19 

-02) 

12. 

2.19 

-02 

-  2.48 

-02) 

13. 

2.48 

(-02) 

-  5.25 

(-02) 

14. 

5.25 

(-02) 

-  1.11 

(-01) 

15. 

1.11 

-01 

-  1.58 

-Ol) 

16. 

1.58 

(-01) 

-  5.50 

(-0l) 

17.  5. 50 (-01)  -  l.ll(-OO) 


18.  1.11(+00)  -  1 .83 (+00) 


19. 

1 . 83  ( 

400) 

- 

2.31 

'+00) 

20. 

2.31  ( 

400 

- 

2.39 

+00) 

21. 

2.39( 

'400 

- 

3.01 

+0°) 

22. 

3.01 

>00 

- 

4.07 

+00 

23. 

4.07 

'+00 

- 

4.72 

+00 

24. 

4.72( 

(+00) 

4.97 

+00) 

25. 

4.97| 

400 

- 

6.38 

+00) 

26. 

6.38| 

400 

- 

7.41 

+00) 

27. 

7.41 

400 

- 

8.19 

+0°) 

28. 

8.19 

(+00) 

- 

9.05 

+0°) 

29. 

9.05 

+00 

- 

1.00 

+01) 

30. 

1.00 

+01 

- 

1.11 

+01) 

31. 

1.11 

+01 

- 

1.22 

+01 ) 

32. 

1.22 

+01 

- 

1.28 

+01) 

33. 

1.28 

+01 

- 

1.38 

+01 ) 

34. 

1.38 

+01 

- 

1.42 

+0l) 

35. 

1.42 

i+0l) 

- 

1.49 

(+0l) 

36. 

1.49 

f+01 

- 

1.69 

(♦oil 

37. 

1.69 

(+01 

- 

1.96 

(+01) 

1. 

1.001 

-11) 

- 

4.14 

-07) 

2. 

4.14 

-07 

- 

1.13 

'-06) 

3. 

1.13 

-06) 

- 

3.06 

'-06 

4. 

3.061 

'-06) 

- 

1.07 

-05) 

5. 

1.07 

-05 

1  - 

2.90 

'-05) 

6. 

2.90 

-05 

1  - 

1.01 

'-04 

7. 

1.01 

'-04 

1  - 

2.75 

-03) 

8. 

2.75 

'-03 

- 

5.83 

-03) 

9. 

5.83 

'-03 

- 

1.23 

-03) 

10. 

1.23 

(-03 

- 

3.35 

-03 

11. 

3.35 

-03 

- 

1.03 

-02) 

12. 

1.03 

-02 

- 

2.19 

-02) 

13. 

2.19 

-02 

- 

2.48 

-02) 

14. 

2.48 

-02 

- 

3.43 

-02) 

15. 

3.43 

-02 

- 

5.25 

-02 

18. 

5.25 

-02 

- 

1.11 

^-01) 

17. 

1.11 

-°! 

- 

1.58 

(-01 

18. 

1.58 

-01 

- 

2.47 

-Oi 

19. 

2.47 

-01 

- 

3.68 

-Oi 

20. 

3.68 

-01 

- 

5.50 

-Oi 

21. 

5.50 

-01 

- 

6.39 

-Oi 

22. 

6.39 

-01 

- 

7.43 

-01 

23. 

7.43 

-oi 

- 

8.21 

-Oi 

24. 

8.21 

-01 

- 

9.62 

-Oi 

25. 

9.62 

-01 

- 

1.11 

+00) 

26. 

1.11 

+0° 

- 

1.42 

+oa) 

27. 

1.42 

+0° 

- 

1.83 

+00 

28. 

1.83 

(+0°) 

- 

2.31 

+00) 

30. 

2.31 

+00 

- 

2.39 

+00 

31. 

2.39 

+00 

- 

3.01 

+0°) 

32. 

3.01 

+00 

• 

4.07 

+00 

33. 

4.07 

+00 

- 

4.72 

+ 

0 

0 

34. 

4.72 

+00 

- 

4.97 

+00) 

35. 

4.97 

+00 

- 

6.38 

+0°) 

36. 

6.38 

00 

- 

7.41 

+00) 

37. 

7.41 

+U0 

- 

8.19 

+00 

37. 

8.19i 

+00 

- 

9.05 

+00) 

38. 

9.051 

+00 

- 

1.00 

+01) 

39. 

l.OOi 

(+01 

- 

1.11 

+oij 

40. 

1 .11 1 

>01 

- 

1.22 

+0! 

41. 

1.22i 

♦01 

- 

1.25 

+01) 

42. 

1.25 

♦01 

)  - 

1.381 

[+01 

43. 

1.38 

(+01 

)  - 

1.42i 

+01 

44. 

1.42 

(+01  ( 

!  - 

1 .49' 

(+01) 

45. 

1.49 

(+or 

)  - 

1.69 

(+oi) 

46. 

1.69 

(+01  ( 

)  - 

1.96 

(+01) 
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Table  3.  ATR  source  and  detector  energy  group  boundaries  for  prompt 
gamma  rays,  and  detector  energy  group  boundaries  for  secondary 
gamma  rays  (MeV)  with  DLC-31  and  DLC-130  group  boundaries  for 
comparison . 


ATR  DLC-31  DLC-130 

Lower  Upper  Lower  Upper  Lower  Upper 

Group  Energy  Energy  Group  Energy  Energy  Group  Energy  Energy 


1.  1 .00 (-02)  -  4. 50 (-02)  1.  1.00 

(-02' 

-  2.00 

(-02' 

1.  1.00 

(-02 

-  2.00 

(-02) 

2.  2.00 

[-02* 

-  3.00 

-02 

2.  2.00 

-02* 

-  3.00 

-02) 

3.  3.00 

-02* 

-  4.50 

-02 

3.  3.00 

-02 

-  4.50 

-02) 

2.  4.50(-02)  -  1.00(-01)  4.  4.50 

-02* 

-  7.00 

-02 

4.  4.50 

-02* 

-  7.00 

'-02) 

5.  7.00 

-02 

-  1.00 

-oi 

5.  7.00 

-02* 

-  1.00 

-01 ) 

3.  1.00 

(-01' 

-  1.50 

(-or 

6.  1.00 

-01* 

-  1.50 

-oi 

8.  1.00 

-Oi 

-  1.50 

-01) 

4.  1.50 

-01 

-  3.00 

-01 

7.  1.50 

-or 

-  3.00 

-Oi 

7.  1.50 

-m 

-  3.00 

-on 

5.  3.00 

-oi 

-  4.50 

-°! 

8.  3.00 

-01 

-  4.50 

-°i 

8.  3.00 

-0i 

-  4.50 

-01) 

6.  4.50 

-01 

-  7.00 

-01 

9.  4.50 

-01 

-  7.00 

-0i 

9.  4.50 

-Oi 

-  7.00 

-01) 

7.  7.00 

-°i 

-  1.00 

>00 

10.  7.00 

-01 

-  1.00 

>00 

10.  7.00 

-01 

-  1.00 

>00) 

8.  1.00 

>00 

-  1.50 

>00 

11.  1.00 

>00 

-  1.50 

>00 

11.  1.00 

>00 

-  1.50 

>00) 

9.  1.50 

>00' 

-  2.00 

>00 

12.  1.50 

>00 

-  2.00 

>00' 

12.  1.50 

>00 

-  2.00 

>oo) 

10.  2.00 

>00 

-  2.50 

>00 

13.  2.00 

>00 

-  2.50 

>00 

13.  2.00 

>00 

-  2.50 

>00 

11.  2.50 

>00 

-  3.00 

>00 

14.  2.501 

>00 

-  3.00 

>00' 

14.  2.50| 

>00 

-  3.00 

>00) 

12.  3.00 

>0° 

-  4.00( 

'>00] 

15.  3.001 

>00 

-  4.00| 

'>00 

15.  3.00| 

>00 

-  4.G0| 

>00) 

13.  4.00 

>00' 

-  5.00| 

'>00 

16.  4 . 00 < 

'>00 

-  5.00 

'>00 

16.  4.00| 

>00' 

-  5.00( 

'>00) 

14.  5.00| 

>00 

-  6.00, 

>00 

17.  5 . 00 ( 

>00 

-  6.00( 

>00 

17.  5.00 

'>00 

-  6.00( 

>00) 

15.  6.001 

>00' 

-  7.00( 

'>00 

18.  6.00( 

>00 

-  7.C0( 

>00' 

18.  8.00 

'>00 

-  7.00( 

>00) 

16.  7.00| 

>00 

-  8.00( 

>00 

19.  7.00( 

>00 

-  8 . 00  ( 

>00 

19.  7.00( 

>00 

-  8.00( 

>oo) 

17.  8.00| 

>00' 

-  1.00( 

>01 

20.  8.00( 

>00 

-  1 .00 

>01 

20.  8.00( 

>00* 

-  1.00( 

>01) 

18.  1.00( 

:>oij 

-  1 ,20( 

>01) 

21.  1.00( 

>01) 

-  1.40( 

>01) 

21.  1 .00i 

>01 

-  1 . 20  ( 

>01) 

22.  1.20i 

>01) 

-  1 . 40 ( 

>01) 

23.  1.401 

>01) 

-  2.00( 

>01) 
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Table  4.  Energy  distributions  for  ATR5  neutron  sources. 


Fission  Weapon 

Energy  Boundaries  *N-S0URCE(1) 

Croup  (MeV)  (Fraction  in  Group) 


1 

1.07 (-5) 

- 

2.90(-5) 

0.0 

2 

2. 90 (-5) 

- 

1.01 (-4) 

0.0 

3 

1.01(-4) 

- 

1.23 (-3) 

0.0 

4 

1 . 23(-3) 

- 

2. 19(-2) 

0.01649 

5 

2 . 19(— 2) 

- 

l.ll(-l) 

0.20617 

6 

0.111 

- 

0.158 

0.01799 

7 

0.158 

- 

0.550 

0.15129 

8 

0.550 

- 

1.11 

0.21587 

9 

1.11 

- 

1.83 

0.14678 

10 

1.83 

- 

2.3i 

0.10173 

11 

2.31 

- 

3.01 

0.03871 

12 

3.01 

- 

4.07 

0.05480 

13 

4.07 

- 

4.97 

0.01177 

14 

4.97 

- 

6.36 

0.01832 

15 

6.36 

- 

8.19 

0.01274 

16 

8.19 

- 

10.0 

0.00734 

17 

10.0 

- 

12.2 

0.0 

18 

12.2 

- 

15.0 

0.0 

Thermonuclear 
*N-SQURCE(2) 
(Fraction  in  Group) 


0.0 

0.00200 

0.05719 

0.34417 

0.10964 

0.01142 

0.09058 

0.08500 

0.06200 

0.02592 

0.02608 

0.02600 

0.01700 

0.01800 

0.01470 

0.01410 

0.02560 

0.07060 
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Table  5. 


Energy  distribution  for  the  ATR5  prompt 
fission  gamma-ray  source. 


Group 


Prompt  Fission 

Energy  Boundaries  *G-SQURCE(1) 

(MeV)  (Fraction  in  Group) 


1 

0.01 

- 

0.045 

0.03259 

2 

0.045 

- 

0.10 

0.01644 

3 

0.10 

- 

0.15 

0.04881 

4 

0.15 

- 

0.30 

0.10321 

5 

0.30 

- 

0.45 

0.13571 

6 

0.45 

- 

0.70 

0.20256 

7 

0.70 

- 

1.00 

0.16332 

8 

1.0 

- 

1.5 

0.14073 

9 

1.5 

- 

2.0 

0.06429 

10 

2.0 

- 

2.5 

0.03743 

11 

2.5 

- 

3.0 

0.02225 

12 

3.0 

- 

4.0 

0.02109 

13 

4.0 

- 

5.0 

0.00746 

14 

5.0 

- 

6.0 

0.00265 

15 

6.0 

- 

7.0 

0.00092 

16 

7.0 

- 

8.0 

0.00038 

17 

8.0 

- 

10.0 

0.00016 

18 

10.0 

- 

12.0 

0.0 

18 
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the  *  DOSE  output.  Previously,  the  ATR  user  could  tries  through  the  use  of  summary  tables  which 

only  obtain  access  to  the  results  of  'FLUXWT  en-  ar9  pad  of  the  'PRINT  command  output. 
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SECTION  3 

UNIFORM  AIR  TRANSPORT  DATA  BASE 


Modification  of  ATR  has  included  the  replacement 
of  tne  source  detector  energy-deferential  trans¬ 
port  data  modules  for  neutrons  secondary  gam¬ 
ma  rays  and  prompt  gamma  rays.  Angular 
fluence  data  have  also  been  modified  to  be  con¬ 
sistent  with  the  scalar  fluence  information, 

3.1  DATA  BASE  CALCULATION. 

The  transport  data  base  for  the  ATR  model  is 
comprised  of  downscatter  spectra  in  22  neutron 
and  18  gamma-ray  fluence  energy  groups  from 
18  neutron  and  18  gamma-ray  source  groups. 
These  data  have  been  produced  using  the  ANISN, 
one-dimensional,  discrete  ordinates  transport 
code  (Ref.  20).  Calculations  have  been  per¬ 
formed  using  an  S40  Gaussian  quadrature  and 
cross  sections  from  the  VITAMIN-E  coupled 
neutron-gamma-ray  Library  having  a  P5  Legendre 
scattering  approximation  for  neutrons  and  Pg  for 
gamma  rays. 

The  VITAMIN-E  cross  section  i.brary  has  a  format 
consisting  of  coupled  neutron-gamma-ray  multi¬ 
group  data  in  1  74  neutron  and  38  gamma-ray  en- 
ergy  groups.  The  energy  structure  of  VITAMIN-E 


neutron  cross  sections  is  comprised  of  nearly 
equal  lethargy  increments  from  2G  MeV  to  0  1  eV 
with  one  additional  group  below  0. 1  eV.  The  con¬ 
version  to  ATR  22  neutron  -  18  gamma-ray  group 
format  was  accomplished  by  .  dilating  source 
values  for  all  VITAMIN-E  groups  constituent  to  a 
specified  ATR  source  group,  such  that  the  sum 
of  the  sources  was  one  particle,  performing  the 
transport  calculation  in  uniform  air  using  the  full 
VITAMIN-E  group  structure  and  summing  fluence 
data  within  ATR  energy  bounds.  The  in-group 
source  weighting  for  gamma-ray  source  groups 
was  based  on  a  1/E  source  shape.  The  in-group 
source  weighting  for  neutrons  was  based  on  the 
ATR  internal  thermonuclear  source  spectrum  In 
DLC-31  format  (37  neutron  groups). 

The  calculations  were  carried  out  in  spherical  ge¬ 
ometry  to  an  optical  depth  of  550  grams/cm2  of 
moist  air  as  specified  in  Table  8.  having  a  density 

of: 

1.115  '  10'3  g/cc  dry  air 

6.330  *  10  s  g/CC  water  vapor 

1.121  *  10-3  g/cc  total  moist  air 


Table  8.  Data  base  moist  air  element  specifications. 


□QmanL_i..Madium 

^QlOtllla 

Cium£ifli_DanaitxlAi 

H 

W*t 

6  317F  02 

4.236F  07 

O  - 

Wrtter 

5  01  ?F  01 

2.1  17E  07 

0  - 

Air 

2  30?E*01 

9.722E  06 

N  • 

Air 

7  5ME+01 

3.625E-05 

A  r  • 

Air 

1  283E  *00 

2.1  70E-07 

22 


Thus,  the  calculations  extended  to  a  radius  of  ap¬ 
proximately  5  km.  Only  the  data  to  500  gram/cc 
(approximately  4.5  km)  were  used  to  establish 
the  data  base  for  ATR.  in  order  to  allow  foi  pertur¬ 
bations  caused  by  the  vacuum  at  the  outer 
boundary  of  the  calculation. 

3.2  DATA  BASE  PARAMETERIZATION. 

It  would  be  possible  to  create  an  air  transport 
model  simply  by  saving  and  performing  interpola¬ 
tion  operations  on  a  complete  set  of  differential 
fluence  data  from  one-dimensional  discrete  ordi¬ 
nate  calculations.  However,  such  a  system  would 
be  very  large  and  unwieldy.  Similarly,  those  data 
could  be  fit  with  high  order  polynomial  functions. 
Again,  the  size  of  the  coefficient  data  array  would 
be  very  large. 

The  modeling  approach  used  in  ATR  is  to  use  effi¬ 
cient  functional  forms  to  fit  the  data  base  and  to 
separate  variable  relationships  so  that  a  minimum 
number  of  computations  have  to  be  performed 
to  reconstruct  the  required  fluence  data.  Pre¬ 
vious  versions  of  ATR  have  incorporated  separate 
functional  fits  of  total  fluence.  energy-differential 
fluenco.  and  angle-differential  fluence  data,  This 
approach  was  taken  because  the  total  fluence 
was  easily  fit  with  few  coefficients,  while  the  dif¬ 
ferential  fluence  models  could  be  based  on  data 
normalized  to  total  fluence.  That  approach  al¬ 
lowed  energy-  and  angle-differential  models  to  be 
constructed  at  specific  distances,  with  interpola¬ 
tion  to  intermediate  distances,  rather  than  as 
continuous  functions  ol  energy. 

The  functional  form  used  in  previous  versions  of 
ATR  to  fit  the  data  as  described  above  is  a  poly¬ 
nomial  series  comprised  of  haif-integer  and  loga¬ 
rithmic  terms,  known  as  Webster  Functions.  Even 
though  these  functions  provided  great  flexibility 
in  fitting  complex  data,  it  was  found  that  the  dif¬ 
ferential  data  could  not  be  adequately  fit  with  one 
polynomial  series.  However,  it  was  discovered 
that  the  ratios  between  the  first  order  functional 
fit  and  the  data  Itself  were  more  amenable  to  the 
fitting  procedure.  Thus. previous  versions  of  ATR 
contain  one  or  two  iterative  correction  routines 
to  improve  the  quality  of  the  fit.  likewise,  when 
revisions  have  been  made  to  the  code,  additional 
corrective  layers  have  been  applied  rather  than 
refitting  the  data  base  from  scratch. 

There  has  been  one  previous  maior  revision  to 
the  ATR  data  base,  at  which  time  the  original  data 
were  modified  to  reproduce  a  sel  computed  us¬ 
ing  improved  cross  sections.  It  was  decided  ihat 


adding  still  another  layer  of  modifications  would 
increase  the  running  time  of  the  code,  further 
complicate  its  already  complicated  interna!  struc¬ 
ture,  and  possibly  provide  an  unreliable  fit  to  the 
data.  Thus,  it  was  decided  to  fit  the  one¬ 
dimensional  transport  data  base  starting  from 
scratch  in  the  current  revision  ol  ATR. 

The  approach  to  data  parameterization  in  the  cur¬ 
rent  revision  of  ATR  is  changed  from  that  used 
in  previous  versions.  Because  the  primary  quanti¬ 
ty  used  in  ATR  for  calculation  of  range-dependent 
do3e  or  kerma  is  energy-differential  fluence, 
there  was  no  particular  reason  to  fit  the  total  flu¬ 
ence  separately,  except  for  the  sake  of  conve¬ 
nience.  Thus,  the  approach  was  taken  to  directly 
parameterize  energy-differential  fluence  as  a 
function  of  distance. 

As  revised.  ATR  computes  energy-differential  flu¬ 
ence  using  two  models.  First,  the  code  calculates 
the  uncollided,  energy-differential  fluence  using 
the  total  cross  section  applicable  to  each  source 
energy  group.  Second,  the  code  computes  the 
scattered  fluence  contribution.  The  model  used 
to  fit  scattered  fluence  data,  including  in-source 
group  scatter,  is  based  on  the  fact  that  the  func¬ 
tion 

Function  =  F (pir)  x  (e-f‘r)  /  r2  (2) 

where  r  is  distance. 

(i  is  the  inverse  relaxation  length,  and 

F(jir)  is  a  polynomial, 

which  is  often  used  to  approximate  the  variation 
of  kerma  oi  fluence  with  distance,  taxes  the  form 
of  a  Taylor  series  on  r  with  powers  in  half-integer 
steps. 

It  was  found  that  a  series  truncated  at  powers  of 
t  3.  including  twelve  functional  terms  and  a  con¬ 
stant,  could  reproduce  the  energy-  differential 
fluence  data  with  great  accuracy  within  one  or 
two  percent  at  each  location.  Through  multiple 
trials  a  subset  of  six  terms  of  the  original  thirteen 
was  found  which  could  reproduce  the  data  nearly 
as  well,  certainly  within  the  t  5  percent  limits 
which  were  the  target  criteria.  Thus,  a  software 
system  was  oevised  which  performed  a  least 
squares  fit  to  all  possible  combinations  of  six  of 
the  available  thirteen  terms,  and  provided  the  re¬ 
sults.  in  terms  of  the  resulting  coefficients,  in  or¬ 
der  of  descending  coefficient  of  determination 
value,  Top  combinations  of  coefficients  were  ten¬ 
tatively  chosen  for  Inclusion  in  the  system.  In 
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those  cases  where  several  combinations  were 
ranked  as  having  equal  values  (to  seven  signifi¬ 
cant  figures)  of  the  coefficient  of  determination, 
that  combination  providing  the  maximum  sum  of 
all  partial  correlation  coefficients  was  chosen 
The  sets  of  coefficients  were  then  individually 
tested  to  insure  that  each  reorcduced  its  data  set 
to  ♦  5  percent  or  oettei  at  almost  ail  locations 


Examples  of  fit  fluence  plotted  together  with  the 
data  from  wh-ch  they  were  derived  are  presented 
in  Figures  7.  8  and  9  for  neutron,  secondary  gam¬ 
ma  ray  and  prompt  gamma  ray  kerma.  respec¬ 
tively.  from  a  fission  source. 


The  angular  ‘luence  parameterization  was  re¬ 
tained  from  ATR  4 


3.3  TRANSPORT  PERTURBATION  BY  AIR 
MOISTURE  VARIATION. 


The  atmosphere  is  nominally  considered  to  con¬ 
tain  nitrogen,  oxygen  and  argon  m  descending  or¬ 
der  of  ahundance.  with  the  last  mentioned  being 
present  at  levels  only  a  little  over  a  percent  of  the 
totai  mass  The  argon  13  important  cause  of  its 
large  photoelectric  cross  section  which  helps 
make  the  atmosphere  opaque  to  photons  below 
approximately  TO  koV  Those  three  constituents 
make  up  'he  so-called  dry  atmosphere,  which 
has  long  been  the  standard  for  use  <n  radiation 
transport  calculations  However.  th«  atmosphere 
also  contains  moisture  m  amounts  which  vary 
over  time,  ranging  from  a  few  tenths  of  a  percent 
to  as  much  as  one  and  a  half  pan,*,''*  ny  weight. 
Consequently  the  atmosphere  generally  contains 
up  to  a  few  tenths  of  percent  h-.drogen  atoms  try 
weight  or  up  to  mote  than  three  percent  by  num 
her . 


As  described  m  section  2  1  .  ATRS  has  heen  given 
the  capability  to  accept  atmospheric  humidity  as 
a  meteorological  variable  This  capability  is  based 
on  scaling  factors  which  have  been  incorporated 
into  ATRS  Those  scaling  factors  have  been  de¬ 
veloped  from  a  set  of  ad|Oint  ANISN.  one dimen¬ 
sional  iranapofl  calculations,  performed  to  a  dry 
air  optical  depth  of  *>5  gram  per  square  centime¬ 
ter.  using  DLC  31  cross  sections  (S40.  P-p.  hav¬ 
ing  moisture  contents  as  follows: 


WaiahLParcsnt 

0.0(Moisture  content  default  if  no 
moisture  specification  is  input) 

0.1 
0  2 
0.4 

0.5645(ATR5  Uniform  Air  Data 
Base) 

0.8 

1.2 

1.6 

2.0 

The  calculations  provided  tissue  kerma  as  a  func¬ 
tion  of  distance  from  each  of  37  neutron  and  21 
gamma  ray  source  energy  groups  These  weie 
collapsed  into  the  18  neutron.  t8  gamma  ray 
group  source  energy  structure  of  ATR5  using  the 
thermonuclear  source  shape  for  neutrons  and 
secondary  gamma  rays  and  inverse  energy  tor 
gamma  rays,  as  described  in  section  3  2.  The  re¬ 
sults  for  each  source  energy  were  htted  as  a 
function  of  optical  depth  t.s  ratio?  of  kerma  calcu¬ 
lated  for  a  particular  mois'ure  content  and  that 
calculated  for  the  ATR5  data  base.  As  applied  in 
ATR5.  the  natural  logarithm  of  the  moisture  con- 
tent  correction  factor  was  expressed  ns  a  func¬ 
tion  of  oohcal  depth  using  the  same  polynomial 
form  used  to  fit  the  uniform  air  fluence  data  base 
That  *nrm  is  a  polynomial  senes  m  halt  integer 
power  steps  from  a  negative  three  10  a  positive 
three  Trials  were  run  to  determine  the  minimum 
number  of  coefficients  required  to  fit  the  ratios 
to  (letter  than  live  percent  accuracy  Those  do- 
tnrm-ned  that  four  tern  s  (plus  a  constant)  were 
required  to  fit  the  neutron  and  secondary  gamma 
my  data,  while  onl'-  Ihree  terms  (plus  a  constant) 
were  required  to  fit  the  gamma  ray  data  As  in 
the  case  of  the  uniform  an  fluences  while  each 
lit  consists  of  four  or  five  terms  they  are  not  nec¬ 
essarily  the  same  four  or  live  terms  for  each 
source  energy 

Mo’Sture  corrections  m  the  form  of  kerma  ra'ios 
for  various  moisture  contents  a.e  depicted  m  Fig¬ 
ures  10.  t  1  and  12  for  fission  neutrons  and  sec¬ 
ondary  gamma  rays  and  fusion  neutrons, 
respectively  It  should  he  noted  that,  as  aryonne 
has  little  etlect  on  the  transport  of  neutrons. 


24 


B»«2-5dFT  TISSUE  KEfMIH  (KM2-RflOS> 


SLANT  RANGE  (KM) 


PROMPT  GRMMR  KERHfl  FROM  TN  NEUTRON  SOURCE  3  IN  INF  .HR 
RATIOS  Of  RTR5  TO  ANISN  ORTA 


Figure  9.  Cmm»  rsy  tissue  kers*  f ros  n  boosted  fission  gnats* 
r»y  source  in  uniform  sir,  ATH5  nmi  its  data  base. 


Moisture  Correction  (Kerma/Kerma) 


Figure  11.  Secondary  gamma  ray  kerma  correction  for  atmos¬ 
pheric  moisture  content  in  uniform  air,  fission 
neutron  source. 
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Moisture  Correction  (Kerma/Kerma ) 


f 


Figure  12.  Neutron  keraa  corrections  for  atnospheric  aoisture 
content  in  unifor*  air,  fusion  (12.2  to  14.9  MeV) 
source . 
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moisture  has  little  effect  on  the  transport  of 
source  gamma  rays.  However,  moisture  has  a 
considerable  effect  on  the  transport  of  neutrons. 
The  nature  of  this  effect  is  to  strip  the  low  energy 
(less  than  1  MeV)  neutrons  from  the  transported 
spectrum.  A  fission  weapon  source  contains  a 
large  number  of  such  neutrons.  Hence,  its  trans¬ 
port  is  significantly  effected  initially,  much  more 
so  than  the  fusion  source.  After  its  initial  loss  of 
low  energy  source  neutrons,  the  transported  fis¬ 
sion  spectrum  stabilizes,  aided  by  the  oxygen 
cross  section  minimum  at  approximately  2.3 
MeV. 

The  prompt  gamma  rays  and  secondary  gamma 
rays  from  fusion  neutrons  are  essentially  unaf¬ 
fected  by  atmospheric  moisture.  However,  the 
effect  of  increasing  moisture  on  secondary  gam¬ 
ma  rays  from  a  fission  neutron  source  is  to  de¬ 


crease  the  distance  from  the  source  at  which 
such  gamma  rays  are  produced.  This  increases 
the  gamma  ray  kerma  at  short  distances,  as 
shown  in  Figure  1 1 .  At  very  large  distances  the 
effect  becomes  less  important  because  the 
source  becomes  more  and  more  equivalent  to 
a  point,  regardless  of  where  the  gamma  rays 
were  produced. 

Note  that  the  moisture  transport  correction  used 
in  ATR5  is  based  on  ratios  of  kerma  for  each 
source  energy  and  is  not  detector  energy  - 
dependent.  Thus,  the  spectrum  is  unchanged 
from  that  contained  in  the  data  base  which  con¬ 
tains  approximately  0.6%  moisture  by  weight. 
This  corresponds  to  a  relative  humidity  of  approx¬ 
imately  50%  on  a  cool  day  ( 1 5  c>  and  is  typical 
of  conditions  which  are  to  be  found  in  central  Eu¬ 
rope  in  the  spring  (Ref.  40). 
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SECTION  4 


AIR-GROUND  INTERFACE  PERTURBATION  MODEL 


Radiation  is  scattered  and  absorbed  in  the  ground 
to  a  different  extent  than  m  air.  Thus,  the  pres¬ 
ence  of  the  air-ground  interface  causes  radiation 

fluences  to  be  perturbed  in  energy  and  angle 
from  those  which  would  otherwise  exist  at  the 
same  distance  from  the  source  in  uniform  air.  At 
the  present  time,  the  ATR  code  does  not  correct 
the  differential  fluence  for  the  presence  of  the 
ground.  Rather,  it  corrects  the  free  field  tissue 
kerma,  an  integral  quantity,  based  on  ratios  of 
two-dimensional  calculations  with  and  without  the 
ground. 

It  is  important  to  note  that  differential  fluences 
produced  by  ATR  for  detector  locations  near  the 
air-ground  interface  are  those  for  uniform  air.  ad¬ 
justed  by  scalar  multipliers,  which  are  functions 
of  source  energy,  horizontal  distance  and  burst 
height,  in  order  to  obtain  the  correct  free  field  tis¬ 
sue  kerma.  Values  for  all  remaining  kerma  or 
dose  quantities  are  obtained  in  ATR  by  applying 
the  fluence-to-dose  conversion  factors  described 
previously  to  the  adjusted  uniform  air  fluence. 

The  first  correction  factors  placed  in  ATR  to  ac¬ 
count  for  the  presence  of  the  Air-ground  inter¬ 
face  were  based  on  the  Fronch-Mooney  first-last 
collision  estimation  method  (Ref.  24).  normalized 
to  the  two-dimensional  calculations  of  Shaker 
(Ref.  53).  which  pertained  to  a  neutron  source 
at  a  height  of  50  ‘set.  The  next  improvement  was 
based  on  the  extensive  two-dimensional  calcula¬ 
tions  of  Pace  (Ref.  44),  which  provided  data  for 
source  heights  to  300  meters,  although  the 
source  data  were  limited  to  fission  and  fusion 
spectra.  The  most  recent  improvement  prior  to 
the  current  e  fort  used  the  results  of  adjoint 
two-dimensional  calculations  raticed  to  one¬ 
dimensional  results  to  cotam  a  full,  source  ener¬ 
gy-dependent  array  of  ad|ustment  factors  iRef. 
26).  However,  these  data  suffered  from  the  fact 
that  the  ground  used  in  the  two-dimensional  cal¬ 
culations  was  unrealistically  dry.  causing  the  neu¬ 
tron  and  secondary  gamma  ray  dose  immediately 
above  the  ground  to  be  unrealistic  as  well. 

The  air-over-ground  transport  corrections  devel¬ 
oped  for  use  in  thecurrent  update  of  ATR  are  also 
based  on  two-dimensional  adiomt  calculations, 
performed  this  time  with  a  credible  amount 
ground  moisture,  13%  by  weight,  which  is  repre¬ 


sentative  of  that  contained  by  mid-European  and 
United  States  soils  under  moderate  weather  con¬ 
ditions  (Ref.  6).  The  elemental  constituents  and 
density  data  for  the  ground  and  atmosphere  used 
in  the  calculations  are  shown  in  Table  9.  The  at¬ 
mosphere  has  a  graded  density  based  on  a  mod¬ 
el  of  spring  conditions  in  central  Europe,  as 
shown  in  Figure  13. 

The  calculations  were  performed  using  the  DOT 
two-dimensional  discrete  ordinates  transport 
code  (Ref.  25)  with  a  tissue  kerma  adjoint 
source,  the  37  neutron,  21  gamma  ray  group 
DLC-31  cross  section  set  with  P3  legendry  scat¬ 
tering  in  a  Sg  quadrature  set  modified  to  S40  in 
the  polar  direction  to  minimize  streaming  along 
discrete  quadrature  angles  (ray  effects).  Calcu¬ 
lations  were  performed  to  a  usable  source  height 
of  1000  meters  and  a  usable  horizontal  distance 
of  2000  meters.  As  an  additional  precaution 
against  the  effects  of  quadrature  streaming,  the 
uniform  air  calculations  used  to  develop  the  cor¬ 
rection  factors  were  also  performed  using  the 
DOT  code  with  an  identical  atmosphere  profile 
mirrored  at  the  location  of  the  air  ground  inter¬ 
face. 

The  results  of  the  calculations,  as  ratios  of  the 
air-ground  to  uniform  air  kerma  for  a  fission 
source  are  shown  in  Figures  14,  15  and  16  for 
neutrons,  secondary  gamma  rays  and  prompt 
gamma  rays,  respective^-.  Equivalent  data  for  fu¬ 
sion  (12.2  to  14.9  MeV)  neutrons  are  given  for 
neutron  and  secondary  gamma  ray  kerma,  re¬ 
spectively,  m  Figures  17  and  18.  The  data  are 
very  well  behaved  except  within  a  forty-five  de¬ 
gree  cone  below  the  source.  Beyond  that  inner 
region  it  is  clear  that  the  relationship  between 
ground-perturbed  and  uniform  air  transport  is  lin¬ 
ear  with  respect  to  source  height  and  horizontal 
distance.  This  (act  has  been  used  to  extrapolate 
the  results  of  these  calculations  to  greater  dis¬ 
tances.  Within  the  lorty-five  degree  cone,  the  re¬ 
lationship  between  the  uniform  air  and 
air-over-ground  results  is  irregular,  probably  indi¬ 
cating  instabilities  in  the  calculation  of  the  latter 
case.  The  irregularities  do  decrease  with  increas¬ 
ing  burst  height.  Also,  they  were  more  pro¬ 
nounced  when  kerma  ratios  were  taken  with 
respect  to  ANISN  (one-dimensional)  rather  than 
DOT  (two-dimensional)  calculation  results. 
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Table  9.  Air  and  ground  elemental  constituents  (weight  55), 


Dry  Air  Constituents 

Dry 

Ground 

Level  1 

Level  2 

Level  3 

Elements 

Constituents 

0-400m 

400-800m 

800-1100m 

H 

0.14 

C 

1.43 

N 

0.14 

75.72 

75.72 

75.72 

0 

50.63 

24.28 

24.28 

24.28 

Na 

0.79 

Mg 

1.14 

A1 

4.60 

Si 

32.14 

K 

1.65 

Ca 

2.30 

Mn 

0.14 

Fe 

2.25 

Ti 

0.36 

P 

0.12 

Density  (g/cc) 

1.40 

1.155-3* 

1.116-3 

1.074-3 

Free  Moisture 
Density  (g/cc) 

0.21 

6.750-6 

6.148-6 

4.844-6 

Total  Density 
(g/cc) 

1.61 

1.162-3 

1.122-3 

1.079-3 

•Read  as  1.155M0*3. 
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Atmosphere 


Prototype 


Model 


Ground  \ 


Density 

(g  Ice) 

Elevation 

Elevation 

Dry  Air 

Water 

(m) 

(m) 

1.051-3* 

4.03-6 

1500 

1.065-3 

3.94-6 

1407 

1410- 

Atmosphere 
Region  3 

p  dry  Air  1.074-3 

p  water  4.84-6 

1.097-3 

6.45-6 

1086 

iioo- 

Atmosphere 
Reg.on  2 

p  dry  Air  1.116-3 

1.123-3 

5.83-6 

827 

0  water  6.15-6 

1.141-3 

6.35-6 

723 

710- 

* 

Atmosphere 
Region  1 

1.157-3 

6.93-6 

437 

p  dry  Air  1.155-3 
"p  water  6.75-6 

1.181-3 

7.05-6 

310 

3101 

• 

• 

< 

p>* 

r 

Ground 

309  J 

L 

•Read  as  1.051  *10* 3 


Figure  13.  Atmospheric  density  nnd  moisture  profile. 
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The  ratios  of  air-over-ground  to  uniform  air  tissue 
kerma  were  fit  with  a  multi- variable  function, 
including  the  following  terms.  Note  that  all  dis¬ 
tance  quantities  (slant  range,  R.  and  burst  height. 
H)  are  expressed  as  optical  depth  (mass  per  unit 
area) . 

TR'-L  ’EAMS  f QR  FiTT-,0  ;a-OvEa- 

GRGUND  CCRRECTiUN  FAt^^mS 

R  in  quarter  power  steps  from  1  4  to  2 

H  in  quarter  power  steos  from  1 .4  to  2 

(R'H)  in  half  power  steps  from  12  to  2 

9  grazing  angle,  in  radians 

SIN  0 

COS  0 

R'SIN  0 

R-COS  0 

CONSTANT 

Trials  were  run  to  identify  the  most  important  co¬ 
efficients.  fifteen  of  which  were  required  to  repro- 
duce  the  data  everywhere  to  within  five  percent 
or  better.  Those  fifteen  terms  are  as  follows; 


FINAL  TERMS  USED  TO  FIT  AiR-OVER- 
GROUND  CORRECTION  FACTORS 


1 .  H  1/4 

9.  R  2 

2.  H'  2 

10.  0 

3.  H  1 

n  CCS  - 

4.  H  2 

12.  SN  0 

5.  R  1  4 

13.  R'COS  0 

6.  R  1  2 

14.  /P-H)1  2 

7.  R  1 

15.  CONSTANT 

8.  R  3/2 

the  relative  importance  of  these  f  fteen  terms  in 
fittmg  the  data  varies  with  location.  At  short  dis¬ 
tances  (less  than  approximately  300  meters  burst 
height  and  500  meters  horizontal  distance  in  data 
base  air)  the  H  (burst  height)  terms  dominate  the 
value  produced  by  the  function,  Reyond  those 
confines  the  value  .s  dominated  by  the  slant 
range  and  grazing  angle.  At  the  !imts  of  the  data 
the  rslue  's  dominated  by  the  grazing  angle 
a'cne  Thus,  the  extracc'aticn  of  the  correction 
factor  bevond  the  oonf  res  c*  me  data  base  is 
accomyushed  on  the  basis  of  grazing  angle, 
wh .ph  13  used  to  determine  the  applicable  correc¬ 
tion  factor  taken  from  the  data  base  limit. 


SECTION  5 

DELAYED  GAMMA  RADIATION 


ATR  calculates  delayed  radiation  propagation  tak¬ 
ing  into  account  the  time-dependent  fission  prod¬ 
uct  gamma  ray  emissions,  fireball/ weapon  debris 
height  and  air  density  surrounding  the  rising  fire¬ 
ball.  The  fireball  rise  end  air  density  treatments 
in  ATR5  are  based  on  empirical  observations  and 
a  one  dimensional  model,  respectively,  and  are 
unchanged  from  those  used  in  ATR4.  The  bal¬ 
ance  of  the  model  u3ed  in  ATR5  differs  from  that 
used  in  ATR4  as  follows: 

1 .  The  ATR5  source  term  accounts  for  differ¬ 
ences  between  spectra  and  photon  emis¬ 
sion  rates  between  fissionable  isotopes 
U235,  U238  and  Pu239.  whereas  the 
ATR4  model  is  applicable  to  U235  only. 

2.  The  ATR5  delayed  gamma  ray  transport 
model  uses  the  same  data  bese  as  that 
used  for  the  prompt  gamma  rays,  calcu¬ 
lating  scalar  fluence  as  its  primary  quanti¬ 
ty.  whereas  ATR4  has  a  separate  photon 
transport  data  base  consisting  of  a  ker- 
ma-optical  depth  (g/cm2)  relationship. 

3 .  The  ATR5  model  has  been  used  to  calcu¬ 
late  dose  and  dose  rate  values  measured 
at  atmospheric  tests  and  has  been  em¬ 
pirically  adjusted  to  provide  the  best 
agreement  with  available  test  data. 

Delayed  radiation  lluences  are  calculated  in  ATRG 
according  to  the  expression: 

29  18  3 

6,(R)  a  Y  '  FFR  *  £  j  4>|(r.pr.hs)  *  /  (dt  I(FF,) 
kat  1*1  At*  i=»1 

•  (FFY,)  •  S„(t)  (3) 

where  <t>|(R)  is  the  total  gamma  ray  fluence  in 
energy  group  j.  at  a  given  location  R,  rel¬ 
ative  to  the  burst  point, 

Y  is  the  total  device  yield  (user  input) 

FFR  is  the  fraction  of  the  total  yield  attrib¬ 
utable  solely  to  fission  (user  input) 

♦|w(r,pf.hs)  Is  the  gamma  ray  lluence  m 
energy  group  |  per  source  photon  in 
source  group  1 .  determined  according  to 


the  parameters  r,  pr  and  hs  and  the  mid¬ 
point  of  time  increment  k.  where  r  is  the 
source-detector  slant  distance,  pr  is  the 
product  of  r  and  the  average  density  (p) 
over  its  length  and  hs  is  the  source  height 
above  the  ground. 

(FFj)  is  the  number  of  fissions  per  unit 
yield  of  fissile  isotope  i. 

(FFY|)  is  the  fraction  of  the  fission  yield 
due  to  isotope  i.  and 

S*(t>  is  the  time-dependent  source  rate 
in  1  energy  groups  due  to  isotope  i. 

Values  for  quantities  V.  FFR  and  FFY  must  be  in¬ 
put  by  the  user  (see  section  2) .  The  time  incre¬ 
ment  fluence  (<j>p<t).  per  source  group,  is 
calculated  using  the  same  data  base  as  that  for 
prompt  gamma  radiation  transport.  The  parame¬ 
ters  r,  pr,  and  hs  are  calculated  within  ATR5. 
based  on  initial  conditions  of  the  burst,  using 
models  described  in  the  ATR4  Users  Guide  (DNA 
3995F)  (Ref.  33).  ATR5  integrates  the  fission 
product  gamma  ray  dose  rate  over  the  period 
from  zwo  to  sixty  seconds  numerically,  using 
nineteen  time  steps.  The  variation  of  the  dose 
rate  with  time  depends  on  weapon  yield,  burst 
height,  target  height  and  distance.  Therefore,  in 
order  to  employ  so  few  t.me  steps  the  code  itself 
chooses  the  duration  of  each  step.  It  does  so  on 
the  basis  of  the  cumulative  uncollided  fluence. 
which  it  calculates  in  fifty  equal  logarithmic  time 
steps  from  one  tenth  *o  sixty  seconds,  plus  points 
at  0.001  and  0.01  seconds.  The  code  selects  the 
times  at  which  the  uncollioed  fluence  is  closest 
to  (but  larger  than)  the  cumulative  value  in 
one-fifteenth  increments  The  last  time  step  is 
further  divided  into  five  ..mailer  increments,  in 
equal  logarithmic  intervals,  due  to  the  sharp  dose 
rate  gradient  prevalent  at  late  times. 

The  number  of  fissions  per  kiloton  of  fission  yield 
(FFi)  is  stored  m  the  code  and  need  not  be  en¬ 
tered  by  the  user.  The  values  are: 


laotupfl 

Fiasiona/kt 

U235 

1.467F  *  23 

U238 

1.475E  f  23 

Pu239 

1.419E  ♦  23 

4 1 


5.1  FISSION  PRODUCT  GAMMA  RAY  SOURCE. 

The  fission  product  gamma  ray  source  term  Sh ft ) 
is  a  function  of  isotope  <i) .  energy  (l)  and  time 
The  most  comprehensive  measurements  cf  such 
source  spectra  have  been  mcde  by  Fisher  and 
Engle  (Ref,  22).  Dickens,  et  at.  (Ref,  14.  15)  and 
AK,,ama  and  An  -'Ref  2>  The  1962  measure¬ 
ments  included  mean  times  after  last  fission  as 
early  as  0.35  seconds  for  U23r  U238  and  Pu239 
Dickens,  et.al..  measured  da. a  as  early  as  ap¬ 
proximately  two  seconds  after  thermal  fission  of 
U235  and  Pu239.  Akiyama  and  An  measured  fast 
fission  -  related  data  for  several  isotopes  but  at 
relatively  late  times,  a  few  tens  of  seconds  after 
fission. 

The  energy  emission  rates  measured  by  Fisher 
and  Engle  are  not  in  agreement  with  those  of  the 
other  two  groups  for  times  after  fission  at  which 
no  substantial  differences  shoulo  remain  between 
fast  and  thermal  hssion.  This  discrepancy  is  im¬ 
portant  because  only  the  Fisher  and  Engle  data 
provide  fission  product  spectra  applicable  -o  ail 
three  fissile  nuclides  at  times  of  greatest  interest 
Fortunately,  the  differences  between  the  three 
sets  of  data  have  been  reconciled  as  oart  of  the 
!J  S  -  Japan  Join*  Assessment  of  Atomic  Bomb 
Radiation  Dosimetry  in  Hiroshima  and  Nagasaki 
(Ref.  58).  The  Fisher  and  Engle  spectral  data 
have  been  adjusted  based  on  those  of  Dickens, 
e*  at.,  with  a  consistent  application  of  that  adjust¬ 
ment  reproducing  the  Akiyama  and  An  U238  en¬ 
ergy  emission  rate  at  40  seconds.  The  adiusted 
spectra  normalized  to  one  gamma  ray  per  fis¬ 
sion  per  second,  are  given  in  Tables  10.  1  1  and 

Isciaxia 

IJ235  fO-’K  11.42  sec) 

(It  42-,t-  60  sec) 

1)238  (O.t  2.19  sec) 

(2. 19'.t-  0  25  sec) 

(8  25 et-'  15  sec) 

( 1 5^t'  60  sec) 

Pu239  (0-  t'  60  sec) 


12  for  U235,  U238  and  Pu239.  respectively.  The 
development  and  testing  of  these  source  data 
are  described  in  detail  in  Appendix  E  of  this  re- 

por'. 

In  the  process  of  correcting  the  Fishier  and  Engle 
spectra  using  the  data  of  Dickens,  et  a!  .  the  en- 
e>g.  emission  rate  v<as  e*fec  -.ft,  re "mmaiized 
to  that  applicable  to  thermal  re./'en  •  ss-on.  This 
problem  is  alluded  to  in  Apcendix  E,  wnich  de¬ 
scribes  the  correction  process  and  the  subse- 
auent  test-ng  of  the  corrected  spectra  and 
energy  emission  ra,  _s.  Subsequent  to  the  com- 
cietion  cf  the  work  deser  bed  m  Aopendix  E.  En¬ 
gland  and  Schenter  (Ref.  19)  created  a  data  base 
from  a  combination  of  U.S.  (ENDF  B  V  and  VI) 
and  Japanese  (JNDC86)  sources,  consisting  of 
fission  product  inventories  for  each  fissile  nuclide 
and  the  beta  and  gamma  ray  energy  emission 
rates  for  each  fission  product.  These  data  are 
available  for  thermal,  fast  and  14  MeV  neutron 
f-ssion  and  reproduce  the  measurements  of  Dick¬ 
ens.  et  al..  and  Akiyama  and  An  within  a  few  per¬ 
cent.  Therefore,  the  energy  emission  rates  for 
U235.  U238  and  Pu239  adopted  for  use  in  ATR5 
are  those  taken  from  England  and  Schenter  for 
fast  fission  and  fit  to  the  functional  form: 

Energy  Emission  Pate  =  A 

iTat 

(MeV  per  fission  per  second'  (4) 

where  t  is  the  t.me  in  seconds.  Tne  coefficients 
used  m  the  emission  rate  model  are: 


Coeff'ciont 

JL 

JL 

0  606 

0.833 

4  818 

7  234 

1  191 

1  039 

2  40Q 

2  1)69 

0  8  28 

-t .  105 

•0.44? 

0.589 

0  340 

0  618 
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Table  10.  Fission  product  gamma  ray  source  spectra  for  U235  at  mean  times 
aftei*  fission  of  0.35,  1.5,  4.75,  11.5  and  40  seconds. 


Energy  Upper  Gamma  rays/MeV-Gamma  ray 

Group  Energy  Time  (sec)  : 

No.  (MeV)  0.35  1.50  4.75  11.50  40.00 

1  0.045a  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0.000E+00  0.000E+00 

2  0.100  O.OOOE+OO  0.000E+00  C.OOOE+OO  O.OOOE+OO  O.OOoE+OO 

3  0.150  4.933E-01b  3.675E-01  4.307E-01  4.210E-01  3.769E-01 

4  0.300  1.038E+00  8.693E-01  9.983E-01  1.010E+00  9.071E-01 

5  0.450  6 . 942E-01  6.720E-01  7.311E-01  7.382E-01  7.324E-01 

6  0.700  6.764E-01  7.565E-01  8.226E-01  7.801E-01  6.182E-01 

7  1.000  5.391E-01  5.263E-01  4.511E-01  4.589E-01  4.965E-01 

8  1.500  3 . 230E-01  3.388E-01  3.050E-01  3.252E-01  3.648E-01 

9  2.000  1.G05E-01  1.715E-01  1.728E-01  1.664E-01  1.868E-01 

10  2.500  1.066E-01  1.100E-01  1.081E-01  1.065E-01  1.209E-01 

11  3.000  6.688E-02  7.002E-02  6.091E-02  6.624E-02  8.265E-02 

12  4.000  3.840E-02  4.177E-02  3.993E-02  3.671E-02  3.763E-02 

13  5.000  1.396E-02  1.189E-02  1.035E-02  1.091E-02  1.230E-02 

14  6.000  2.797E-03  3.618E-03  3.420E-03  3.571E-03  3.628E-03 

15  7.000  9.842E-04  8.458E-04  7.598E-04  5.387E-04  4.233C-04 

16  8.000  1.399E-04  1.400E-04  1.396E-04  1.399E-04  1.398E-04 

17  10.000  0.000E+00  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

18  12.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 


aLower  energy  bound  0.010  MeV 
bRead  as  4.933*10*1 


Table  11  Fission  product  gamma  ray  source  spectra  for  U238  at  mean  times 
fission  of  0.35,  1.5,  4.75,  11.5  and  40  seconds. 


Energy 

Group 

No. 

Upper 

Energy 

(MeV) 

Time  (sec) 
0.35 

Gamma  i  aj  s/MeV-Gamma  ray 

1.50  4.75  11.50 

40.00 

1 

0.045a  0.000E+00 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

2 

0.100 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

3 

0.150 

5.568E-0lb  3.742E-01 

5.271E-01 

4 . 505E-01 

3.881E-01 

4 

0.300 

1 . 108E+00 

9.559E-01 

1 . 125E+00 

1.049E+00 

1.009E+00 

5 

0.450 

7.432E-01 

7. 151E-01 

7.608E-01 

8.385E-01 

7.773E-01 

6 

0.700 

6.223E-01 

8.674E-01 

7.825E-01 

7.532E-01 

6.814E-01 

7 

1.000 

4.881E-01 

4.934E-01 

4.347E-01 

4.335E-01 

4.891E-01 

8 

1.500 

3.039E-01 

3.313E-01 

2.826E-01 

2.956E-01 

3.311E-01 

9 

2.000 

1.564E-01 

1.852E-01 

1.685E-01 

1.663E-01 

1.605E-01 

10 

2.500 

1.253E-01 

1 . 151E-01 

1.076E-01 

1 . 121E-01 

1 .21  IE-01 

11 

3.000 

8. 130E-02 

8.057E-02 

6.249E-02 

6.765E-02 

8.479E-02 

12 

4.000 

3.971E-02 

4.338E-02 

4.011E-02 

3 . 923E-02 

3 . 392E-G2 

13 

5.000 

1.387E-02 

1 . 180E-02 

9.619E-03 

1 . 124E-02 

9 . 236E-03 

14 

6.000 

3 . 882E-03 

3.409E-03 

3. 587E-03 

3.951E-03 

3. 109E-03 

15 

7.000 

1.409E-03 

8.844E-04 

7.119E-04 

8. 154E-04 

3.538E-CM 

16 

8.000 

1.398E-04 

1.399E-04 

1.402E-04 

1.398E-04 

1.403E-04 

17 

10.000 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

IS 

12.000 

0.00OE+00 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

0.000E+C0 

aLower  energy  bound  0.010  MeV 
*>Read  as  5.568*10-1 
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Table  12.  Fission  product  gamma  ray  source  spectra  for  U239  at  mean  times 
after  fission  of  0.35,  1.5,  4.75,  11.5  and  40  seconds. 


Energy 

Group 

No. 

tpper 

Energy 

(MeV) 

Time  (sec) 
0.35 

Gamma  rays/MeV-Gamma  ray 

1.50  4.75  11.50 

40.00 

1 

0.045a  0. 000E+00 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

2 

0.100 

0 . OOOE+OO 

0. 0O0E+00 

0. OOOE+OO 

0  OOOE+OO 

0. OOOE+OO 

3 

0.150 

5.425E-01*>  3.966E-01 

4.7G9E-01 

4.377E-01 

3.764E-01 

4 

0.300 

1 . 094E+00 

9.689E-01 

1 . 107E+00 

1 . 145E+00 

9.238E-01 

5 

0.450 

7 . 172E-01 

7 . 5G0E-01 

7:595E-01 

8. 111E-01 

7.964E-01 

0 

0.700 

7.448E-01 

8.944E-01 

9.002E-01 

8 . 086E-01 

0.484E-O1 

7 

1.000 

5.210E-01 

4.882E-01 

4.412E-01 

4.426E-01 

4.859E-01 

8 

1.500 

2.982E-01 

3. 194E-01 

2.792E-01 

2 . 967E-01 

3.430E-01 

9 

2.000 

1.490E-01 

1.496E-01 

1.098E-O1 

1.G79E-01 

1.839E-01 

10 

2.50C 

1.014E-01 

9.263E-02 

9 . 549E-02 

9 . 96GE-02 

1.300E-01 

11 

3.000 

7.294E-02 

5.896E-02 

5 . 283E-02 

5 . 935E-02 

8.725E-02 

12 

4.001 

3.249E-02 

3.024E-02 

3.11GE-02 

2 . 809E-O2 

3 . 345E-02 

13 

5.000 

1.250E-02 

9.317E-03 

6 . 433E-03 

6.850E-03 

7.245E-03 

14 

0.000 

2. 182E-03 

1 . 805E-03 

1 . 888E-03 

1 . 850E-03 

2.019E-03 

15 

7.000 

5.915E-04 

5.696E-04 

5.417E-04 

3.981E-04 

3.390E-O4 

10 

8.000 

1.398E-04 

1 .400E-04 

1.396E-04 

1 . 397E-04 

1.403E-C4 

17 

10.000 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

18 

12.000 

0. OOOE+OO 

O.OOCE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

aLower  energy  bound  0.010  MeV, 
“Read  as  5. 425* 10“ 1 
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These  coefficients  reproduce  the  England  and 
Schenter  gamma  ray  energy  enTssion  rate  data 
within  two  percent  or  less  over  the  time  periods 
indicated. 

The  reported  emission  rate  spectra  have  been 
assigned  to  time  regimes,  the  limits  of  which  are 
the  time  increment  boundaries  closest  to  the 
mean  of  the  logarithm  of  adjacent  reported 
times,  as  follows: 

Time  (sec)  of 

Iima..iaesJ  Applicable  Spectrum 
0  -  0.075  0.35 

0.75  -  3.25  1.50 


ATR5  model  lay  in  its  inability  to  account  for  the 
ground-reflected  shock  wave  from  the  burst.  The 
ground-reflected  shock  wave  fills  in  the  low  densi¬ 
ty  region  which  includes  the  rising  fireball,  provid¬ 
ing  more  mass  for  the  fission  product  gamma 
rays  to  penetrate  and.  consequently,  a  lower 
dose  rate,  as  described  in  Appendix  E  of  this  re¬ 
port.  This  effect  is  particularly  well  defined  for  a 
burst  at  110m  scaled  height,  as  shown  in  Figure 
20.  The  two-dimensional  model  begins  to  fill  in 
the  line  of  sight  from  source  to  detector  at  ap¬ 
proximately  1  second,  similar  to  the  effect  shown 
in  the  measured  data,  though  the  two-dimension¬ 
al  model  ever-estimates  the  effect.  The  data 
from  ATR5,  based  on  a  one-dimensional  hydro- 
dynamic  model,  show  no  such  effect  and  over¬ 
predict  the  dose  rate  at  later  times. 


3.25  - 

7.75 

4.75 

7.75  - 

22.5 

11.50 

22.5  - 

60.0 

40.00 

5.2  EMPIRICAL  CORRECTIONS. 

The  ATR5  model  was  used  to  calculate  delayed 
radiation  as  measured  at  several  atmospheric 
weapon  tests.  These  included  yields  from  several 
kilotons  to  several  megatons  and  scaled  burst 
heights  (HOB(m) /yield  to  the  one-third  power) 
from  zero  to  over  one  hundred  meters.  The  mea¬ 
surements  included  time-dependent  data  from 
active  dosimeters  and  total  dose  data  from  pas¬ 
sive  dosimeters,  such  as  film  badges.  The  film 
badge  data  were  used  only  when  the  vield  was 
sufficiently  large  that  delayed  radiation  compo¬ 
nent  was  dominant  or  when  the  calculated 
prompt  and  secondary  components  were  con¬ 
firmed  by  time  dependent  measurements,  and 
could  be  subtracted  from  the  total  to  obtain  the 
delayed  component. 

The  calculations  and  measurements  for  each 
shot  were  compared  over  all  distances  and  the 
mean  of  all  comparisons  for  each  shot  was  com¬ 
pared  with  those  for  all  the  other  shots.  Upon 
comparing  the  calculations  and  measurements 
as  a  function  of  scaled  height  of  burst,  as  shown 
in  Figure  19,  It  was  evident  that  the  ATR5  model 
deviated  from  the  measured  data  as  a  strong 
function  of  the  scaled  burst  height.  Further  com¬ 
parisons  against  the  results  of  a  model  which 
containrd  the  same  source  and  transport  com¬ 
ponents  as  ATR5  but  included  a  more  sophisti¬ 
cated,  two-dimensional  hydrodynamic  model 
indicated  that  the  primary  shortcoming  of  the 


Unfortunately,  reliable  delayed  gamma  ray  dose 
rate  data  are  not  available  for  scaled  burst 
heights  greater  than  1 10.4  meters.  However,  re¬ 
sults  of  the  two-dimensional  delayed  radiation 
model  combined  with  detailed  calculations  of  the 
prompt  and  secondary  radiation  components 
compare  well  with  time-integral  measurements  of 
total  initial  gamma  ray  dose  recorded  for  bursts 
having  scaled  burst  height  as  large  as  182.3  me¬ 
ters  (Nagasaki)  and  235.2  meters  (Hiroshima). 
Thus,  the  two-dimensional  model  was  used  as 
tne  standard  by  which  to  assess  the  deviation  of 
the  one-dimensional  model  at  scaled  burst  height 
greater  than  1 10.4  meters.  As  can  be  seen  from 
Figure  19  the  one-  and  two-dimensional  models 
come  into  agreement  at  a  scaled  burst  height  of 
235.2  meters.  It  is  not  surprising  that  this  should 
occur,  since  at  such  large  heights  the  reflected 
shock  wave  has  little  or  no  effect  on  the  rising 
fireball. 

Comparisons  of  experimental  and  two-dimen¬ 
sional  model  dose  values  with  those  produced 
using  the  one-dimensional  model  ATR5  were 
used  to  produce  an  empirical  correction  factor 
based  on  scaled  burst  height.  Linear  interpolation 
is  used  to  obtain  correction  factor  values  be¬ 
tween  break  point  values,  as  shown  below. 
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two-diaentional  aodel,  as  related  to  scaled  bexg] 


flUEN/SEC  **  TIME 


:d  Shockwave 
.ine  of  Sight 
Seconds 


L  i.  i  i 


S) 


10 


10 


delayed  radiation  exposure 
ed  hurst  height  event  taken 
tance . 


The  correction  factors  for  zero  and  17.1  meter 
scaled  burst  height  were  set  to  unity  because  the 
comparisons  of  ATR5  with  available  data  pro¬ 
duced  very  large  uncertainty  bands  at  those  burst 
heights,  which  included  the  unity  value.  Thus,  the 
data  did  not  conclusively  indicate  that  the  prob¬ 
lem  of  reflected  shock  for  very  low  burst  heights 
produced  any  effect  beyond  that  already  ac¬ 
counted  for  in  the  model.  (For  very  low  burst 
heights  the  one-dimensional  model  includes  a 
provision  for  increasing  the  yield  for  the  purpose 
of  calculating  the  hydrodynamic  perturbation  of 
the  surrounding  air). 

Even  with  the  scaled  burst  height-based  empirical 
correction,  problems  remain  in  ATR,  oroisably 
having  to  do  with  other  limitations  of  the  hydrody¬ 
namic  model,  such  as  its  inability  to  desense  the 
movement  of  the  gamma  ray  source  within  the 
fireball  or  the  formation  of  the  toru3  cloud  at  late 


times,  as  described  in  Appendix  E.  These  prob¬ 
lems  are  evidenced  by  rarge-dependent  discre¬ 
pancies  between  calculated  and  measured  dose 
values.  These  are  particularly  evident  for  surface 
bursts  and  very  low  air  burst  ,3hOB  <  50  meters) . 
with  the  corrected  (SHOB  effect)  model  underes¬ 
timating  the  measured  values  by  twenty  to  thirty 
percent  within  a  few  fireball  radii  and  over  esti¬ 
mating  them  by  a  similar  amount  at  distances  of 
many  fireball  radii.  Unfortunately,  there  is  not  a 
strong  enough  correlation  with  fireball  radius  or 
any  othsr  readily  observable  parameter  to  create 
a  scalar  correction  factor  to  minimize  the  error 
in  a  systematic  and  reliable  way.  It  is  recom¬ 
mended  that  future  research  be  directed  toward 
developing  a  reliable  high  order  model  of  the  phe¬ 
nomena  affecting  delayed  radiation  transmission 
and  that  such  a  model  be  used  as  the  basis  for 
deriving  further  improvements  for  ATR,  since  the 
available  measured  data  are  very  limited. 
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SECTION  6 


SAMPLE  CALCULATIONS 


This  section  contains  samples  of  calculations 
which  exemplify  the  new  capr-c'ities  added  to 
ATR5.  The  input  and  output  for  each  calculation 
are  described  in  detail,  particularly  as  regards 
ATP5  modifications.  Additional  explanations  of 
ATR  input  and  output  can  be  obtained  from  the 
documenration  for  ATR4  (DNA  3995F). 

6.1  SAMPLE  CALCULATION  1.  USE  OF  THE 
•MET  COMMAND. 

This  calculation  illustrates  the  use  of  two  varia¬ 
tions  of  the  'MET  command  to  obtain  the  same 
neutron  values.  In  the  first  entry  the  'MET  com¬ 
mand  is  used  to  specify  the  Temperature  (20  c) . 
Pressure  (1000  mb)  and  relative  humidity  (90%). 


After  an  'EXC  command,  which  causes  the  cal¬ 
culation  to  be  executed,  the  'MET  command  is 
re-entered,  proviomg  data  whicn  supersedes  that 
contained  in  the  previous  calculations  but  leaving 
all  other  commands  unchanged.  The  revised 
'MET  entry  provides  a  moist  air  density  (1.179 
E-3  g/oc)  and  a  percent  moisture  by  weight 
(1.3043),  which  are  identical  to  those  obtained 
using  the  data  from  the  first  'MET  entry.  The  re¬ 
sults  of  the  two  entries  are  virtually  identical,  as 
should  be  expected.  Note  that  ATR5  calculates 
a  ground  level  from  the  data  provided  in  the  'MET 
command.  That  ground  level  will  supersede  any 
which  is  entered  by  the 'GROUND  command  un¬ 
less  no  'MET  input  i3  provided. 


SAMPLE  PROBLEM  1  -  INPUT 


'ATR  VERSION  5  -  -  -  CAIC  (1  R  1989  REVISION) 
'TITLE  SAMPLE  PROBLEM  1  -  AT  ERSION  5 
'N-SOURCE(I) 

'N-YIELD  10 
'N-NORM  2.0E  +  23 
'HS.M.  100 
'HT.M, 1 

' MET ( 1 )  20  1000  90 
•RH.M.  500(500)2000 
'DOSE/N/ 

'EXC 

*MET  (3)  1.179E-3  1.3043 

•EXC 

•STOP 

•FIN 
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SAMPLE  PROBLEM  1  -  OUTPUT 


ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  1  -  ATR  VERSION  5 

*********************************** 


NEUTRON  SOURCE  INTERNAL  FISSION 
NORMALIZATION-2. 000E+23  NEUTRON  /KT,  YIELD=1 .OOOE+Ol  KT 
TOTAL  OUTPUT=2 . OOOE+24  NEUTRON 

SOURCE  SPECTRUM 

ENERGY  (MEV  )  N  N/MEV  ENERGY  (MEV  )  N  N/MEV 

1 . 07E-05-2 . 9GE-05  O.OOE+OO  O.OOE+OO  1 .83E+00-2.31E+00  2.03E+23  4.24E+23 
2 . 90E-05- 1.01 E-04  O.OOE+OO  O.OOE+OO  2.31E+00-3.01E+00  7.74E+22  1.11E+23 
1 .01E-04-1 .23E-03  O.OOE+OO  O.OOE+OO  3.01E+00-4.07E+00  1.10E+23  1.03E+23 
1 . 23E-03-2 . 19E-02  3.30E+22  1.60E+24  4.07E+00-4.97E+00  2.35E+22  2.62E+22 
2 . 1 9E-02- 1 . 1 IE-01  4.12E+23  4.63E+24  4.97E+00-6.38E+00  3.66E+22  2.60E+22 
1 . 11E-01-1 .58E-01  3.60E+22  7.66E+23  6.38E+00-8.19E+00  2.55E+22  1.41E+22 
1 . 58E-01 -5 . 50E-01  3.03E+23  7.72E+23  8.13E+00-1 .OOE+Ol  1.47E+22  8.11E+21 
5.50E-01-1 . 11E+00  4.32E+23  7.71E+23  1 .00E+01-1 .22E+01  O.OOE+OO  O.OOE+OO 
1 . 11E+00-1 .83E+00  2.94E+23  4.08E+23  1 .22E+01-1 .49E+01  O.OOE+OO  O.OOE+OO 

*********************************** 

ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  I  -  ATR  VERSION  5 


GROUND  LEVEL 

.500KM, 

59.819GM/CM**2, 

1.041KFT, 

.311MILES 

HORIZ.  RANGE  RH- 

.500KM, 

57.782GM/CM**2, 

1.640KFT, 

.311MILES 

•SLANT  RANGE 

RS- 

.510KM, 

60.096GM/CM**2, 

1.072KFT, 

.317MILES 

TARGET  ALT. 

HT- 

.001KM, 

. 117GM/CM**2, 

.003KFT, 

.001MILES 

SOURCE  ALT. 

HS- 

. 100KM, 

11.613CM/CM**2, 

.328KFT, 

.062MILES 

•SLANT  ANGLE 

AN- 

-11. 200DEGREES  (COS-  .98096) 

•CALCULATED  FROM  OTHER  COORDINATES 
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NEUTRON 


) 


DOSE  VS.  HORIZ.  RANGE  (FLADS 
DOSE  1=  ANS  STANDARD- REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3=  MID-PHANTOM 
DOSE  4=  CONCRETE 
DOSE  5=  AIR 

DOSE  6=  NON- IONIZING  SI. 

DOSE  7=*  IONIZING  SILICON 
DOSE  8=  MID-HEAD 
DOSE  9=  1MEV  EQ.  FLUENCE 
DOSEIO-  USER  SUPPLIED 

HORIZ.  RANGE  DOSE  1  DOSE  2  DOSE  3  DOSE  4  DOSE  5  DOSE  6 

5 . C00E+02  M  4.58E+06  3.44E+04  i.l8E+04  3.53E+03  5.41E+03  4.355+02 

1 . 000E+03  M  7.42E+03  5.P1E+02  2.08E+02  5.82E+01  9.83E+01  7.06E+00 

1 . 500E+03  M  2.12E+02  1.65E+01  0.O3E+OO  1.75E+00  2.91E+00  2.10E-01 

2 . OOCE+03  M  9.21E+00  7.34E-01  2.59E-01  7.94E-02  ] .27E-01  9.43E-03 

DOSE  7  DOSE  8  DOSE  9  D0SE10 
4.04E+02  1.59E+04  1.32E+13  O.OOE+OO 
7.05E+00  2.79E+02  2.14E+11  O.OCE+OO 
2.28E-01  8.17E+00  0.34E+O9  O.OOE+OO 
1.07E-02  3.55E-01  2.85E+08  O.OOE+OO 


* ‘EXECUTION  COMIT.ETED 


8.2  SAMPLE  CALCULATION  2.  USER- 
SUPPLIED  'FLUXWT. 

A  gamma-ray  source  is  used  to  illustrate 
'FLUXWT  input  and  output  in  ATR5.  The 
"FLUXWT  command  allows  the  user  to  provide 
his  own  fluence-to-dose  factors.  These  are  en¬ 


tered  after  the  units  are  given  (up  to  16  charac¬ 
ters  not  including  the  bracketing  commas)  in 
energy  order  low  to  high.  In  ATR5,  the  resuiting 
dose  values  are  printed  as  part  of  the  'DOSE 
command  in  position  8  for  gamma  rays  and  posi¬ 
tion  10  for  neutrons. 


SAMPLE  PROBLEM  2  -  INPUT 

■ATR  VERSION  5 - SAIC  (1  APR  1£89  REVISION) 

■TITLE  SAMPLE  PROBLEM  2  -  ATR  VERSION  5 

*  G-SOURCE  ( 1 ) 

■G-YIELD  10 

*  G-NCRM  7.5E+22 

*  GROUND,  M. 500 
•HS.M.50 
•HT.M.1 

•RH.M,  500(500)2000 
■DOSE/G/ 

'  FLUXWT/G/, ROENTGENS. 7. 185E-1 1  3.616E-11  5.464E-11  1.079E-10 
1.993E-10  3.066E-10  4.433E-10  6.054E-10  7.820E-10  9.06SE-10 
1 . 1 10E-09  1.255E-09  1.487E-09  1.707E-09  1.890E-09  2.192E-09 


2.448E-09  2.823E-09 


•EXC 

‘STOP 

■FIN 
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SAMPLE  PROBLEM  2  -  OUTPUT 


ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  2  -  ATR  VERSION  5 

I****.****************************** 

GAMMA  SOURCE  INTERNAL  FISSION 
NORMALIZATION-7.SOOE+22  GAMMA  /KT,  YIELD=1 .OOOE+Ol  KT 
TOTAL  OUTPUTS . 500E+23  GAMMA 

SOURCE  SPECTRUM 

ENERGY  (MEV  )  G  G/MEV  ENERGY  (MEV  )  G  G/MEV 

1 . 00E-02-4 . 50E-02  2.44E+22  0.98E+23  2.00E+00-2.60E+00  2.81E+22  5.81E+22 
4 . S0E-02- 1 . OOE-O 1  1.23E+22  2.24E+23  2. 50E+00-3. OOE+OO  1.07E+22  3.34E+22 
1 .OOE-Ol-1 .50E-01  3.G0E+22  7.32E+23  3. 00E+00-4. OOE+OO  1.58E+22  1.58E+22 
1 . 50E-01 -3 . OOE-Ol  7 .74E+22  5.10E+23  4 . 00E+00-5 . OOE+OO  5.59E+21  5.59E+21 
3.00E-01-4.50E-01  1.02E+23  8.79E+23  5. OOE+OO-0. OOE+OO  1.99E+21  1.99E+21 
4 . 50E-01 -7 . OOE-Ol  1.52E+23  6.08E+23  6. OOE+OO-7. OOE+OO  8.94E+20  0.94E+2O 
7. OOE-Ol -1 .OOE+OO  1 .22E+23  4.08E+23  7. 00E+00-8. OOE+OO  2.83E+20  2.83E+20 
1 . OOE+OO- 1 .50E+00  1 .06E+23  2.11E+23  8. OOE+OO- 1 .OOE+Ol  1.18E+20  5.90E+19 
1 . 50E+00-2 . OOE+OO  4.82E+22  9.04E+22  1 .OOE+Ol-l .20E+01  0. OOE+OO  0. OOE+OO 


ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  2  -  ATR  VERSION  5 


GROUND  LEVEL 

. 500KM , 

69.787GM/CM**2, 

1.640KFT, 

.31 1MILES 

HORIZ.  RANGE  R1I- 

.600 KM, 

58.O06GM/CM**2, 

1.640KFT, 

.31 1MILES 

•SLANT  RANGE 

RS- 

.602 KM, 

58.482GM/CM**2, 

1.048KFT, 

.312MILES 

TARGET  ALT. 

HT- 

.001  KM, 

.117GM/CM**2, 

.003KFT, 

.001MILES 

SOURCE  ALT. 

HS- 

.060 KM, 

6. 821  CM/CM* *2, 

. 184KFT, 

.031MILES 

•SLANT  ANGLE 

AN- 

-6 . 697DEGREES  (COS-  .996: 

>3) 

•CALCULATED  FROM  OTHER  COORDINATES 


GAMMA  DOSE  VS.  HORIZ.  RANGE  (RADS  ) 

DOSE  1=  ANS  STANDARD-REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3=  MID-PHANTOM 
DOSE  4=  CONCRETE 
DOSE  5=  AIR 

DOSE  6=  IONIZING  SILICON 
DOSE  7=  MID-HEAD 
DOSE  8=  USER  SUPPLIED 

HORIZ.  RANGE  DOSE  1  DOSE  2  DOSE  3  DOSE  4  DOSE  5  DOSE  6 

5 . 000E+02  M  1 . 72E+03  1 . 30E+03  7 . 95E+02  1 . 71E+03  1 . 18E+03  1 . 83E+03 

1 .000E+03  M  4.55E+01  3.64E+01  2.32E+01  4.34E+01  3.31E+01  4  57E+01 

1.500E+03  M  2.90E+00  2.41E+00  1.57E+00  2.69E+00  2.18E+00  2.81E+00 

2.000E+03  M  2.83E-01  2.40E-01  1.59E-01  2.61E-01  2.18E-01  2.71E-01 

DOSE  7  DOSE  8 
9 . 25E+02  1 . 36E+03 
2.64E+01  3.80E+01 
1 . 77E+00  2.51E>00 
1.78E-01  2.50E-01 


** EXECUTION  COMPLETED 
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8.3  SAMPLE  CALCULATION  3.  USE  OF  THE 
*FP  COMMANDS. 

ATR5  allows  the  user  to  specify  the  fractions  of 
tiss.on  yield  due  to  the  three  primary  fissile  mate¬ 
rials  U-235.  U-238  and  Pu-239.  These  fractions 
are  input  via  the  * FP-t  command  and  are  used 

pr  iv  m  the  'airiuiaton  of  the  ''ssion  product 
ideiaved)  component.  Fatiuie  to  pmvae  a  ‘FO-i 


command  when  calculating  delayed  radiation  re¬ 
sults  in  U-235  fission  by  default.  In  the  case  of 
the  sample  problem,  a  hypothetical  weapon  hav¬ 
ing  a  total  yield  of  100  kt.  of  which  95%  is  due 
to  f  ssion.  is  specified  using  the  "FP-Y  and  "FP-F 
commands,  respectively.  The  "FP-I  command 
v  oi.  -ves  that  95%  of  the  yield  is  aue  to  U-235. 
c.%  s  hue  to  U-238  and  no  yeid  is  cenved  f'om 


SA‘-;°LE  PROBLEM  3  -  INPUT 


•ATH  VERSION  5  -  -  -  SAlC  (1  APR  1989  REVISION) 


•TITLE  SAMPLE  PROBLEM  3  -  ATR  VEPSiON  5 

•FP-Y  100 

'FP-F  0.95 

•FP-I  0.95  0.05  0.0 

■HS.M.250 

•HT.M,  1 

*  GROUND.  M.  500 
•RH.M.FOO 

•RH.M,  1000(500)2500 
•DOSc  FP/ 

•EXC 

•STOP 

•FIN 
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SAMPLE  PROBLEM  3  -  OUTPUT 


ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  3  -  ATR  VERSION  5 

*********************************** 


FISSION  PRODUCT  YIELD=1 . 00E+02XT ,  FISSION  FRACTION*  9.50E-01 

*********************************** 

ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  3  -  ATR  VERSION  5 


*********************************** 


GROUND  LEVEL 

.500KM, 

59.787GM/CM**2, 

1.640KFT, 

.311MILES 

HORIZ.  RANGE 

RH- 

l.OOOKM, 

113.897GM/CM“2, 

3.281KFT, 

.621MILES 

'SLANT  RANGE 

RS= 

1.031KM, 

118.805GM/CM“2f 

3.381KFT, 

. 640MILES 

TARGET  ALT. 

HT- 

.OC1KM, 

.117GM/CM“2, 

.003KFT, 

.001MILES 

SOURCE  ALT. 

HS= 

.250KM, 

28.823GM/CM“2, 

.820KFT, 

. 155MILES 

'SLANT  ANGLE 

AN= 

-13.982DEGREES  (COS-  .97037) 

‘CALCULATED  FROM  OTHER  COORDINATES 

i  *  *  *  *  *  * 

*  * 

*  *  *  *  * 

********* 

****** 

*  *  *  *  * 

FISSION  PRODUCT  DOSE  VS 

.  HORIZ.  RANGE 

(RADS 

) 

DOSE  1= 

ANS  STANDARD-REM 

DOSE  2* 

SOFT  TISSUE 

DOSE  3= 

MID-PHANTOM 

DOSE  4- 

CONCRETE 

DOSE  5» 

AIR 

DOSE  6- 

IONIZING  SILICON 

DOSE  7- 

MID-HEAD 

DOSE  8- 

USER  SUPPLIED 
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HORIZ.  RANGE  DOSE  1  DOSE  2 
1.000E+03  M  *  1 . 04E+04  7 . 99E+03 
1 .500E+03  M  8 . 92E+02  7.12E+02 

2 . OOOE+03  M  9. 15E+01  7.54E+01 

2 . 500E+03  M  1 . 11E+01  9.31E+00 

DOSE  7  DOSE  8 
5.72E+03  0 . OOE+OO 
5 . 15E+02  0 . OOE+OO 
5.51E+01  0. OOE+OO 
6.85E+00  0. OOE+OO 


DOSE  3  DOSE  4  DOSE  5  DOSE  6 
4.95E+03  1.03E+04  7.27E+03  1.09E+04 
4.52E+02  8.52E+02  6.46E+02  8.96E+02 
4.88E+01  8.51E+01  6.83E+01  8.88E+01 
6.12E+00  1.02E+01  8.43E+00  1.06E+01 


******************  +  ***  +  ****  +  ******* 

“EXECUTION  COMPLETED 
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8.4  SAMPLE  CALCULATION  4,  CALCULATING 
TOTAL  DOSE  (‘DOSE/T/). 

ATR5  has  completely  integrated  the  fission  prod¬ 
uct  (delayed)  radiation  module  into  the  code  for 
the  purpose  of  calculating  dose,  as  specified  us¬ 
ing  the  'DOSE  command.  The  commands: 

•DOSE/N/ 

•DCSE/G/ 

*  DOSE/NG/ 

•DOSE/FP/ 


provide  only  the  values  computed  for  each  com¬ 
ponent.  respectively.  The  command  'DOSE/GG/ 
provides  the  total  gamma-ray  dose  for  all  sources 
specified  including  fission  production  radiation, 
regardless  of  whether  the  individual  dose  compo¬ 
nents  are  requested.  The  'DOSE/GG/  command 
should  not  be  used  when  only  on  source  of  gam¬ 
ma  rays  (including  neutrons)  is  stipulated.  Like¬ 
wise,  the  ’DOSE/T/  command  provides  a  total  of 
all  components  resulting  from  the  sources  stipu¬ 
lated,  including  fission  product  radiation,  regard¬ 
less  of  whether  the  dose  components  are 
requested  individually. 


SAMPLE  PROBLEM  4  -  INPUT 

*ATR  VERSION  5  -  -  -  SAIC  (1  APR  1989  REVISION) 


•TITLE  SAMPLE  PROBLEM  4  -  ATR  VERSION  5 
•N-SOURCE(2) 

•N-YIELD  200 
’N-NORM  2.0E+23 

*  G-SOURCE  ( 1 ) 

'G-YIELD  200 
'G-NORM  7.5E+22 
*FP-Y  200 

•FP-F  0.9 
•FP-I  0.0  0.2  0.8 

•  GROUND,  M,  500 
•HS.M.300 
*HT,M,1 

•RH,M,  1500(500)3000 
'DOSE/N/ 

•DOSE/G/ 

•DOSE/NG/ 

•DOSE/FP/ 

•DOSE/GG/ 

•DOSE/T/ 

•EXC 

'STOP 

'FIN 
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SAMPLE  PROBLEM  4  -  OUTPUT 


ATP.  PROBLEM  NUMBER  1 


SAMPLE  PROBLEM  4  -  AIR  VERSION  5 


************** 


******************* 


NEUTRON  SOURCE 


INTERNAL  THERMONUCLEAR 


NORMALIZATION-2. 000E*23  NEUTRON  /KT,  YIELD=2.000E+02  KT 
TOTAL  0UTPUT=4 . OOOE+25  NEUTRON 


ENERGY (MEV  ) 

1 . 07E-05-2 . 90E-05 
2.90E-05-1.01E-04 
1 .01E-04-1 .23E-03 
1 . 23E-03-2 . 19E-02 


SOURCE  SPECTRUM 
N  N/MEV  ENERGY (MEV  ) 

O.OOE+OO  O.OOE+OO  1 .83E+00-2.31E+00 
8.00E+22  1.11E+27  2.31E+00-3.01E+00 
2.29E+24  2.03E+27  3.01E+00-4.07E+00 
1.38E+25  6.66E+26  4.07E+00-4.97E+00 


1.04E+24 


N/MEV 
2. 16E+24 


1.04E+24  1 .49E+24 
1 .04E+24  9.81E+23 


2. 19E-02-1 . 11E-01  4 . 39E+24  4.92E+25  4 . 97E+00-8 . 38E+00 


6 . 80E+23 
7 . 20E+23 


1 . 11E-01-1 .58E-01  4 . 57E+23  9.72E+24  6.38E+00-8.19E+00  5.88E+23 


1.58E-01-5.50E-01 
5.50E-01-1 . 11E+00 


3.62E+24  9.24E+24  8.19E+00-1.00E+01  5.64E+23 
3.40E+24  0 .07E+24  1 .OOE+Ol-1 .22E->01  1.02E+24 


7.56E+23 
5. 11E+23 
3.25E+23 
3. 12E+23 
4.65E+23 


1 . 1 1E+00-1 . 83E+00  2.48E+24  3.44E+24  1 .22E+01-1 .49E+01  2.82E+24  1.05E+24 


60 


*********************************** 


ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  4  -  ATR  VERSION  5 

*********************************** 

GAMMA  SOURCE  INTERNAL  FISSION 
NORMALIZATION-7. S00E+22  GAMMA  /KT,  YIELD=2 . 0O0E+02  KT 
TOTAL  OirTPUT=l .  500E+25  GAMMA 

SOURCE  SPECTRUM 

ENERGY (MEV  )  G  G/MEV  ENERGY (MEV  )  G  G/MEV 

1 . 00E-02-4 . 50E -02  4.89E+23  1.40E+25  2.00E+00-2.50E+00  5.61E+23  1J.2E+24 
4 . 50E-02- 1 . 00E-01  2.47E+23  4.48E+24  2. 50E+00-3. OOE+OO  3.34E+23  6.68E+23 
1 .OOE-OI-l .50E-01  7.32E+23  1.48E+25  3. OOE+OO-4. OOE+OO  3.16E+23  3.16E+23 
1.50E-01-3.00E-01  1.55E+24  1.03E+25  4. OOE+OO-5. OOE+OO  1.12E+23  1.12E+23 
3 . 00E-01 -4 . 50E-01  2.04E+24  1.3QE+25  5. 00E+00-6. OOE+OO  3.S7E+22  3.97E+22 
4 . 50E-01-7 . 00E-01  3.04E+24  1.22E+25  6 . OOE+OO-7 . OOE+OO  1.39E+22  1.39E+22 
7 . 00E-01 - 1 . OOE+OO  2.45E+24  8.17E+24  7 . 00E+00-8 . OOE+OO  5.66E+21  5.66E+21 
1 . OOE+OO- 1 .50E+00  2.11E+24  4.22E+24  8  OOE+OO-1 .00E+01  2.36E+21  1.18E+21 
1.50E+00-2. OOE+OO  9.64E+23  1.93E+24  1 .00E+01-1 .20E+01  0. OOE+OO  0. OOE+OO 

*********************************** 

ATR  PROBLEM  NUMBER  1  SAMPLE  PROBLEM  4  -  ATR  VERSION  5 

*********************************** 

FISSION  PRODUCT  YIELD-2. 0OE+02KT,  FISSION  FRACTION-  9.00E-01 


ATR  PROBLEM  NUMBER 


SAMPLE  PROBLEM  4  -  ATR  VERSION  5 


GROUND  LEVEL 

.500KM, 

59.787GM/CM’  2, 

1.640KFT, 

HORIZ.  RANGE 

RH= 

1.500KM, 

170.018GM/CM”2, 

4.921KFT, 

'SLANT  RANGE 

RS= 

1.530KM, 

175.904GM/CM”2, 

5.0I8KFT, 

TARGET  ALT. 

HT= 

.OOIKM, 

.  117GM/CM”2, 

.003KFT, 

SOURCE  ALT. 

HS= 

.300 KM, 

34.504GM/CM”2, 

.984KFT, 

SLANT  ANGLE 

AN= 

-11.273DEGREES  (COS=  .98071) 

’CALCULATED  FROM  OTHER  COORDINATES 


+  *«>  +  ****  +  +  ***  +  +  **  +  *  +  +  *  +  *  +  *  +  + 


.311MILES 
.932MILES 
. 950MILES 
.OOIMILES 
.  186MILES 


NEUTRON  DOSE  VS.  HORIZ.  RANGE  (RADS  ) 

DOSE  1*  ANS  STANDARD-REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3=  MID-PHANTOM 
DOSE  4-  CONCRETE 
DOSE  5=  AIR 

DOSE  6-  NON-IONIZING  SI. 

DOSE  7*  IONIZING  SILICON 
DOSE  8=  MID-HEAD 
DOSE  9-  1MEV  EQ.  FLUENCE 
D0SE10=  USER  SUPPLIED 

HORIZ.  RANGE  DOSE  1  DOSE  2  DOSE  3  DOSE  4  DOSE  5  DOSE  0 

1 . 500E+03  M  1 . 82E+04  1.52E+03  5.50E+02  1.84E+02  2.05E+O2  1.96E+01 

2 . OOOE+03  M  1 . 18E-I-03  9.89E+01  3.50E+01  1.17E+01  1.70E+01  1.28E+00 

2.500E+03  M  8.78E+01  7.34E+00  2.56E>00  8.S8E-01  1.25E+00  9.51E-02 

3. OOOE+03  M  8.93E+00  5.80E-01  2.01E-01  C.70E-02  9.88E-02  7.53E-03 


DOSE  7  DOSE  8  DOSE  9  DOSEIO 
3.70E+01  7 . 55E+02  5.92E+11  O.OOE+OO 
2.22E+00  4.83E+01  3.86E+10  O.OOE+OO 
1.53E-01  3 . 55E+00  2.88E+09  O.OOE+OO 
1.13E-02  2.79E-01  2.28E+08  O.OOE+OO 


*  *  *  *  * 


GAMMA 


************************** 

DOSE  VS.  HORIZ.  RANGE  (RADS  ) 

DOSE  1=  ANS  STANDARD-REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3-  MID-PHANTOM 
DOSE  4*  CONCRETE 
DOSE  5-  AIR 

DOSE  tt-  IONIZING  SILICON 
DOSE  7-  MID-HEAD 
DOSE  8-  USER  SUPPLIED 


HORIZ.  RANGE 
1 . 500E+03  M 
2 . 000E+03  M 
2 . 500E+03  M 
3 . OOOE+03  M 


DOSE  1  DOSE  2 
7.05E+01  5.84E+01 
7.21E+QO  6.11E+00 
9.42E-01  8.08E-01 
1.42E-01  1.23E-01 


DOSE  3  DOSE  4 
3.80E+01  8.56E+01 
4.04E+00  6.65E+00 
5.41E-01  8.66E-01 
8.32E-02  1.31E-01 


DOSE  5  DOSE  6 
5.29E+01  6.84E+01 
5.53E+00  6.91E+00 
7.31E-01  8.98E-01 
1.11E-01  1.36E-01 


DOSE  7  DOSE  8 
4.29E+01  O.OOE+OO 
4.52E+00  O.OOE+OO 


0.O1E-O1  O.OOE+OO 
9.21E-02  O.OOE+OO 
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*********************************** 


NEUTRON  GAMMA 

DOSE  VS.  HORIZ.  RANGE 

(RADS 

) 

DOSE  1= 

AN'S  STAND.ARD-REM 

DOSE  2= 

soft  tissue 

DOSE  3= 

MID-PHANTOM 

DOSE  4= 

CONCRETE 

DOSE  5= 

AIR 

DOSE  6= 

IONIZING  SILICON 

DOSE  7= 

MID-HEAD 

DOSE  8= 

USER  SUPPLIED 

HORIZ.  RANGE 

DOSE  1 

DOSE  2  DOSE  3 

DOSE  4 

DOSE  5 

DOSE  6 

1 . 500E+03  M 

1.42E+03 

1.26E+03  8.72E+02 

1.32E+03 

1.14E+03 

1.38E+03 

2 . 000E+03  M 

2.01E+02 

1.78E+02  1.24E+02 

1.88E+02 

1.61E+02 

1.94E+02 

2 . 500E+03  M 

3.35E‘01 

2.99E+01  2.07E+01 

3.11E+01 

2.70E+01 

3.24E+01 

3.000E+03  M 

8  24E+00 

5.58E+00  3.89E+00  5.81E+00  5.04E+00 

6.05E+00 

DOSE  7  DOSE  8 

9.67E+02  O.OOE+OO 
1.37E+02  O.OOE+OO 
2.29E+01  O.OOE+OO 
4.27E+00  O.OOE+OO 
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COUPLED  GAMMA 


HORIZ.  RANGE 
1 .500E+03  M 
2.000E+03  M 
2 . 500E+03  M 
3 . 000E+03  M  9 

DOSE  7  DOSE  8 
2 . 17E+03  O.OOE+OO 
2 . 68E+02  O.OOE+OO 
3.90E+01  O.OOE+OO 
6.50E+00  O.OOE+OO 


DOSE  VS 

.  HORIZ.  RANGE 

DOSE  1= 

AN'S  STANDARD-REM 

DOSE  2= 

SOFT  TISSUE 

DOSE  3= 

MID-PHANTOM 

DOSE  4= 

CONCRETE 

DOSE  5= 

AIR 

DOSE  6= 

IONIZING  SILICON 

DOSE  7- 

MID-HEAD 

DOSE  8= 

USER  SUPPLIED 

DOSE  1 

DOSE  2  DOSE 

(RADS 


DOSE  4  DOSE  5  DOSE  6 
3.55E+03  2.93E+03  1.92E+03  3.38E+03  2.66E+03  3.55E+03 

4.20E+02  3.57E+02  2.38E+02  3.92E+02  3.23E+02  4.09E+02 

5.99E+01  5.18E+01  3.51E+01  5.55E<0i  4.69E+01  5.78E+01 

.83E+00  8.80E+00  5.88E+00  9.10E+00  7.78E+00  9.47E+00 
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*********************************** 


TOTAL  DOSE  VS.  HORIZ.  RANGE  (RADS  ) 

DOSE  1=  ANS  ST.AND.ARD-REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3=  MID-PHANTOM 
DOSE  4=  CONCRETE 
DOSE  5=  AIR 

DOSE  6=  IONIZING  SILICON 
DOSE  7=  MID-HEAD 
DOSE  8~  USER  SUPPLIED 

DOSE  1  DOSE  2  DOSE  3  DOSE  4  DOSE  5  DOSE  6 
2. 18E+04  4.45E+03  2.47E+03  3.56E+03  2.92E+03  3.59E+03 

1.60E+03  4 .55E+02  2.73E+02  4.03E+02  3.40E+02  4.11E+02 

1 . 48E+02  5.92E+01  3.76E+0!  5.64E+01  4.81E+01  5.79E+01 

1 .68E+01  9. 18E+00  S.08E+00  9.17E+00  7.83E+00  9.48E+00 

DOSE  7  DOSE  8 
2 . 92E+03  O.OOE+OO 
3. 14E+02  O.OOE+OO 
<*.25E+01  O.OOE+OO 
6.78E+00  O.OOE+OO 


HORIZ.  RANGE 
1.500E+03  M 
2 . 000E+03  M 
2.500E+03  M 
3 . OOOE+G3  M 


FISSION  PRODUCT  DOSE  VS.  HORIZ.  RANGE  (RADS  ) 

DOSE  1=  ANS  STANDARD-REM 
DOSE  2=  SOFT  TISSUE 
DOSE  3=  MID-PHANTOM 
DOSE  4=  CONCRETE 
DOSE  5=  AIR 

DOSE  6=*  IONIZING  SILICON 
DOSE  7=  MID-HEAD 
DOSE  8=  USER  SUPPLIED 


HORIZ.  RANGE 
1 .500E+03  M 
2 . 0O0E+03  M 
2.500E+03  M 
3 . 000E+03  M 


DOSE  1  DOSE  2  DOSE  3  DOSE  4  DOSE  5  DOSE  6 

2.05E+03  1.01E+O3  1.01E+03  1.99E+03  1.40E+O3  2.11E+03 

2.12E+02  1.72E+02  1.10E+02  1.98E+02  l.fE+02  2.08E+02 

2.55E+01  2.12E+01  1.38E+01  2.35E+01  1.92E+01  2.45E+01 

3 . 44E+00  2.90E+00  1.91E+00  3.10E+OO  2.03E+OO  3.28E+00 


DOSE  7  DOSE  8 
1 . 10E+O3  O.OOE+OO 
1  25E+02  O.OOE-tOO 
1 . 55E+01  O.OOE+OO 
2. 14E+00  O.OOE+OO 


^execution  completed 
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APPENDIX  A 
HISTORY  OF  ATR 


Considerable  effort  was  devoted  to  the  under¬ 
standing  of  weapons  radiation  transport  in  the  at¬ 
mosphere,  especially  near  the  air-ground 
interface,  in  the  decade  of  the  60' s.  This  effort 
was  a  combination  of  experimental  and  analytical 
studies,  a  large  portion  of  which  were  supported 
by  DNA  (then  DASA) .  All  the  experimental  work 
was  performed  at  the  Nevada  Test  Site  and  in¬ 
cluded  the  two  tower  experiments,  BREN  (Ref. 
28)  and  HENRE  (Ref.  27) ,  the  former  using  a  bare 
fast  reactor  source  and  the  latter  using  an  accel¬ 
erator-driven  14  MeV  neutron  source.  There  was 
also  interest  in  trying  to  predict  analytically  the 
radiation  dose  and  fluence  levels  measured  in  the 
NTS  atmospheric  shots  of  the  50' s. 

The  DNA-sponsored  efforts  culminated  in  publi¬ 
cations  by  Straker  (Refs.  53  and  56)  and  by  L.G. 
Mooney  and  R.L.  French  (Ref.  10).  The  Mooney 
and  French  work  was  the  most  comprehensive 
ever  attempted  to  completely  describe  the  initial 
radiation  intensity  from  actual  tests,  selected 
from  Teapot,  Plumbbob  and  later  series.  It  used 
the  Straker  air-over-ground  data,  adjusted  for 
height-of-burst  using  the  method  of  French  (Ref. 
24).  and  incorporated  a  method  for  estimating 
the  contribution  of  debris  radiation  to  the  total. 
Mooney  and  French  developed  their  debris  radi¬ 
ation  predictive  model  based  on  the  work  of 
Loewe  (Ref.  32). 

The  work  of  Mooney  and  French  showed  quite  en¬ 
couraging  agreement  with  test  shot  data.  Unfor¬ 
tunately.  no  attempt  was  made  to  incorporate 
their  methods  into  a  general  tool  for  initial  weapon 
radiation  transport  prediction.  However,  a  pro¬ 
posal  to  do  exactly  that  was  submitted  to  DNA 
(DASA  (RARP))  by  Science  Applications,  Inc.  in 
1970.  Oddly,  this  proposal  did  not  arise  from  re¬ 
view  of  the  RRA  work,  but  from  needs  perceived 
by  Dr.  James  Lonergan  of  SAIC,  pertaining  to 
in-flight  radiation  vulnerability  analysis  of  strategic 
missiles.  Dr.  Lonergan 's  concept  was  to  couple 
the  energy  and  angle-dependent  response  calcu¬ 
lated  in  the  adjoint  mode  for  important  internal 
missile  components  with  the  results  of  uniform 
air  transport  calculations  of  the  type  published  by 
Straker.  To  do  this  efficiently,  it  was  necessary 
that  the  air  transport  data  be  placed  In  a  format 
suitable  for  easy  retrieval  by  computer.  SAIC.  un¬ 
der  the  resulting  contract  and  follow-ons  and  un¬ 
der  an  additional  study  funded  by  Ballistic 


Research  Laboratories  (BRL),  produced  five 
successive  versions  of  that  prediction  tool,  which 
was  given  the  name  A«r  Iransport  of  filiation 
(ATR) .  The  documentation  for  these  versions  of 
ATR  are  listed  in  Table  13. 

The  code  released  as  Version  1  was  an  interim 
version  and  included  only  the  neutron  and  sec¬ 
ondary  gamma-ray  portion  of  the  data  base.  Ver¬ 
sion  2  was  actually  the  first  complete  version  of 
ATR  as  originally  proposed.  It  had  functional  fits 
of  the  Straker  energy  and  angular-differential  uni¬ 
form  air  neutron  and  secondary  gamma-ray 
transport  data  and  an  additional  data  base  for 
gamma  ray  and  x-ray  transport,  generated  in  a 
like  format  to  the  Straker  data.  It  also  incorpo¬ 
rated  an  air-over-ground  transport  correction  ap¬ 
proach  for  integral  dose  based  on  the  French 
method  and  normalized  to  the  Straker 
air-over-ground  calculated  data.  There  were  no 
plans  at  the  time  to  incorporate  a  debris  radiation 
model  in  the  code.  The  reason  for  this  is  that  radi¬ 
ation  transport  research  of  the  early  70 's  was 
preoccupied  with  strategic  applications,  particu¬ 
larly  the  penetration  of  missile  silos.  Works  such 
as  that  by  Mynatt,  et  al.  (Ref.  42),  indicated  that 
the  most  important  components  of  initial  weap¬ 
ons  radiation  for  this  application  were  the  neu¬ 
trons,  which  produced  gamma  rays  in  the  silo 
cover,  and  the  high-energy  secondary  gamma 
radiation,  particularly  that  produced  above  the 
silo.  The  moderately  energetic  (-1 .0  MeV)  debris 
radiation  is  much  less  penetrating  through  such 
facilities,  particularly  if  incident  at  an  oblique 
angle. 

In  1972,  new  neutron  cross  section  evaluations 
for  nitrogen  and  oxygen,  sponsored  by  DNA 
(RARP) ,  wero  published  by  Young  and  Foster  of 
LASL  (Refs.  61  and  62).  These  cross  sections 
predicted  half  of  the  secondary  gamma  radiation 
from  14  MeV  neutrons  in  nitrogen  called  for  by 
the  earlier  Straker  set.  The  reason  given  was  that 
the  earlier  set  had  underestimated  the  competi¬ 
tion  from  charged  particle  emission  at  high  neu¬ 
tron  energies.  At  the  same  time,  the 
Young-Foster  data  also  predicted  higher  trans¬ 
mitted  dose  of  high  energy  neutrons  and  greater 
secondary  gamma-ray  production  by  low  energy 
neutrons  that  did  the  Straker  set.  The  appear¬ 
ance  of  the  new  cross  sections  caused  consider- 
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Table  13.  ATR  documentation  summary. 


REPORT 


SAI-71-565-IJ 
November  1971 


DNA2803I 
May  1972 


DNA3144A 
April  1973 

DNA3279T 
August  1974 


DNA3362Z 
August  1974 

DNA3819F 
July  1975 


0NA4061 
January  1976 

DNA3395F 
January  1976 

BRl  CR  343 
August  1977 


CONTENT 


Paper  given  at  RSIC  Workshop  on  Radiation  Transport 
in  Air.  Describes  logic  of  code  and  seme  data  base 
development. 

First  report  describing  basic  concepts  of  ATR  and 
the  first  distributed  version  of  the  code;  describes 
neutron  and  secondary  gamma-ray  data  base  genera¬ 
tion.  Includes  air/ground  and  exponential  air 
correction  factors. 

Users  manual  for  ATR-2  version  of  the  code;  does 
not  describe  data  base  generation. 

Describes  data  base  generation  for  photons  (prompt 
gamma  rays  and  x-rays)  and  the  prompt  gamma-ray 
air/ground  correction  factors. 

Summary  of  the  capabilities  of  the  ATR  code  with 
updates  to  ATR-2. 

Describes  ATR-3  including  new  data  base  using  DMA 
cross  section  library,  new  air/ground  correction 
factors,  low  energy  x-rays,  and  new  REGROUP  routine. 

Describes  TDATR,  the  time-dependent  prompt  photon 
and  secondary  ganma-ray  version  of  ATR. 

Describes  fission  product  model  and  summarizes  total 
capability  of  the  ATR-4  code. 

Describes  ATR-4. 1  with  energy-dependent  air/ground 
correction  factors.  Work  supported  by  BRL. 
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able  argument.  It  was  found  that  they  predicted 
the  results  measured  in  laboratory  experiments 
but  did  not  predict  the  results  of  the  HENRE  ex¬ 
periment.  The  opposite  was  found  to  be  true  of 
the  Straker  set  (Ref.  54). 

Accepting  the  laboratory  experiments  as  the 
more  reliable  comparison.  DNA  and  the  Cress 
Section  Evaluation  Working  Group  (CSEWG)  ac¬ 
cepted  the  Young-Foster  data  to  supplant  that  of 
Straker  in  the  Evaluated  Nuclear  Data  File  (ENDF) 
library.  Subsequently,  it  was  decided  by  DNA  to 
sponsor  an  update  of  ATR  to  include  air  transport 
data  generated  using  the  Young-Foster  cross 
sections.  SAIC  produced  the  homogeneous  air 
transport  results  in  a  format  identical  to  the  origi¬ 
nal  Straker  calculations  and  Oak  Ridge  National 
Laboratory  produced  a  set  of  source  altitude- 
dependent  integral  dose  calculations  (Ref.  44), 
all  under  DNA  sponsorship.  SAIC  then  incorpo¬ 
rated  these  data  into  ATR,  applying  the 
air-over-ground  correction  factors  for  14  MeV 
neutron-produced  integral  dose  to  all  neutron 
source  energies  above  5  MeV  and  those  for  the 
fission  weapon  source  to  all  neutron  source  ener¬ 
gies  below  5  MeV.  Improvements  were  also  made 
in  the  x-ray  transport  at  the  time.  The  results  were 
published  as  Version  3  of  ATR. 

During  the  period  in  which  Version  3  of  ATR  was 
being  produced,  a  new  interest  in  debris  radiation 
developed.  The  driving  force  behind  this  interest 
was  the  possibility  that  low  altitude  ABM  missiles 
might  have  to  fly  through  the  debris  cloud  of 
another  such  missile.  Both  Radiation  Research 
Associates  (RRA)  and  SAIC  were  sponsored  by 
DNA  to  study  the  debris  radiation  problem.  The 
SAIC  product  was  the  Integral  Dose  Environment 
Analysis  (IDEA)  Code  (Ref.  55),  which  used  the 
ATR3  data  base  and  “pr"  scaling  to  estimate  the 
transport  of  radiation  from  all  weapon  source 
components  through  a  perturbed  atmosphere. 
The  time  dependent  hydrodynamic  behavior  of 
the  atmosphere  surrounding  a  burst  was  de¬ 
scribed  using  a  complex  model  (LAMB)  (Ref.  43) 
produced  by  the  Air  Force  Weapons  Laboratory. 
Near  the  end  of  ATR  3  development,  it  was  de¬ 
cided  to  incorporate  a  debris  radiation  model  on 
ATR.  However,  it  was  found  that  the  hydrodynam¬ 
ic  and  cloud  rise  model  used  in  IDEA  was  too  long 
running  and  complex  to  use  as  it  stood.  There¬ 
fore,  it  was  decided  to  use  the  simpler  models 
developed  by  RRA  (Ref.  41).  These  were  essen¬ 
tially  the  same  models  used  to  compute  debris 
radiation  dose  in  the  RRA  report  on  weapon  tests 
referred  to  previously.  Air-over-ground  correc¬ 
tion  factors  for  the  debris  radiation  dose  were 


based  on  the  French  “First-Last  Collision  Meth¬ 
od,  "  as  were  those  for  prompt  gamma  rays.  This 
was  known  not  to  be  very  satisfactory  for  gamma 
ray  applications,  however,  no  suitable  body  of 
data  similar  to  the  Pace  calculations  existed  for 
gamma-ray  transport  in  an  air-over-ground  con¬ 
figuration.  The  resulting  code  became  Version  4 
of  ATR,  the  first  version  of  the  code  capable  of 
calculating  all  of  the  most  important  initial  radi¬ 
ation  components.  The  capabilities  of  ATR4  are 
summarized  in  Figure  21. 

Version  4  of  ATR,  released  in  January  of  1976 
was,  in  fact,  the  conceptual  embodiment  of  all 
aspects  of  the  RRA  study  of  test  shot  dose  distri¬ 
bution  analysis  published  almost  exactly  6  years 
earlier.  Some  nuclear  data  changes  had  oc¬ 
curred  in  the  interim,  but  that  original  concept  of 
estimating  the  radiation  dose  distribution  result¬ 
ing  from  any  device,  based  on  a  single  compre¬ 
hensive  set  of  pre-calculated  data  came  to 
fruition  in  ATR4. 

In  April  of  1973.  a  report  containing  a  large  com¬ 
pilation  of  air-over-ground  transport  calculations 
was  published  by  the  Naval  Weapons  Evaluation 
Facility  (Ref.  5).  It  was  thought  at  the  time  that 
those  results  might  be  used  to  produce  the 
air-over-ground  correction  factors  for  use  in 
ATR3  as  then  contemplated.  However,  upon  re¬ 
view  by  DNA  and  SAIC.  it  was  found  that  the 
NWEF  data  contained  many  errors  and  inconsis¬ 
tencies  which  made  it  unsuitab’e  for  such  use. 
Among  these  were  the  use  of  Pre-Young-Foster 
cross  sections,  insufficient  problem  dimensions 
and  oscillations  in  the  fine-mesh  dose  distribution 
data.  The  results  of  this  review  were  forwarded 
to  NWEF.  In  spite  of  this,  a  group  at  Los  Alamos 
Scientific  Laboratory  (LASL),  including  T.W. 
Dowler  and  H.A.  Sandmeier,  used  the  NWEF  data 
base  to  calculate  dose  distributions  from  actual 
tactical  nuclear  devices.  The  resulting  report  was 
reviewed  by  DNA,  which  determined  that  it  was 
necessary  to  rebut  its  findings.  At  first,  it  was 
thought  that  ATR3,  then  in  preparation  might  be 
used  for  this  purpose.  However,  it  was  finally  de¬ 
cided  by  DNA  that  such  a  rebuttal  must  be  done 
on  the  basis  of  state-of-the-art  transport  technol¬ 
ogy  in  order  to  have  the  greatest  effect.  There¬ 
fore,  6. A.  Straker,  by  then  at  SAIC,  was 
sponsored  to  perform  a  complete  set  of  two-di¬ 
mensional  neutron,  secondary  gamma-ray 
prompt  gamma-  ray  and  debris  gamma-ray  cal¬ 
culations  applicable  to  tactical  device  detona¬ 
tions.  It  was  further  decided  that,  since  the 
number  of  response  functions  of  interest  was  lim- 
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ited,  whereas  the  number  of  sources  was  poten¬ 
tially  limitless,  the  problems  would  be  run  in  the 
adjoint  mode.  This  would  provide  integrated  dose 
data  at  any  range  of  interest  for  an  arbitrary  iso¬ 
tropic  source  at  any  altitude .  As  part  of  this  effort, 
it  was  decided  to  investigate  the  offset  of  soil 
composition,  particularly  moisture  content,  on 
radiation  transport.  A  one-dimensional  analysis 
of  this  effect  was  performed,  revealing  that  such 
a  dependence  did  exist.  Therefore,  pains  were 
taken  to  perform  the  2-D  calculation  using  an 
area-weighted  central  German  soil. 


The  results  of  the  DNA-sponsored  tactical  nuclear 
weapon  radiation  environment  study  (Ref.  26) 
were  published  in  1976.  Shortly  thereafter  it  was 
proposed  by  SAIC  to  DNA  that  the  full  source- 
energy-dependent  air-over-ground  correction 
factors  for  all  radiation  types  available  from  the 
new  calculations  be  incorporated  in  ATR.  The 
proposal  was  not  accepted  by  DNA.  However, 
Ballistic  Research  Laboratories  (BRL)  did  offer  to 
sponsor  the  project.  The  result  of  this  was  ATR4. 1 
(Ref.  34),  which  differs  from  ATR4  as  described 
ir.  Table  14. 


Table  14.  Comparison  of  transport  data  bases  used  for  air- 

over-ground  correction  factors  in  ATR4  and  ATR4.1. 


ATR4 

ATR4.1 

Description 

DOT  calculations  by 

Pace,  ORNL 

TWEEDEE  data  base  from  Adjoint 
DOT  calculations  by  Grltzner, 
Straker 

Source  Energy 

Dependence 

Two  sources 
o  14  MeV 

o  Degraded  fission 

23  neutron  source  energy  bands 
from  15  MeV  to  thermal.  19 
photon  source  energy  bands 
from  10  MeV  to  0.01  MeV 

HOB 

Three  source  heights 

Radiation  Components 

Tissue  dose  for  neutrons 
and  secondary  ganrmas 

Tissue  dose  for  neutrons, 
secondary  gammas,  end  prompt 
gammas 

Energy  and  Angular 
Dependent  Effects 
at  Detector 

No 

No 

Soil 

NTS 

Dry  (no  free  water)  area 
weighted  West  German 

Target  Altitude 

Continuous  from  0  to  1000  m 

1  m 

Ground  Range 

M  .3  km 

■vi.3  km 

APPENDIX  B 


EQUATIONS  TO  CALCULATE  METEOROLOGICAL  VARIABLES 
PRESENT  IN  A  TIST  ATMOSPHERE 

The  following  paragraphs  briefly  describe  a  pro-  Equations  and  constants  used  in  this  procedure 
cedure  to  compute  the  more  important  meteoro-  for  calculating  these  variables  were  obteined 
logical  variables  present  in  a  moist  atmosphere.  from  References  29.  23  and  50. 


Equations  Used 

(1)  Latent  Heats  (L) 

(a)  Vaporization  (T  £  -20  °C)  =  597.2  -  .543t  °C  (cal/g)  (5) 

(b)  Sublimation  (T  <  -20  °C)  =  677.0  -  .0356t  °C  (cal/g)  (6) 

(2)  Dew  Point  Temperature  (Td)  /  Frost  Point  Temperature  (Tp)  for  T  <  -  20  ”C 

Tp  =  C,  /  {C2  -  In  [(RH/100)  x  exp  (C2.  C,/T))}  (K)  (7) 

where:  Ci  =  L  /  Rv  (8) 

C2  =  L  /  273.155  x  Rv  (9) 

and  Ry  is  the  gas  constant  for  water  vapor  =  0.11025  cal/g  K. 

(3)  Relative  Humidity  (RH%) 

RH  =  exp  [C,  x  (t/T  -  1/Td)J  x  100  T.  Td  in  K  (%)  (10) 

(4)  Saturation  Vapor  Pressure  (es) 

e,  *  6.11  x  exp  (C2  -  C)/T)  (mb)  (11) 

(5)  Saturation  Mixing  Ratio  (Ws) 

W,  =  «  x  e,  /  (P  -  e,)  (ppm)  multiply  by  103  to  get  g.'kg.  (12) 

where:  P  is  pressure  (mb) 

and  «  =  (My  /  Md)  =  (18.016  g/mole  /28 . 966  g/mole)  =  .622  (13) 

(6)  Mixing  Ratio  (w) 

w  =  RH  x  W,  /  100  (ppm)  ( i  4) 

also  w  =  pv  /  pd  (15) 


where  pv  =  density  of  water  vapor  present 

Pd  =  density  of  dry  air  containing  the  vapor 

(7)  Vapor  Pressure  (e) 

e  =  w  x  P  /  (t  t  w)  (mb)  (16) 

(8)  Virtual  Temperature  (T*) 

T*  =  T  x  ((1  +  w/«)  /  (1  +  w))  (K)  (17) 
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(9)  Density  of  Moist  Air  (pm) 

Pm  =  P  /  (Rd  x  T*)  =  3.4838  x  10‘4  (P  /  T*)  (g/cm3)  (18) 

where  Rd  is  the  gas  constant  for  dry  air  =  .06857  cal/gm 

(10)  Density  of  Dry  Air  (pd) 

I'd  =  Pm  Ml  +  w)  (g/cm3)  (19) 

(11)  Absolute  Humidity  pv  =  pm  -  pd  (g  cm-3)  (20) 

(12)  Specific  Humidity  (q) 

q  =  pv  /  pm  from  (6)  we  see  that  q  :  w  (21) 

(13)  Potential  Temperature  (0) 

0  =  T  x  (  1000  /  (P  -  e)Kd  (K)  (22) 

where:  T  is  in  K 

and  Kd  =  Rd  /  Cp  (23) 

and  Cp  =  specific  heat  capacity  at  constant  pressure  =  2.4  cal/g  K. 

(14)  Equivalent  Potential  Temperature  (0e) 

0e  =  ©  x  exp  [L  x  W,  /  Cp  x  T)  (K)  (24) 

where  t  is  in  K. 
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APPENDIX  C 


GENERATION  AND  TESTING  OF  ATR5  RESPONSE  FUNCTIONS 


NEUTRON  RESPONSE  FUNCTIONS. 


Neutron  response  functions  span  the  energy 
range  from  1.00  x  10'11  MeV  to  14.9  MeV  using 
22  energy  groups.  Energy-specific  KERMA  values 
were  collapsed  into  this  22  group  format  using  a 
fluence  spectrum  that  consisted  of  a  300  K 
Maxwell-Boltzman  distribution  for  neutron  ener¬ 
gies  below  5  KT  (0.129  eV)  and  a  1/E  distribution 
for  neutron  energies  above  5  KT.  This  ffuence 
spectrum  is  representative  of  the  spectrum  that 
exists  at  distances  of  500  m  io  5000  m  in  air  from 
a  typical  fission  source. 


The  ANS  Standard  response  function  is  based  on 
the  neutron  flux-to-dose  rate  factors  published  by 
the  American  Nuclear  Society  (Ref.  3). 

The  Soft  Tissue  response  function  is  based  on 
Kerr's  model  (Ref.  35). 

The  Mid-Phantom  and  Mid-Head  response  func¬ 
tions  are  based  on  the  work  described  in  Appen¬ 
dix  D. 

The  Concrete  response  function  is  for  the  con¬ 
crete  composition  given  in  Table  15.  The  elemen¬ 
tal  KERMA  factors  were  generated  from  the  DNA 
Working  Cross  Section  Library  (Ref.  1 3)  using  the 
MACK  code  (Ref.  1). 


Table  ’5.  Composition  of  concrete*. 


Element 

Hydrogen 

Oxygen 

Sodium 

Magnesium 

Aluminum 

Silicon 

Sulfur 

Potassium 

Calcium 

Iron 

Total 


wVoinht  Fraction 

9.0056 

0.4S(i3 

0.0171 

0.0024 

0.0456 

0.3158 

0.0012 

0.0192 

0.0826 

0.0122 

1.0000 


•  Thif  cone'*,*  <  t*>m*o  Whn*-O'ort»l*m  tyo'eal  cofnoo*i<mo  w  tt*«*r*ne*  5 1  and  "NB9  eoncr*'*  w  Aatatanca  ‘>0  it  hat  a  penalty  of  !  3$ 

QiCITO 
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The  Dry  Air  response  function  is  for  the  air  com¬ 
position  given  in  Table  16.  The  elemental  KERMA 
factors  were  generated  from  the  DNA  Working 
Cross  Section  Library  using  the  MACK  code 

The  Non-Ionizing  Silicon  and  Ionizing  Silicon  re¬ 
sponse  functions  are  based  on  the  work  reported 
by  Rogers  et  al.  (Refs  45  and  46)  and  their  tabu¬ 
lated  results  (Ref.  59).  The  1  MeV  Equivalent  Flu- 
ence  response  function  is  based  on  the 
recommendations  of  Scott  (Ref.  49). 

GAMMA  RAY  RESPONSE  FUNCTIONS. 

Gamma  ray  response  functions  span  the  energy 
range  from  0.01  MeV  to  12.0  MeV  using  18  ener¬ 
gy  groups.  Energy-specific  KERMA  values  were 
collapsed  into  this  18  group  format  using  a  flu- 
ence  spectrum  that  consisted  of  an  E3  distribu¬ 
tion  for  photon  energies  below  0.045  MeV  and  a 
1/E  distribution  for  photon  energies  above  0  045 
MeV.  This  fluence  spectrum  is  representative  of 
the  spectrum  that  exists  at  distances  of  500  m 
to  5000  m  in  air  from  a  prompt  gamma  fission 
source. 

The  ANS  Standard  response  function  is  based  on 
the  gamma-ray  !lcx-to-dose-rate  factors  pub¬ 
lished  by  the  American  Nuclear  Society  (Ref.  3). 


The  Soft  Tissue  response  function  is  based  on 
Kerr’s  model  (Ref.  35). 

The  Mid-Phantom  and  Mid-Head  response  func¬ 
tions  are  based  on  the  work  described  in  Appen¬ 
dix  C. 

The  Concrete  response  function  is  based  on  Hub- 
bell's  work  (Ref.  31).  The  concrete  composition 
is  given  in  Table  15. 

The  Dry  Air  response  function  is  based  on  Hub- 
bell’s  work  (Ref.  31).  The  air  composition  is  given 
in  Table  16. 

The  Ionizing  Silicon  response  function  is  based 
on  Hubbell's  work  (Ref.  31). 

COMPARISON  OF  OLD  AND  NEW 
RESPONSE  FUNCTIONS. 

Table  17  provides  representative  neutron  and 
gamma-ray  fluences  for  this  comparison  of  re¬ 
sponse  functions.  The  neutron  fluence  is  for  a 
distance  of  1000  m  from  a  1  KT  fission/fusion 
(0. 5/0.5)  neutron  source:  the  fluence  was  calcu¬ 
lated  using  ATR4  with  subsequent  regrouping  to 
the  ATR5  neutron  energy  group  format. The  gam¬ 
ma-ray  fluence  is  for  a  distance  of  1000  m  from 


Table  16.  Composition  of  dry  air'. 


We'ght  Fraction 

Elsmani 

For  Gamma  KERMA 

For  Neutron  KERMA 

Oxygen 

0.232 

0.245 

Nitrogen 

0  755 

0.755 

Argon 

0  C 1 3 

0  000 

Total 

1  0U0 

1  000 

Th#  f*  <  Kir  compotfUon  uaad  for  gamma  ray  KTHMA  if  »>nal  >n#0  hy  ,n  11  at  r#rxv*#ri  m  a  TO  Rloo#  •(•m«n»al 

P#otroo  ¥  F  HMA  factor*  Mf  not  avaflaof#  for  Arg*^  for  fh#  neutron  Kf  I1MA  of  Or  /  air  Argon  fraction  wA*  ASfignrKl  *0  G-yg«n 
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Table  17.  Representative  neutron  and  gamma-ray  fluences 
used  to  illustrate  difference  in  kerma  factors 
between  ATR4  and  ATR5*. 


Neutron  Fluence  Gamma-Ray  Fluence 


Group  No. 

ATR4  Group 

ATR5  Group 

ATR4  Group 

ATR5  Group 

1 

2.94E-13 

2.94E-13 

1.54E-14 

1.29E-14 

2 

3.25E-13 

3.25E-13 

6J6E-14 

6.78E-14 

3 

3.89E-13 

3.89E-13 

4.15E-14 

2.82E- 14 

4 

5.43E-13 

5.43E-13 

1.89E-14 

3.21E-14 

5 

4.58E-13 

4 .58E-13 

9.49E-15 

1.30E-14 

6 

5.88E-13 

5.88E-13 

I.18E-14 

1.25E-14 

7 

8.55E-13 

1.13E-12 

7.24E-15 

7.98E-15 

8 

9.04E-13 

1.12E-12 

4.S4E-15 

8.70E-15 

9 

2.26E-12 

1.78E-12 

3.08E-15 

6.29E-15 

10 

1.38E-12 

1.60E-13 

4.78E-15 

5.01E-15 

11 

7.39E-13 

1.22E-12 

4.13E-15 

4.21E-15 

12 

3.65E-13 

7.39E-13 

5.01E-15 

6.64F-15 

13 

1.79E-13 

3.65E-13 

4.21E-15 

5.14E-15 

14 

3.74E-14 

1.65E-13 

6.34E-15 

4.20E-15 

15 

1.09E-13 

1.60E-13 

5.14E-15 

3.42E-15 

16 

1 .04E-13 

1.04E-13 

6.08E-15 

1.46E-15 

17 

9.83E-14 

9.83E-14 

3.01E-15 

7.39E-16 

18 

1.09E-13 

1 . 10E-13 

7.39E-16 

0.00E-00 

19 

7.24E-14 

7.13E-14 

20 

4 . 31 E- 14 

4. 3 IE- 14 

n  i 

4*  * 

4.79E-14 

4.79E-14 

22 

1.21E-13 

1.16E-13 

*  The  neutron  fluence  was  calculated  using  ATR4,  and  then  regrouped  into 
ATIto  format.  The  gamma  ray  fluence  was  calculated  using  ATR5,  and  then 
regrouped  into  ATR4  format. 
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£  1  kT  prompt  gamrra  fission  source:  this  fluence  the  fluences  from  Table  17  and  the  old  (ATR4) 

was  calculated  using  ATR5  with  subsequent  re-  and  the  new  (ATR5)  response  'unctions, 

grouping  to  the  ATR4  gamma  ray  energy  group 
format. 

Table  19  shews  me  change  ;n  calculated  re- 

Table  18  shows  the  responses  calculated  using  sponses  for  comparable  response  'unchons. 


Table  18.  Neutron  and  gamma-ray  kerma  and  dose  values 
calculated  from  ATR4  and  ATR5  with  identical 
fluence  values*. 


Neutron  Responses 

ATR4  ATR5 


Henderson  Tissue 

=  8.28E-21 

Snyder-Neufsld 

=  1.51E-20 

Tissue 

=  9.14E-21 

Soft  Tissue 

=  9.32E-21 

Mid-Phantom 

=  3.65E-21 

Mid-Phantom 

=  3.33E-21 

Concrete 

-  1.22E-21 

Concrete 

=  1.24E-21 

Air 

-  1.64E-21 

Air 

-  1.58E-21 

Non-Ionizing  Silicon 

=  1.13E-22 

Non-Ionizing  Silicon 

=  1.18E-22 

Ionizing  Silicon 

-  3.00E-22 

Ionizing  Silicon 

~  3.10E-22 

ANS  Standard 

-  1.09E-19 

Mid-Head 

-  4.52E-21 

1  MeV  Equil .  Flu. 

=*  3.56E-12 

Gamma  Ray 

Responses 

ATR4 

ATR5 

Henderson  Tissue 

=  7.18E-23 

Concrete 

-  8.30E-2 3 

Concrete 

=  7.53E-23 

Air 

-  6.85E-23 

Air 

=  6.30E-23 

Si  1  icon 

-  7.93E-23 

Ionizing  Silicon 

=  7.84E-23 

.ANS  Standard 

=  8.10E-23 

Soft  Tissue 

=  8.96C-23 

Mid-Phantom 

4.69E-23 

Mid-head 

=  5.22E-23 

*  Since  the  ATR4  and  ATR5  fluences  are  based  on  one  fluence  regrouped  into 
the  other  energy  group  format,  the  fluences  are  essentially  equivalent. 
Accordingly,  this  table  shows  the  effect  of  changing  the  response 
function. 
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Table  19.  Changes  in  selected 
f luence-to-kerma  or 

Response 

Tissue 

Mid-Phantom 

Concrete 

Air 

Non-Ionizing  Silicon 
Ionizing  Si  1  icon 


kerma  and  dose  values  due  to 
-dose  factor  changes*. 

Percent  Change 

Neutron 

Gamma  Ray 

+2% 

NA 

-9% 

NA 

+2% 

-0K 

-4* 

-8% 

+4S 

NA 

*3% 

-1* 

*  ("ATR5"  -  "ATR4")  /  "ATR4". 

So,  for  example,  if  ATR4  and  ATR5  were  used  to  calculate  Concrete  kerma,  at 
a  range  of  1000  ■  from  a  1  kT  prompt  gamma  source,  the  ATR5  kerma  would  be 
9^  less  than  the  ATR4  kerma,  assuming  the  transported  fluences  were  the 
same . 
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APPENDIX  D 


KERMA  IN  AN  ANTHROPOMORPHIC  PHANTOM 


This  Appendix  describes  the  procedure  used  to 
calculate  neutron  and  gamma-ray  Mid-Phantom 
and  Mid-Head  energy-angle  differential  response 
functions.  This  task  was  performed  under  a  DNA 
sponsored  effort,  Contract  No.  DNAOQ1-83-C- 
0252. 

The  transport  of  externally  incident  fluence  to 
points  at  the  mid-head  and  mid-thorax  of  an  an¬ 
thropomorphic  phantom  has  been  assessed  us¬ 
ing  the  VCS  code  system  (Ref.  47)  which 
incorporates  the  MORSE  Monte  Carlo  transport 
code  (Ref.  57) .  The  MORSE  code  performs  trans¬ 
port  calculations  in  three  dimensions. 

Computations  were  performed  using  a  three- 
dimensional  adult  male  anthropomorphic  phan¬ 
tom  developed  by  Cristy,  (Ref.  7)  having  a  total 
mass  of  approximately  74  Kg.  This  phantom  and 
its  internal  lung  and  skeletal  detail  are  depicted 
in  Figure  22.  Specifications  of  elemental  constitu¬ 
ents,  a3  recently  published  by  Kerr,  (Ref.  35)  are 
provided  in  Table  20  along  with  revised  density 
values  (Ref.  8). 

Transport  computations  have  been  performed 
using  the  DNA  Few  Group  Cross  Section  Library 
(DLC-31)  with  P3  Legendre  expansion.  Calcula¬ 
tions  have  been  performed  in  the  adjoint  mode. 
In  this  mode  the  physics  of  the  radiation  transport 
are  reversed  by  inverting  the  energy  dependent 
quantities  and  reflecting  the  angular  variables 
through  the  origin.  This  causes  neutrons  and 
gamma  rays  in  the  computation  to  gain  energy 
in  scattering  and  causes  gamma  rays  to  "pro¬ 
duce”  neutrons. 

The  adjoint  approach  is  employed  by  starting  a 
particle  at  a  detector  point  and  following  it  until 
it  leaks  from  the  system  along  a  direction  ft'  and 
at  a  leakage  surface  defined  by  its  normal  vector 
n.  The  correlated  descriptions  of  the  starting  and 
leaking  particles  are  saved. 


The  system  leakage  fluence  (n  •  n’)4>*(E’,  O’) 
may  be  viewed  as  having  the  qualities  of  detector 
weighted  by  the  transport  through  the  intervening 
media.  It  may  be  convoluted  with  any  fluence 
<j>(E’,  D')  defined  as  incident  on  the  system  from 
which  the  leakage  is  scored  to  obtain  the  detector 
response  for  that  particular  incident  fluence.  as 
in  the  expression: 

Total  Kerma  =  /  J  dfi’dE>(E\fl')(n-  n')<j>*(E’, 
n1)  n’  E* 


(25) 

at  some  range  r  from  the  radiation  source.  This 
process  is  shown  schematically  in  Figure  23. 

To  obtain  an  energy  differential  response  func¬ 
tion,  the  angular  integration  is  carried  out  in  ad¬ 
vance,  yielding  the  energy  dependent  (isotropic) 
transport-weighted  Kerma  factors,  <j>‘(E). 

Computations  of  transport  in  the  anthropomor¬ 
phic  phantom  have  been  performed  for  two  de¬ 
tector  locations.  The  first  location  is  mid-head, 
which  is  specified  to  be  the  geometric  center  of 
the  ellipsoid  head  of  the  phantom.  The  second 
location  is  mid-thorax,  which  is  specified  to  be 
at  the  geometric  center  of  the  elliptical  torso 
cross  section,  midway  between  the  top  of  the 
pelvis  and  the  top  of  the  torso.  Starting  histories 
were  sampled  uniformly  from  58  (37  neutron,  21 
gamma-ray)  energy  groups.  A  total  of.  60,000 
starting  histories  were  used  for  the  mid-head  lo¬ 
cation,  whiie  a  total  of  45,300  were  used  for  the 
mid-thorax  location.  Due  to  splittings  of  histories 
in  the  Monte  Carlo  process,  approximately  twice 
the  number  of  starting  histories  were  actually  tal¬ 
lied  as  having  leaked  from  the  system,  a  few 
more  for  the  mid-thorax  location,  a  few  less  for 
the  mid-head  location.  The  results  of  these  calcu¬ 
lations  may  be  considered  to  be  highly  reliable, 
statistically,  in  terms  ol  total  Kerma  transmission. 


87 


ELEVATION  30.0 


flX  I  M'JTH  60.0 

SCALING  20.00 


Figure  22.  Adult  male  reference  man. 


Table  20.  Elemental  composition  for  various  components 
of  reference  man. 


Element  Lung 

H  10.21 

C  10.24 

N  2.91 

0  75.63 

Na  0.19 

Mg  0.007 

P  0.080 

S  0.23 

Cl  0.27 

K  0.20 

Ca  0.009 

Fe  0.037 

Density  (g/cm3)  0.296 


Percent  by  Weight 


Skeleton 

Soft  Tissue 

7.28 

10.51 

24.64 

22.63 

3.06 

2.34 

46.88 

63.69 

0.32 

0.011 

0.11 

0.013 

5.03 

0.13 

,  0.31 

0.20 

0.14 

0.14 

0.15 

0.20 

12.07 

0.024 

0.008 

0.006 

1.40 

1.04 

89 


I®. 


I* 


i  | 

■»  i 


.f 


APPENDIX  E 

DELAYED  GAMMA  RADIATION  TRANSPORT  METHODOLOGY  DEVELOPMENT 


INTRODUCTION. 

The  improvement  of  the  ATR  delayed  radiation 
model  has  been  based  on  the  study  of  empirical 
data  from  the  atmospheric  tests  conducted  in  Ne¬ 
vada  and  Pacific  test  sites  and  the  study  of  alter¬ 
native  models.  This  Appendix  describes  the 
development  and  study  of  an  alternative  model 
to  that  contained  in  ATR.  Study  o'  this  alternative 
model  yielded  significant  insights  into  possible 
shortcomings  of  the  ATR  model,  particularly  the 
hydrodynamic  model.  Unfortunately,  the  alterna¬ 
tive  hydrodynamic  model  was  too  complex  for  di¬ 
rect  inclusion  in  ATR.  However,  it  did  provide 
additional  data  from  which  to  develop  the  scaled 
height-of-bur3t  correction  factor  described  in  the 
body  of  this  report.  In  addition,  it  Drovided  the  fis¬ 
sion  product  radiation  source  terms  incorporated 
in  ATR5. 

Delayed  gamma  radiation  is  that  component  of 
initial  nuclear  weapon  radiation  (radiation  emitted 
within  the  first  60  seconds  fo'lowing  initiation) 
emitted  by  products  of  fission  and  neutron  activa¬ 
tion  which  are  parts  of  the  debris  of  the  nuclear 
explosion.  These  debris  initially  reside  within  the 
fireball,  a  luminous  sphere  of  hot  gases  gener¬ 
ated  by  the  interaction  with  the  air  of  x-rays  ra¬ 
diated  from  the  hot  debris  within  a  millionth  of  a 
second  after  initiation. 

The  fireball  also  defines  the  edge  of  a  region  of 
very  low  density  air,  a  density  well,  created  by 
the  expansion  of  the  very  hot  gases  within.  The 
high  pressure  region  created  at  the  edge  of  that 
density  well  forms  a  shock  wave  which  breaks 
away  from  the  fireball  and  travels  outward,  super¬ 
imposed  on  the  ambient  air.  The  net  effect  of 
these  perturbations  is  to  decrease  the  amount  of 
air  on  a  line  from  source  to  target,  thereby  en¬ 
hancing  the  transmission  of  delayed  radiation 
over  that  which  would  occur  in  unperturbed  air. 
This  effect  is  often  referred  to  as  hydrodynamic 
enhancement. 

The  density  well  Is  buoyant  and  begins  to  rise 
within  a  fraction  of  a  second  of  its  creation.  Its 
momentum  is  augmented  by  the  returning  shock 
wave,  now  reflected  off  the  ground.  The  rise  of 
the  fireball,  which  includes  the  radiation  source. 


causes  the  distance  between  source  and  target 
to  increase,  thereby  causing  the  radiation  trans¬ 
mission  to  decrease  once  the  expansion  of  the 
fireball  has  slowed  or  stopped. 

As  the  fireball  rises  the  outer  surface  cools  and 
is  subject  to  drag  forces  from  the  surrounding  at¬ 
mosphere.  That  causes  the  spherical  fireball  to 
become  a  toroid-shaped  cloud  of  condensed 
gases,  containing  the  weapon  debris. 

Throughout  this  process  the  weapon  debris,  in¬ 
cluding  the  fission  products,  are  in  motion  relative 
to  the  fireball.  First,  the  debris  are  propelled  to¬ 
ward  the  top  of  the  fireball  with  an  initial  impulse 
provided  by  the  reflected  shock.  Eventually  the 
debris  pierce  the  top  of  the  rising  fireball  and  flow 
down  the  side,  with  most  becoming  incorporated 
in  the  torus  itself  or  in  the  cloud  skirt  below  the 
torus.  Those  effects  become  less  pronounced  as 
the  scaled  burst  height  increases,  where  the 
scale  burst  height  is  defined  to  be  the  burst  height 
divided  by  yield  to  the  one  third  power  (HOB/ 

y1/3). 

This  complex  system,  consisting  of  a  radioactive, 
hence  time-varying,  source  of  gamma  rays  and 
similarly  time-varying  geometry,  in  which  the 
amount  of  air  between  source  and  target  is 
changing,  even  as  the  source  rises  from  its  origi¬ 
nal  location,  must  be  modeled  in  some  detail  in 
order  to  properly  describe  the  intensity  of  the 
delayed  radiation  component.  The  motivation  for 
doing  so  is  that  delayed  gamma  radiation  com¬ 
prises  the  majority  of  the  gamma  ray  field 
strength  from  large  yield  weapons. 

Three  technical  areas  must  be  addressed  in  the 
modeling  process.  These  are: 

Fission  Product  Gamma  Ray  Source 

Gamma  Ray  Transport 

Nuclear  Weapon  Hydrodynamics 

Data  required  for  modeling  the  processes  of  in¬ 
terest  in  these  areas  have  been  acquired  through 
laboratory  and  field  measurements.  However, 
these  data  are  by  no  means  complete.  Here  fol¬ 
lows  a  discussion  of  such  modeling. 
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FISSION  PRODUCT  GAMMA  RAY  SOURCE. 

Measurements  of  gamma  and  beta  radiation 
source  rates  as  functions  of  time  aftei  fission 
have  been  made  in  the  laboratory.  Maienscnein, 
et  al.,  (Refs.  38  and  39),  reported  measure¬ 
ments  made  in  the  late  1950's  and  early  1960's 
at  Oak  Ridge  National  Laboratory  in  terms  of 
spectra,  total  number  and  total  energy  emission 
rates.  In  the  early  1960's,  Fisher  and  Engle  (Refs. 
21  and  22),  of  Los  Alamos  National  Laboratory 
performed  similar  experiments  using  a  fast  reac¬ 
tor  source,  which  allowed  them  to  obtain  results 
at  very  early  times  (.35  sec).  They  also  per¬ 
formed  such  measurements  for  several  fission¬ 
able  materials  including  U-235,  Pu-239  and 
U-238. 

Dickens  et  al.,  (Ref.  14  and  15),  have  published 
results  of  recent  measurements  made  ORNL,  in¬ 
cluding  those  of  spectra,  again  using  thermal 
neutron  fission.  Also,  Japanese  scientists  Akiya- 
ma  and  An  (Ref.  2)  have  published  results  of  their 
measurements  using  a  fast  reactor  source.  Th9 
Japanese  data  are  in  terms  of  total  energy  rate 
only,  and  are  limited  to  times  of  twenty  seconds 
or  greater. 

Results  of  several  measurements  in  terms  of  en¬ 
ergy  emission  rates  at  forty  seconds,  the  earliest 
time  common  to  ail  the  available  measurements, 
are  provided  in  Table  21 .  There  is  a  limit  to  what 
can  be  construed  from  these  data  because  of 
disparities  in  the  uncertainties  quoted  by  the  indi¬ 
vidual  measurers.  The  Maienschein  et  al.,  data 
are  assigned  a  standard  deviation  of  approxi¬ 
mately  twenty-five  percent,  those  of  Fisher  and 
Engle  twelve  percent  (U-235  and  U-238)  and 
twenty  one  percent  (Pu-239)  end  those  of  Dick¬ 
ens  et  al, ,  and  Akiyama  and  An  less  than  five  per¬ 
cent.  Nevertheless,  the  data  having  the  lowest 
uncertainties  are  consistent  with  each  other  and 
are  not  inconsistent  with  the  value  provided  by 
Maienschein  et  al.,  whereas  the  data  of  Fisher 
stand  apart  from  the  high  precision  data.  Unfortu¬ 
nately.  the  latter  provide  the  only  available  spec¬ 
tral  data  for  very  early  times  and  for  U-238 
fission.  Therefore,  it  i3  desirable  to  reconcile 
them  to  those  of  (he  other  measurers.  To  accom¬ 
plish  this,  the  spectra  of  Dickens,  et  al.,  and  of 
Fisher  and  Engle  are  compared  for  U-235  and 
Pu-239.  common  to  both  experiments  (Table  22, 
columns  three  and  four).  It  Is  found  that  the  two 
sets  of  spectra  differ  from  each  other  by  the 


same  energy-dependent  ratio  (last  column)  af*er 
removal  of  a  constant  correction  term  in  the  Fish¬ 
er  and  Engle  Pu-239  measurerr  mt  (.65,  column 
5).  That  correction  is  consistent  with  the  Fisher 
and  Engle  disclaimer  concerning  their  estimation 
of  the  number  of  fissions  ;n  the  Pu  sample.  The 
energy  dependent  relationship  found  to  exist  be¬ 
tween  the  Dickens,  ct  al.,  and  the  Fisher  and  En¬ 
gle  U-235  and  Pu-239  may  also  be  applied  to  the 
latter’s  U-238  measurements  The  result,  in 
teims  of  total  gamma  ray  energy  emission  rate 
agrees  within  a  few  percent  with  the  value  re¬ 
ported  by  Akiyama  and  An  at  40  seconds  post¬ 
irradiation  time.  Therefore,  it  is  concluded  that 
a  systematic  error  must  be  present  in  the  Fisher 
and  Engle  data  but  that  it  can  be  removed 
through  the  use  of  the  Dickens,  et  al  -  related, 
energy-dependent  corrections.  The  resulting 
source  rate  models  for  three  representative  en¬ 
ergy  groups  along  with  supporting  experimental 
data  are  shown  as  functions  of  time  in  Figures  24, 
25  and  26.  Figure  26  also  indicates  some  addi¬ 
tional  empirical  corrections  represented  by  the 
dotted  lines  made  to  the  Fisher  and  Engle  data 
to  make  their  temporal  behavior  self-consistent 
and  also  more  consistent  with  those  of  Dickens 
et  al.  The  resulting  source  rate  models  are  not 
only  consistent  with  recent  experiments  but  are 
also  inclusive  of  early  times  and  data  for  U-238. 
Source  rate  spectra  tor  five  times.  .35  through 
40  seconds,  are  given  for  U-235.  U-238  and 
Pu-239  in  Tables  23,  24  and  25,  respectively. 
Source  rates  for  times  other  than  those  shown 
are  obtained  by  1  n-1  r«  interpolation.  For  times 
earlier  than  .35  seconds  the  source  rate  is  ob¬ 
tained  by  In- In  extrapolation,  based  on  the  earli¬ 
est  two  times,  .35  and  1 .5  The  contribution  from 
sources  beyond  30  seconds  is  negligible  There¬ 
fore  ,  no  extrapolation  to  times  greater  than  40  se¬ 
conds  is  performed. 

An  addition  complication  has  not  been  taken  into 
account  in  the  delayed  gamma  radiation  source 
model.  Calculations  using  nuclide-specific  gam¬ 
ma  omission  rates  (Ref.  18)  indicnte  that  fast  fis¬ 
sion  produces  more  gamma  ray  energy  at  early 
times  than  does  thermal  fission  (Table  26).  This 
is  supported  to  some  extent  by  a  comparison  of 
the  new  American  (thermal  fission)  and  Japa¬ 
nese  (fast  fission)  data  in  Table  21 .  Having  in  ef¬ 
fect  normalized  Fisher  and  Engle’s  results  to 
those  of  Dickens,  et  al.,  the  former  are  now  con¬ 
sistent  with  those  for  thermal  rather  than  fast  fis¬ 
sion. 
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Table  21. 


Fission  product  gamma  ray  source  energy  emission 
rate  (0.1  to  5  MeV)  at  40  sec  after  fission. 


Nucl ide 

MeV/fis-sec 

Experimentor 

U-235 

U-2  38 

Pu-229 

Maienschein,  et  al. 

0.0193* 

- 

- 

Fisher  and  Engle 

0.0216 

0.0275 

0.0194 

Dickens,  et  al . 

0.0166 

- 

0.0130 

Akiyama  and  An 

0.0171 

0.0205 

0.0131 

♦Energy  contribution  0.1  <  E  <  0.28  estimated  by  extrapolation. 


Table  22.  Mean  values  of  Dickens  to  Fisher  and  Engle  source 
ratios  by  energy  group  for  U-235  and  Pu-239. 


Group 

No. 

Upper 

Energy 

(MeV) 

V 

U-235 

s 

^FSE 

PU-239 

R239/R235 

VSF&£ 

a 

Pu-239 

Adjusted 

VSFXE 

Mean 

1 

6.419 

•394( . 179 ) b 

. 360( .161) 

.914 

.424 

.409( .170) 

2 

5.636 

.S20( .096) 

. 654( .217) 

.798 

.769 

•795( .164) 

3 

4.918 

.  627 ( .416) 

.630( .103) 

1.005 

.741 

.634( .281) 

4 

4.257 

.866( .095) 

.734 ( .059) 

.848 

.864 

.865( .079) 

5 

3.655 

.958( .048) 

.317 ( .044) 

.853 

.961 

.960( .046) 

6 

3.110 

,852( .076) 

.730( .031) 

.857 

.859 

.856( .058) 

7 

2.620 

.904 ( .046) 

.779( .025) 

.862 

.916 

•910( .037) 

8 

2.189 

•869( .028) 

,738( .026) 

.849 

.868 

•869( .027) 

9 

1.808 

.779( .034) 

.674(.045) 

.865 

.793 

.786( .040) 

10 

1.478 

.787( .049) 

.688( .033) 

.874 

.809 

.798( .042) 

11 

1.195 

.722( .044) 

. 615( .035) 

.852 

.724 

. 723( .040) 

12 

.954 

. 7  51 ( .035) 

. 620( .012) 

.826 

.729 

,740( ,026) 

13 

.749 

•674( .100) 

.515( .084) 

.764 

.606 

•640( .093) 

14 

.575 

. 671 ( .086) 

•569( .082) 

.848 

.669 

. 6  70( .084) 

15 

.428 

. 57  7 ( .106) 

.495( .104) 

.858 

.582 

•580( .105) 

16 

.309 

.484( .038) 

.448(.C94) 

.926 

.527 

.506( .074) 

17 

.213 

.394( .077) 

. 345( .160) 

.876 

.406 

•4U0( .127) 

.137 

mean' 

.854(0.040) 

aRat1o 

of  the  Pu 

-239  source  ratio  to  that  of 

U-235. 

bread 

as  mean  and  (fractional 

standard  deviation  of  the  sample  - 

o  .  /mean), 
n  —  i 

cmean 

of  group  values  E  <  4.257  MeV . 
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Figure  24. 


Fission  product  photon  eaission  rates  vs  tiae, 
<  E  <  2.19  KeV,  for  U-235,  U-239  and  Pu-238. 
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Table  23.  Fission  product  gamma  ray  source  rate  spectra  for 
U-235  at  mean  times  after  fission  of  0.35,  1.5, 
4.75,  11.5  and  40  seconds. 


Upper  Gama  Rays/Fission-second-MeV 

Group  Energy 


No. 

M. 

0.35 

1.5 

4.75 

11.5 

40. 

1 

14.0 

0.0 

0.0 

0.0 

0..0 

0.0 

2 

10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

** 

•J 

o 

• 

CO 

5.50-5* 

2.98-5 

1.51-5 

7.01-6 

2.08-6 

4 

7.0 

3.87-4 

1.80-4 

8.22-5 

2.70-5 

6.30-6 

5 

6,0 

1.10-3 

7.70-4 

3.70-4 

1.79-4 

5.40-5 

6 

5.0 

5.49-3 

2.53-3 

1.12-3 

5.47-4 

1.83-4 

7 

4.0 

1.51-2 

8.89-3 

4.32-3 

1.84-3 

5.60-4 

8 

3.0 

2.63-2 

1.49-2 

6.59-3 

3.32-3 

1.23-3 

9 

2.5 

4.19-2 

2.34-2 

1.17-2 

5.34-3 

1.80-3 

10 

2.0 

6.31-2 

3.65-2 

1.87-2 

8.34-3 

2.78-3 

11 

1.5 

1.27-1 

7.21-2 

3.30-2 

1.63-2 

5.43-3 

12 

1.0 

2.12-1 

1.12-1 

4.88-2 

2.30-2 

7.39-3 

13 

,;o 

2.66-1 

1.61-1 

8.90-2 

3.91-2 

9.20-3 

14 

.45 

2.73-1 

1.43-1 

7.91-2 

3.70-2 

1.09-2 

15 

.30 

4.08-1 

1.85-1 

1.08-1 

5.06-2 

1.35-2 

16 

.15 

1.94-1 

7.82-2 

4.66-2 

2.11-2 

5.61-3 

17 

.10 

0.0 

0.0 

0.0 

0.0 

0.0 

18 

.07 

0.0 

0.0 

0.0 

0.0 

0.0 

19 

.045 

0.0 

0.0 

0.0 

0.0 

0.0 

20 

.030 

0.0 

0.0 

0.0 

0.0 

0.0 

21 

.020 

0.0 

o 

• 

o 

0.0 

0.0 

0.0 

.010 


•read  as  5.5*10"S. 
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Table  24.  Fission  product  gamma  ray  source  rate  spectra  for 
U-238  at  mean  times  after  fission  of  0.35,  1.5, 
4.75,  11.5  and  40  seconds. 


Upper  Gamma  Rays/F1ssion-$econd-MeV 

Group  Energy 


No. 

(MeV) 

0.35 

1.5 

4.75 

11.5 

40. 

1 

14.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

8.0 

1.26-4* 

5.87-5 

2.58-5 

9.84-6 

2.78-6 

4 

7.0 

1.27-3 

3.71-4 

1.31-4 

4.33-5 

7.01-6 

5 

6.0 

3.50-3 

1.43-3 

6.60-4 

2.78-4 

6.16-5 

6 

5.0 

1.25-2 

4.95-3 

1.77-3 

7.91-4 

1.83-4 

7 

4.0 

3.58-2 

1.82-2 

7.38-3 

2.76-3 

6.72-4 

8 

3.0 

7.33-2 

3.38-2 

1.15-2 

4.76-3 

1.68-3 

9 

2.5 

1.13-1 

4.83-2 

1.98-2 

7.89-3 

2.40-3 

10 

2.0 

1.41-1 

7.77-2 

3.10-2 

1.17-2 

3.18-3 

11 

1.5 

2.74-1 

1.39-1 

5.20-2 

2.08-2 

6.56-3 

12 

1.0 

4.40-1 

2.07-1 

8.00-2 

3.05-2 

9.69-3 

13 

.70 

5.61-1 

2.80-1 

1.44-1 

5.30-2 

1.35-2 

14 

.45 

6.70-1 

3.00-1 

1.40-1 

5.90-2 

1.54-2 

15 

.30 

9.9S-1 

4.01-1 

2.07-1 

7.38-2 

2.00-2 

16 

.15 

5.02-1 

1.57-1 

9.70-2 

3.17-2 

7.69-3 

17 

.10 

0.0 

0.0 

0.0 

0.0 

0.0 

18 

.070 

0.0 

0.0 

0.0 

0.0 

0.0 

19 

.045 

0.0 

0.0 

0.0 

0.0 

0.0 

20 

.030 

0.0 

0.0 

0.0 

0.0 

0.0 

21 

.020 

0.0 

0.0 

0.0 

0.0 

0.0 

.010 


•read  as  i.26*10"4. 
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Table  25.  Fission  product  gamma  ray  source  rate  spectra  for 
Pu-239  at  mean  times  after  fission  of  0.35,  1.5, 
4.75,  11.5  and  40  seconds. 


Group 

No. 

Upper 

Energy 

(MeV) 

0.35 

Gamma  Rays/Fission- 

1.5  4.75 

-second-MeV 

11.5 

40. 

1 

14.0 

0.0 

0.0 

o.c 

0.0 

0.0 

2 

10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

8.0 

4.56-5* 

2.63-5 

1.25-5 

5.65-6 

1.64-6 

4 

7.0 

1.93-4 

1.07-4 

4.85-5 

1.61-5 

3.97-6 

5 

6.0 

7.12-4 

3.39-4 

1.69-4 

7.48-5 

2.36-5 

6 

5.0 

4.08-3 

1.75-3 

5.76-4 

2.77-4 

8.47-5 

7 

4.0 

1.06-2 

5.68-3 

2.79-3 

1.16-3 

3.91-4 

8 

3.0 

2.38-2 

1.07-2 

4.73-3 

2.40-3 

1.02-3 

9 

2.5 

3.31-2 

1.74-2 

8.55-3 

4.03-3 

1.52-3 

10 

2.0 

4.86-2 

2.81-2 

1.52-2 

6.79-3 

2.15-3 

11 

1.5 

9.73-2 

6.00-2 

2.50-2 

1.20-2 

4.01-3 

12 

1.0 

1.70-1 

9.17-2 

3.95-2 

1.79-2 

5.68-3 

13 

.70 

2.43-1 

1.68-1 

8.06-2 

3.27-2 

7.58-3 

14 

.45 

2.34-1 

1.42-1 

6.80-2 

3.28-2 

9.31-3 

15 

.30 

3.57-1 

1.82-1 

9.91-2 

4.63-2 

1.08-2 

16 

.15 

1.77-1 

7.45-2 

4.27-2 

1.77-2 

4.40-3 

17 

.10 

0.0 

0.0 

0.0 

0.0 

0.0 

18 

.070 

0.0 

0.0 

0.0 

0.0 

0.0 

19 

.045 

0.0 

0.0 

0.0 

0,0 

0.0 

20 

.030 

0.0 

0.0 

0.0 

0.0 

0.0 

21 

.020 

.010 

0.0 

0.0 

0.0 

0.0 

0.0 

-5 


♦read  as  4.56*10 


Table  26. 


Fast/thermal  fission  product  gamma  ray  energy 
ratios . 


Calculated  Using  ENDF/B-V 


Nuclide 


Time  (sec) 

U235 

Pu239 

0.1 

1.049 

1.072 

1.0 

1.062 

1.062 

5.0 

1.025 

1.037 

10.0 

1.004 

1.029 

20.0 

0.993 

1.020 

Calculated  by  England  and  Schenter 

0.1 

1.038 

0.945 

1.0 

1.067 

0.986 

5.0 

1.052 

1.002 

10.0 

1.032 

0.997 

20.0 

1.020 

1.002 

50.0 

1.018 

1.019 

emission 
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No  action  was  taken  to  correct  tne  source  for  the 
fast-thermal  fission  effect  at  the  time  the  correc¬ 
tions  were  made  to  the  Fisher  ana  Engle  source 
data.  This  was  due  to  questions  regarding  the  ac¬ 
curacy  of  calculations  using  individual  nuclide 
data.  Those  questions  were  apparently  well 
founded,  at  least  for  Pu239,  based  on  the  more 
recent  analysis  of  England  and  Schenter  (Ref. 
19) ,  data  from  which  also  appear  in  Table  26.  T^e 
ATR5  delayed  gamma  ray  energy  emission 
source  rates  are  taken  from  England  and  Schent¬ 
er  data  applicable  to  fast  neutron  fission. 

GAMMA  RAY  TRANSPORT. 

The  capability  of  current  codes  and  cross  sec¬ 
tions  to  model  gamma  ray  transport  in  uniform 
air  was  tested  using  BREN  Co-60  source- 
measurement  data  (Ref.  28)  as  shown  in  Figure 
27.  Those  calculations  were  performed  for  a 
source  343m  above  ground  using  the  two- 
dimensional  discrete  ordinate  code  DOT  (Ref. 
25)  with  the  36  energy  group  Vitamin  C  cross 
section  set  (Ref.  43).  Ideally,  it  would  be  desir¬ 
able  to  use  the  DOT  code  to  calculate  the  delayed 
radiation  transport.  However,  because  the  de¬ 
bris/fireball/  ground  geometry  is  dynamic,  it 
would  require  thirty  or  more  separate  calculations 
to  adequately  represent  the  propagation  of 
delayed  radiation  from  any  given  weapon  burst. 
Therefore,  an  alternative  method  has  been  found 
which  provides  a  more  practical  approach  to  the 
problem.  According  to  this  approach  the  gamma 
ray  transport  is  modeled  in  a  one-dimensional 
system  of  uniform  air.  The  amount  of  air  along 
a  line  from  the  debris  (the  gamma  ray  source) 
to  the  target  at  a  specific  time  is  determined  us¬ 
ing  a  hydrodynamic  model,  which  will  be  dis¬ 
cussed  later  in  this  section.  The  energy  and 
angle-differential  fluence  rate  at  the  target  is  de¬ 
termined  for  the  time-specific  source  rate  spec¬ 
trum  transported  through  an  identical  thickness 
of  air,  using  the  uniform  air  transoort  data  base. 
This  umform-for  perturbed  air  substitution  is  a 
good  approximation  for  gamma  rays,  the  trans¬ 
port  of  which  is  dominated  by  absorption  and  tor 
ward  scattering. 

An  early  version  of  the  one-dimensional  uniform 
air  data  base  was  taken  from  the  Air  Transport 
of  Radiation  (ATR)  code  (Ref.  33).  These  data 
were  shown  to  be  in  error  when  compared  with 
those  generated  using  the  more  recent  Vitamin 
C  cross  sections  as  shown  in  figure  2i  .  The  cur¬ 
rent  data  base  has  been  developed  using  the  lat¬ 
est  cross  section  set.  Vitamin  E  (Ref.  60),  the 


results  from  which  agree  well  with  those  ca'cu- 
lated  using  Vitamin  C,  as  shown  in  Figure  28. 

A  single  one-dimensional,  uniform  air  gamma  ray 
transport  data  base  has  been  developed  to  serve 
prompt  and  delayed  applications.  Calculations  of 
gamma  ray  energy  and  angle-differential  fluence 
as  a  function  of  dislance  through  550  grams/cm2 
or  approximately  5000  meters  have  been  per¬ 
formed  for  sources  in  each  of  the  38  Vitamin  E 
energy  groups,  as  shown  in  Table  27.  The  calcu¬ 
lations  were  performed  using  a  S40  quadrature 
and  a  P3  Legendre  scattering  representation.  The 
results  of  the  calculations  have  been  collapsed 
into  the  18  gamma  ray  energy  group,  structure 
common  to  ATR5.  That  structure  is  also  shown 
in  Table  27.  The  collapsing  procedure  has  been 
performed  using  a  gamma  ray  source  spectrum 
inversely  proportions  to  energy.  The  effect  of  the 
ground  on  gamma  ray  transport  is  accounted  for 
using  the  correction  factor  approach  described 
in  tho  body  of  this  report. 

NUCLEAR  WEAPON  HYDRODYNAMICS. 

A  schematic  representation  of  nuclear  weapon 
hydrodynamics,  exclusive  of  fireball  rise,  is  given 
in  Figure  29.  In  an  air-ground  environment  the  ini¬ 
tial  air  shock  wave  created  by  the  rapid  expansion 
of  hot  gasses  around  the  weapon  reflects  off  the 
ground  and  rises  after  the  initial  wave,  forming 
the  Mach  stem,  in  which  the  initial  and  reflected 
shock  unite. 

Observations  of  field  tests  and  small  scale  high 
explosive  experiments  have  been  used  to  devel¬ 
op  large  empirical  and  semi-empirical  models  of 
weapon  hydrodynamics.  These  are  embodied  in 
computer  codes  which  are  by  and  large  too  cum¬ 
bersome  for  use  in  conjunction  with  transport 
data  for  delayed  gamma  ray  modeling.  However, 
small  codes  which  scale  data  generated  by  more 
complex  codes  have  been  developed.  STLAMB 
(Ref  9)  a  more  recent  variant  of  the  Low  Altitude 
Multi  Burst  (LAMB)  code  developed  by  Needham 
and  Wittwer  (Ref.  43)  is  one  of  these.  This  code 
interpolates  between  and  extrapolates  from  hy¬ 
drodynamic  data  obtained  from  model  calcula¬ 
tions  for  one  kiloton  and  one  megaton  yields  in 
uniform  air.  It  accounts  for  the  presence  of  the 
ground  and.  hence,  the  reflected  shock  by 
means  of  an  image  burst  below  the  ground  plane. 
STLAMB  Is  used  as  the  hydrodynamic  model  in 
the  delayed  gamma  radiation  system  used  to  cal¬ 
culate  survivor  air  density  profiles  as  a  function 
of  time  and  computes  the  amount  of  air  between 
source  and  detector.  It  also  calculates  cloud  rise. 


Figure  27.  Bren  Co-60  dose  at  1  i  detector  height. 
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Figure  28.  Uniform  air  tissue  kerma  from  a  Fisher  and  Engle 
U-235  gamma  ray  source  at  one  second. 


Table  27.  Vitamin  E  and  ATR5  gamma  ray  cross  section 
library  energy  group  structures. 


Upper  Energy 

(MeV) 


Vitamin  E 
Group  No. 


ATR5 

Group  No. 


Lower  Bound 


2.00+1* 

1 

1.40+1 

2 

1.20+1 

3 

18** 

1.00+1 

4 

17 

8.0C+0 

5 

16 

7.50+0 

6 

7.00+0 

7 

15 

6.50+0 

8 

6.00+0 

9 

14 

5.50+0 

10 

5.00+0 

11 

13 

4.50+0 

12 

4.00+0 

13 

12 

3.50+0 

14 

3.00+0 

15 

11 

2.50+0 

16 

10 

2.00+0 

17 

9 

1.66+0 

18 

8 

1.50+0 

19 

1.33+0 

20 

1.00+0 

21 

7 

8.00-1 

22 

7.00-1 

23 

6 

6.00-1 

24 

5.12-1 

25 

5.10-1 

26 

4.50-1 

27 

5 

4.00-1 

28 

3.00-1 

29 

4 

2.00-1 

30 

1.50-1 

31 

3 

1.00-1 

32 

2 

7.50-2 

33 

7.00-2 

34 

6.00-2 

35 

4.50-2 

36 

1 

3.00-2 

37 

2.00-2 

38 

1.00-2 

28 

♦Read  as  2.0*10l. 

**ATR  energies  are  numbered  from 
low  to  high  energy. 
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Figures  30  through  36  show  the  development  of 
the  system  of  hydrodyriamically  perturbed  air 
surrounding  a  burst  at  1 1 0  meters  scaled  height 
of  burst  (Burst  height  +  yield  1/3  )  as  modeled  us¬ 
ing  STLAMB.  Time-dependent  gamma  ray  KER- 
MA  measurements  were  made  on  an  equivalent 
test  (Ref.  11),  which  can  be  used  to  help  validate 
the  delayed  gamma  ray  model.  The  development 
of  the  weapon  hydrodynamic  process  includes  a 
spherical  density  well  which  is  later  perturbed  by 
the  reflected  shock  front.  By  1  second  the  re¬ 
flected  shock  is  shown  to  be  passing  through  the 
approximate  center  of  the  original  fireball  region. 
Thereafter,  the  shock  moves  rapidly  above  the 
now  rising  fireball,  which  eventually  forms  a  toroid 
shaped  system  and  moves  well  above  its  starting 
point. 

Figures  37  through  44  show  a  similar  progression 
for  a  burst  at  73  meters  scaled  height  of  burst, 
also  a  subject  of  gamma  ray  KERMA  rate  meas¬ 
urements,  alcng  with  that  of  a  very  similar  shot 
(71  meters  SHOB).  In  this  case  the  reflected 
shock  is  moving  back  through  the  fireball  less 
than  a  tenth  of  a  second  after  its  initiation.  By  a 
second  it  has  passed  entirely  through  the  fireball. 
It  is  shown  to  proceed  to  form  a  torus  similar  to 
that  of  110m  SHOB  shot. 

These  depictions  of  density  profiles  have  been 
created  with  data  from  STLAMB  and,  thus,  are 
the  basis  for  determining  tne  thickness  and  length 
of  the  transport  path  from  source  to  detector. 
STLAMB  also  computes  the  altitude  o<  tho  highest 
temperature  in  the  system,  which  may  be  asso¬ 
ciated  with  the  debris.  However,  the  debris  is 
probably  distributed  over  a  larger  volume  within 
the  fireball.  In  fact,  observations  of  field  tests 
suggest  that  the  debris  actually  flows  within  the 
fireball,  at  to  the  top,  'ater  down  the  sides  and, 
finally,  into  the  torus  and  skirt.  Unfortunately, 
there  are  not  sufficient  da*a  to  construct  a  model 
of  this  cloud  geometry-debris  location  system  on 
the  basis  or  actual  observations.  Thus,  the  pres¬ 
ent  delayed  radiation  model  has  been  developed 
under  the  assumption  that  the  debris  resides  at 
the  altitude  of  the  hot  spot  but  at  a  point  on  the 
centerline  of  the  system  rather  than  distributed 
as  a  sphere  or  torus.  This  assumption  has  been 
made  to  simplify  the  calculation  process  and  to 
keep  the  running  time  of  the  computer  model  at 
a  reasonable  level.  However,  the  effect  of  mov¬ 
ing  the  source  to  other  locations  has  been  ex¬ 
amined.  Calculations  have  been  made  for  debris 
locations  on  the  centerline  but  high,  0.8  of  the 


fireball  radius  above  the  fireball  center  and  low, 
0.8  of  the  fireball  radius  below  the  fireball  center. 
Calculations  using  these  alternative,  though  not 
necessarily  realistic,  debris  location  have  been 
made  for  comparison  with  experimental  mea¬ 
surements  and  to  determine  the  sensitivity  of 
delayed  gamma  ray  fields  to  such  perturbations. 


Calculations  of  time-varying  exposure  rates 
(roentgen  per  second)  as  measured  by  shielded 
ionization  chambers  (Conrad  devices)  have  been 
performed  using  the  revised  Fisher  and  Engle 
source  values,  transport  methods  as  previously 
described  and  STLAMB  hydrodynamic  modeling. 
The  ionization  chambers  were  located  near  the 
center  of  25  inch  lengths  of  steel  pipe  oriented 
vertically  above  the  ground.  The  outside  diameter 
of  the  pipe  was  5.5  inches  and  the  sides  were 
3/16  inch  thick.  Monte  Carlo  calculations  were 
performed  using  the  pipe-ground  geometry  sys¬ 
tem  to  obtain  the  adjoint  leakage  from  the  detec¬ 
tor,  which  could  be  convoluteu  with  the  free  field 
fluence  to  obtain  the  exposure  as  measured  by 
the  ionization  chambers.  The  energy  dependent 
response  of  the  ionization  chamber  was  taken  to 
be  that  described  by  Ehrlich  (Ref.  16)  and  de¬ 
picted  along  with  those  for  other  detector  types, 
relative  to  fluence-to-roentgen  conversion,  in 
Figure  45.  The  preferential  shielding  of  low  ener¬ 
gy  photons,  along  with  the  rapid  drop  of  detector 
response  at  low  photon  energies,  tends  to  re¬ 
move  the  excess  low  energy  fluence  component 
present  in  the  free  field  as  described  previously. 
Results  have  been  obtained  for  the  source  lo¬ 
cated  at  the  hot  spot  altitude  and  for  the  two  vari¬ 
ants  described  previously. 


Figures  46  through  52  present  exposure  rate 
measurement  and  calculation  results  for  the 
110m  SHOB  shot  for  horizontal  distances  from 
457  to  3261  km.  Figures  53  through  58  present 
results  for  the  73m  SHOB  shot  for  horizontal  dis¬ 
tances  from  457  to  2780  m.  Results  are  pres¬ 
ented  in  terms  of  exposure  rate  x  time.  This  is 
done  to  reduce  the  range  of  numerical  values  for 
plotting.  Data  corresponding  to  the  nominal  or  hot 
spot  location  is  depicted  using  a  heavy  solid  line. 
That  corresponding  to  the  low  source  location  is 
depicted  using  the  dotted  line,  that  for  the  high 
using  the  dashed  tine.  The  light  solid  line  depicts 
a  case  In  which  the  effect  of  the  ground  reflection 
of  the  shock  is  removed,  which  is  the  equivalent 
of  the  model  used  in  ATR  and  which  will  be  dis¬ 
cussed  tater. 
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Relative  Sensitivity 


Figure  45.  Highlights  of  energy  dependence  of  instrument 
sensitivity . 
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Source 


Figure  47.  Measured  and  calculated  expo?  rates  at  914 


-No  Reflected  Shock 


sure  rates  at  1371 


Source  location:  H-Hlgh,  N-Ncx»1n«l,  L-Low,  NRS-No  Reflected  Shock 
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T 1  HE  1SEC0N0S) 

Figure  49.  Measured  and  calculated  exposure  rates  at  1828  m,  110m  SHOD 


Source  Location:  N-Nominal  ,  l-Low,  NRS-No  Reflected  Shock 


Figure  50.  Measured  and  calculated  exposure  rates  at  2286  a,  110m  SH0I3 


N-Noflil nal  ,  l-tow,  NR 


Measured  aad  calculated  exposure  rates  at  326 


Source  Location:  H-Hlgh,  N-Nominal,  l-Low,  NRS-No  Reflected 


Measured  and  calculated  exposure  rates  at  457  m,  73m  S110B 


Nominal ,  l-Low,  NRS-No  Reflected  Shock 


Measured  and  calculated  exposure  rates  at  914  n,  73m  SHOB 


Source  Location:  H-Mgh,  N-Nomlnal,  L-Low,  NRS-No  Reflected  Shock 


Measured  and  calculated  ext*  mre  rates  at  137?  a,  73b  SHOB 


Nominal,  L-Low,  NRS-No  Reflected  Shock 


Source  Location:  H-Hlgh.  N-Nomlnal.  L-Low,  NRS-No  Reflected  Shock 


Source  Location:  H-Htgh,  N-Honlnal.  L-Low,  NRS-No  Reflected  Shock 


The  calculations  follow  the  general  shape  of  the 
measured  data  for  times  beyor.d  a  few  tenths  of 
a  second.  At  earlier  times  gamma  rays  from  air 
and  ground  capture  of  prompt  neutrons  dominate 
the  exposure  rate.  Those  are  not  part  of  the 
delayed  radiation  model.  At  times  between  a  few 
tenths  of  a  second  and  a  few  seconds  the  calcu¬ 
lations  display  the  effect  of  the  reflected  shock 
wave  passing  through  the  fireball.  The  measure¬ 
ments  made  on  the  1 1 0m  SH08  shot  also  show 
this,  while  the  effect  is  hidden  by  the  dominant 
secondary  gamma  rays  for  the  other  shot.  Be¬ 
yond  a  few  seconds  the  calculations  are  lower 
than  the  measurements,  quite  a  bit  lower  in  the 
case  of  the  110m  SHOB  shot. 

Before  examining  the  potential  effects  of  source 
movement  within  the  fireball  system  and  other 
factors  wnich  may  effect  the  calculation- 
measurement  comparison  it  is  useful  to  examine 
the  calculated  and  measured  exposure  rate  data 
in  time  intervals  which  allow  identification  of  data 
trends  in  time  and  space.  The  exposure  values 
calculated  in  four  intervals  using  the  nominal  (hot 
spot)  source  locations  are  presented  as  ratios  to 
the  equivalent  measured  values  in  Table  28. 
There  are  several  types  of  errors  to  look  for  from 
the  evidence  of  these  ratios.  A  source  magnitude 
error  should  cause  a  uniform  discrepancy  at  all 
ranges,  as  should  an  error  in  calculating  the  size 
of  the  density  well.  An  error  in  fireball  rise  should 
cause  a  discrepancy  which  is  range-dependent, 
as  should  an  error  in  specifying  the  source  loca¬ 
tion  relative  to  the  fireball. 

During  the  earliest  period  the  calculation- 
measurement  ratios  for  the  1 10m  SHOB  shot  ex¬ 
hibit  a  range  dependence.  Either  a  faster  fireball 
rise  or  a  rise  of  the  debris  within  the  fireball  would 
reduce  the  discrepancy.  However,  photographic 
evidence  (Ref.  12)  indicates  that  for  a  similar 
shot,  Upshot-Knothol6  CLIMAX,  the  debris  re¬ 
main  in  the  fireball  center  through  approximately 
one  second.  This  leaves  the  tireball  rise  rate  as 
suspect,  particularly  for  large  yields,  a  suspicion 
which  is  supported  by  the  same  photographic  evi¬ 
dence.  During  the  period  from  one  to  three  sec¬ 
onds  the  calculation-measurement  ratios  for  Shot 
the  110m  SHOB  shot  still  exhibit  a  discrepancy 
at  close  range  through  to  a  lesser  extent  than  dur¬ 
ing  the  prior  period.  The  higher  tireball  caused 
by  an  increased  rise  rate  during  the  period  to  one 
second  would  probably  account  for  most  of  the 
discrepancy,  in  addition  photographic  evidence 
suggests  that  the  source  is  rising  relative  to  the 
center  of  the  fireball. 


In  the  case  of  the  73m  and  71  m  SHOB  shots  the 
calculation-measurement  ratios  during  the  first 
two  periods  show  good  agreement  in  the  first  two 
periods,  the  mean  ratio  for  the  first  period  being 
1.02  ±  .12  and  that  for  the  second  being  .88 
±  .10,  where  the  means  are  based  on  the  ratios 
for  both  shots  and  the  uncertainty  is  the  standard 
error  of  the  distribution. 

During  the  periods  beyond  three  seconds  the 
calculation-measurement  ratios  for  all  three  ex¬ 
hibit  systematic  variations  with  distance,  being 
low  close  in  and  high  far  out.  According  to  the 
Figures  46  through  58  those  discrepancies  could 
be  eliminated  removing  the  provision  for 
ground-shock  reflection  from  the  model.  The  re¬ 
sults  of  such  a  modification  are  depicted  as  a  light 
solid  line  in  the  figures.  However,  the  early  time 
structure  observed  in  the  measurements  is  no 
longer  observed  in  the  calculations  after  eliminat¬ 
ing  the  reflected  shock.  Also,  the  resulting  fireball 
rise  rate  is  unrealistically  slow,  as  compared  to 
test  shot  photographs.  However,  this  contributes 
to  the  improved  ability  of  the  unreflected  model 
to  calculate  dose  rates  at  iate  times.  The  hydro- 
dynamic  model  used  in  ATR  is  essentially  equiva¬ 
lent  to  the  non-reflected  shock  variant  of  the 
LAMB  model  described  above.  Like  that  model. 
ATR  over-predicts  the  dose  rate  In  the  important 
period  over  the  first  few  seconds  after  the  re¬ 
flected  shock  wave  has  passed  through  the  fire¬ 
ball.  However,  unlike  the  LAMB  model  the  ATR 
cloud  rise  model  is  independent  of  the  wock 
propagation  process,  being  based  on  atmo¬ 
spheric  test  observations. 

The  comparison  between  (STLAMB)  calculated 
and  measured  exposures  in  the  three  to  twenty 
second  time  periods  can  be  improved  by  lowering 
the  source  relative  to  the  fireball  center.  Photo¬ 
graphic  evidence  suggeststhat  such  as  phenom¬ 
enon  does  occur  after  the  debris  rises  to  the  top 
of  the  fireball,  which  takes  place  at  approximately 
three  seconds.  Table  29  provides  an  estimate  of 
the  effect  of  lowering  the  source  (but  not  the  fire¬ 
ball)  .  Those  data  are  based  on  those  presented 
in  Figures  46  through  58.  As  shown  in  the  Table. 
H  stands  for  high  source,  N  for  nominal  and  L  for 
low.  Intermediate  locations  are  given  as  NG  and 
NL  and  are  simple  averages,  representing  a  mid¬ 
point  between  the  two  cases  for  which  calcula¬ 
tions  have  been  performed.  Note  that  raising  the 
source  for  Shot  Hood  between  one  and  three  se¬ 
conds  reduces  but  does  not  eliminate  the  dis¬ 
crepancy.  indicating  that  a  change  in  cloud 
location  is  probably  necessary  also. 
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Table  29.  Delayed  pww  ray  aodel  exposure  calculation  to 
aeasureaent  ratios  (aoving  source  location)  for 
three  JfTS  shots. 


Shot 

Tlae 

Interval 

Source 

Location 

Ranqe 

(») 

457 

914 

1372 

1289 

2286 

2779 

3261 

SHOB:  110a 

C  -  1 

n 

1.08 

1.34 

1.26 

1.16 

0.97 

0.93 

1.05 

1  -  3 

m 

1.40 

1.13 

0.95 

0.93 

0.78 

0.84 

0.84 

3  -  10 

HI 

* 

0.95 

0.84 

* 

0.85 

0.89 

0.89 

1C  -  20 

L 

• 

1.10 

0.90 

* 

0.76 

0.85 

0.75 

0  •  20 

• 

1.15 

0.99 

• 

0.85 

0.88 

0.90 

SHOB:  73a 

0  -  1 

H 

1.18 

:  .06 

0.89 

1.00 

1.02 

0.91 

# 

1  -  3 

H 

0.82 

0.88 

0.88 

0.97 

0.88 

0.95 

• 

3  -  10 

Nl 

0.(7 

0.74 

0.72 

0.82 

0.75 

0.79 

• 

10  -  20 

L 

1.30 

1.18 

0.95 

0.95 

0.75 

0.68 

• 

0  -  20 

0.91 

0.90 

0.84 

0.92 

0.84 

0.84 

• 

SH03:  71a 

0  -  1 

* 

0.81 

1.11 

1.20 

1.01 

1.06 

• 

• 

1  -  3 

n 

0.62 

0.91 

0.99 

0.88 

0.94 

• 

# 

3  -  10 

m. 

0.62 

0.65 

0.72 

0.67 

0.73 

# 

• 

10  -  20 

L 

1.31 

0.76 

0.71 

0.67 

0.63 

• 

• 

o 

• 

o 

0.70 

0.85 

0.89 

0.80 

0.83 

• 

• 
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Lowering  the  sources  at  iate  times  improves  the 
agreement  between  measurement  and  calcula¬ 
tion  but,  just  as  importantly,  it  reverses  the  nature 
of  the  discrepancy  as  it  vanes  with  range  The 
remaining  disc'epancy  ca~  be  accounted  for  to 
a  large  extent  on  the  bosis  of  the  racial  oisc'ace- 
ment  of  the  source  as  t  moves  down  *Ke  ou's  ce 
pi  the  fireball  and  eventua  l,  becomes  net.  un¬ 
rated  m  the  torus 

At  times  between  five  and  ten  seconds  the  outer 
surface  of  the  rising  fireball  begins  to  coot,  the 
vapors  are  subiect  to  drag  from  the  surrounding 
atmosphere.  The  resulting  downward  flow  around 
the  outside  of  the  fireball  is  matched  by  an  up¬ 


ward  flow  at  its  center.  This  flow  pattern  produces 
a  torus,  within  which  reside  the  debris.  The  radius 
of  the  torus  depends  on  yield.  Thus,  the  torus  for 
the  110m  SHOB  shot  attains  a  radius  approxi¬ 
mately  325  meters  by  ten  seconds,  while  those 
for  71  and  73m  SHCB  shots  are  aoproximately 
200  me’ers  m  radius  T"e  low  about  the  torus 
causes  rac-d  m-«,ng.  c.r.c  •eg  and  return  tc  the 
near  ambient  density  U'  der  such  condit  ons  the 
radial  displacement  of  the  debris  causes  an  in¬ 
crease  in  exposure  which  increases  with  range, 
in  the  case  of  the  1 10m  ShOB  shot  that  increase 
m  terms  of  multiplicative  factors  has  been  esti¬ 
mated  usmg  simple  pomt-kerna!  calculation  tech¬ 
niques.  as  follows: 


Radial  Source  Displacement  Correction  Factor.  110  SHCB 


Time  Period  Range 
( sec )  i  m ) 

914 

1372 

2236 

2779 

3261 

3-10 

1  04 

1.07 

115 

1  19 

1.23 

10-20 

i  13 

1  17 

1  25 

1  30 

l  35 

These  factors  have  been  applied  to  the  110m 
SHOB  shot  source  height-adjusted  calculation- 
measurement  ratios  and  the  results  provided  in 
Table  30.  The  data  in  Table  30  represent  a  defi¬ 
nite  improvement  over  the  original  calculation- 
measurement  ratios,  incorporating  nominal 
source  location  i  nta  only,  as  found  ;n  Table  28. 

OTHER  MODEL  COMPARISONS  WITH 
WEAPON  TESTS. 

There  are  other  comparisons  which  can  be  made 
with  data  from  weapon  tests  in  addition  to  those 
discussed  previously-  jwever,  most  involve  to¬ 
tal  exposure  rather  than  exposure  rate. 

Total  gamma  ray  exposure  data  are  available  for 
device  F.  22  Kt  yield,  air  dropped  and  eroded 
at  a  height  of  437m.  Data  are  also  available  for 
device  0,  21  Kt  yield,  ah’  dropped  and  exploded 
at  a  height  of  432m.  Both  devices  were  similar 
to  Fat  Man  in  design  and  high  explosive  thickness. 


Calculated  and  measured  (Refs.  38  and  37)  gam¬ 
ma  ray  exposure  intensities  are  provided  for  de¬ 
vice  F  in  Figures  59  and  80  and  for  device  D  in 
Figures  61  and  62.  Each  shot  had  two  orthogonal 
lines  ol  film  dosimeters,  one  running  west  from 
the  intended  ground  zero  and  one  running  south. 
The  measured  values  have  been  revised  as  rec¬ 
ommenced  by  EHery  Storm.  LANL  (Ref.  52). 
Those  revisions  include  corrections  tor  betatron 


calibration  energy  (meas.  *  .885)  and  for 
energy-dependence  of  film  sensitivity  (sensitive 
films  (exposure  iess  than  1 0r)  *  .909.  insensitive 
films  (exposure  greater  than  tOr)  *  .840). 

The  calculated  exposure  values  are  generally  10 
to  20%  below  the  revised  measured  values.  How¬ 
ever.  agreement  is  somewhat  better  in  the  first 
2000  meters  than  at  greater  ranges.  Taking  the 
source  location  into  account  is  likely  to  improve 
the  agreement,  particularly  at  long  ranges.  How¬ 
ever.  for  shorter  ranges  the  effect  is  likely  to  be 
minimal  given  the  high  burst  height. 

In  summary,  comparisons  of  results  of  the  bn?t 
available  model  for  delayed  gamma  ray  calcula¬ 
tion  with  results  of  gamtna  ray  exposure  rate 
measurements  made  at  atmospheric  tests  indi 
cate  that  the  model  has  short  comings  which  are 
mai.ify  due  to  incorrect  hydrodynamic  data  or  in¬ 
adequate  modeling  detail  based  on  existing  data. 

The  hydrodynamic  model  used  in  tne  ATM  code 
suffers  from  deficiencies  similar  to  those  attrib¬ 
uted  to  the  more  sophisticated  model.  It  does  not 
account  for  source  motion  relative  to  that  of  the 
fireball,  nor  does  it  account  for  the  effect  of  tor  js 
formation.  However,  as  described  in  the  body  of 
this  report,  the  properties  of  the  cloud  rise  and 
blast  model  incorporated  in  ATR  do  allow  succes¬ 
sful  scaling  based  on  scaled  burst  height  to  obtain 
a  reasonable  value  for  the  time-integral  delayed 
gamma  ray  dose,  although  the  rate  data  may  be 
Significantly  in  error  at  any  giver  time. 


Table  30.  Delayed  gamma  ray  todel  exposure  calculation  to 
measurement  ratios  for  the  110m  SHGB  shot  (moving 
source  with  adjustment  for  Torus  formation. 


Time  Source 


Interval 

Location 

Range:  457 

(m) 

914 

1372 

1829 

2286 

2779 

3261 

0  -  1 

N 

1.08 

1.34 

1.26 

1.16 

0.97 

0.93 

1.05 

1  -  3 

NH 

1.40 

1.13 

0.95 

0.93 

0.78 

0.84 

0.84 

3  -  10 

m. 

* 

0.99 

0.90 

* 

0.98 

1.06 

1.09 

10  -  20 

t. 

* 

1.24 

1.05 

• 

0.95 

1.11 

1.01 

O 

CVJ 

1 

o 

* 

1.16 

1.02 

* 

.91 

.97 

1.00 
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